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I. INTRODUCTION 

Current theoretical cosmology is faced with two fundamental problems: what is inflation and what is dark energy? 
In other words, why and how did both the very early and the very late universe expand with acceleration? What is 
the reason for the similarity of the universe evolution at small as well as at large curvature? 

The scenarios to describe these early-time and late-time accelerations are usually very similar, which is why it is 
quite natural to expect that the same theory /principle lies behind both accelerating cosmological epochs. Indeed, 
there are various proposals to construct an acceptable dark energy model, such as: scalar, spinor, (non-)abelian vector 
theory, cosmological constant, fluid with a complicated equation of state, and higher dimensions. Remarkably, the 
same proposals are also intended to describe inflationary era. As a result, we have a number of competing scenarios 
with which to describe inflation and late-time acceleration. This situation is very reasonable due to the fact that the 
evolution of the cosmological parameters is not defined with precise accuracy. Even their current values are defined 
with at least 3-5% error. Moreover, working in framework of one of the above proposals, we are forced to introduce 
the following extra cosmological components: inflaton, a dark component and dark matter. Even if such a scenario 
seems to be partially successful, it immediately introduce a new set of problems: such as the following: coupling with 
usual matter, (anti-)screening of dark components during the evolution of the universe, compatibility with standard 
elementary particle theories, and consistency of formulation. Thus, another natural proposal is to not introduce extra 
fields to resolve the cosmological problems, which is the most economical solution in the spirit of the Occam's razor 
principle. 

The least understood and most fundamental force of nature is the gravitational interaction. It is expected that 
general relativity is just an approximation that is valid at small curvatures and that is predicted by the Theory of 
Everything (strings/M-theory). In the very early universe, some unknown gravitational theory dictated the evolution 
of the universe. Thus, a very reasonable assumption is the gravitational alternative for a unified description of the 
inflation and dark energy. In other words, the modified gravity, which represents a classical generalization of general 
relativity, should consistently describe the early-time inflation and late-time acceleration, without the introduction of 
any other dark component. A sector of modified gravity that contains the gravitational terms relevant at very high 
energies produced the inflationary epoch. In the course of evolution, the curvature decreases and general relativity 
gives a sufficient approximation in the intermediate universe. With a further decrease of the curvature, when sub- 
dominant terms quickly grow, one sees the change from deceleration to cosmic acceleration. Thus, the early-time as 
well as the late-time cosmic speed-up is caused simply by the fact that some sub-dominant terms of gravitational 
action may become essential at large or small curvatures. Moreover, such an approach may also be considered as 
dynamic solution of the cosmological constant problem. Of course, the complete gravitational action should be defined 
by a fundamental theory, which remains to be the open problem of modern high-energy physics. In the absence of 
fundamental quantum gravity, the modified gravity approach is a phenomenological model that is constructed by 
complying with observational data and data from local tests. 

In addition, just as by-product to be a cosmological model, the modified gravity may provide the explanation for 
dark matter. It may resolve the coincidence problem simply by the fact of the universe expansion. It also may describe 
the transition from deceleration to acceleration of the universe, and it may be useful for high-energy physics problems 
(i.e., unification of all interactions, and hierarchy problem resolution). Even if the current universe is entering the 
phantom phase, modified gravity effectively describes the transition from the non-phantom to phantom era without 
the need to introduce the exotic matter (phantom) with extremely strange properties. 

In this work, we present a detailed review of a number of popular models of modified gravity. Their properties 
and different representations are discussed. Assuming the spatially flat FRW cosmology, we investigate the modified 
gravity background evolution with the aim of achieving a unified description of the universe's inflation with dark 
energy epoch. Using sufficient freedom in the choice of the action as some function of curvature invariants, we develop 
the cosmological reconstruction scheme, which is explicitly applied to obtain the (effective quintessence/cosmological 
constant/phantom) accelerating expansion for the models under discussion. Late-time dynamics is analyzed and it is 
shown that the future universe may end in one of the four known finite-time singularities. For instance, the effective 
phantom superacceleration brings the universe to the famous Big Rip (or Type I) singularity [lj. In order to avoid 
the future singularity one has to add the extra higher-derivative gravitational terms, which are relevant at the early 
universe and help in the unification of early-time inflation with late-time acceleration because such terms induce 
inflation. 

It should be noted that, as a rule, modified gravity equations of motion are higher-derivative differential equations. 
Thus, we consider only the background evolution of such theories in this report; we do not discuss the cosmological 
perturbation theory. The reason for this is very simple. Actually, in almost all existing approaches to cosmological per- 
turbations in modified gravity, the corresponding equations are reduced to second-order differential equations. Clearly, 
this is not a well-justified approach, and the results obtained in this way are questionable. The new higher-derivative 
covariant perturbation theory should be developed to obtain correct answers to the related circle of questions. 
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The work is organized as follows. The second chapter gives the introduction to several popular models of modified 
gravity: i 7 '(i?)-theory, modified Gauss-Bonnet and string- inspired gravity, non-local gravity, non-minimally coupled 
theories, F(R) Hofava-Lifshitz gravity, and power-counting renormalizable covariant gravity. The general properties 
of -F(i?)-theory are discussed in section III Al Its different representations, namely Jordan frame, scalar-tensor form 
(Einstein frame) and fluid representation, are given. The thermodynamics of Schwarzschild-(anti-)de Sitter space 
are briefly discussed. Some local tests, in particular, the Newtonian regime and matter instability, are investigated. 
Several specific models (i.e., theory with negative and positive curvature power terms, lni?-model, and exponential 
theory), are presented. The spatially flat FRW cosmology is investigated. Viable modified gravities that allow the 
unification of the early-time inflation with late-time acceleration are given (for the first work where such a unification 
was proposed, see 0). Note that the literature on F(i?)-gravity is vast, and the earlier reviews on this theory are 
given in refs. @, H. Different aspects of F(i?)-gravity are discussed in refs. 0-0], and various exact solutions are 
presented in refs. Q. The comparison with observational data and local tests is made in refs. [1, 0,0- As is indicated 
in this work, there are a number of viable F(R) models that successfully pass local tests and are in accordance with 
observational data. 

Modified Gauss-Bonnet gravity is considered in section IIIB| and its spatially flat FRW equations are given. The 
de Sitter solution and its thermodynamics are also studied, and the possibility of the unification of inflation with 
dark energy is mentioned. Section III CI is devoted to a brief study of string-inspired (scalar-Gauss-Bonnet) gravity. 
The emergence of an accelerating dark energy epoch in such a model is proposed. The non-local generalization of 
F(i?)-gravity and Gauss-Bonnet gravity is discussed in section III Dl We find the scalar-tensor representation for 
such non-local gravities and investigate their spatially flat FRW cosmology. The possibility of early-time or late-time 
acceleration is again demonstrated. Non-minimally coupled theories are presented in section III El We consider the 
theories where scalar, vector or Yang-Mills Lagrangian is multiplied to function F(R). FRW equations are given and 
possibility of cosmic acceleration is mentioned. In section III Fl we propose modified F(R) Hofava-Lifshitz gravity. The 
original Hofava-Lifshitz theory, based on general relativity, has been introduced as a candidate for renormalizable 
quantum gravity at the price of the violation of the Lorentz invariance. The F(R) generalization of Hofava-Lifshitz 
gravity yields consistent FRW equations that coincide with equations associated with the conventional theory by 
choosing the parameters to be specific values. We study several explicit models that are introduced in section III Al 
in this framework and show that they may again predict the unified description of early-time inflation with dark 
energy. Motivated by the Hofava-Lifshitz proposal, the covariant power-counting renormalizable gravity is introduced 
in section III Gl Two formulations of the theory: one utilizing an extra perfect fluid and the other using the Lagrange 
multiplier constraint, are given. Accelerating FRW cosmology is briefly discussed. In section Hi Dl it is shown that 
arbitrary modified gravity may be interpreted as an effective fluid with an inhomogeneous equation of state. The 
late-time evolution of a specific dark fluid with a power-law equation of state is investigated. It is demonstrated that 
an effective quintessence/phantom fluid may drive the universe to a finite-time future singularity. 

The background evolution of modified gravity is studied in the third chapter. Usually, in theories like general 
relativity, we start from a theory that is defined by the action and solve the corresponding equations of motion to define 
the background dynamics. However, an approach is already too complicated for simple models of modified gravity. 
Thus, using the fact that modified gravity is defined up to some arbitrary function of curvature invariants, the inverse 
problem is analyzed. The corresponding cosmological reconstruction scheme is proposed. Within the framework of 
this scheme, we show how the complicated background cosmology, which complies with observational data, may be 
reconstructed. This provides way to define (at least partially) the reconstructed form of modified gravity, which has 
a particular cosmological solution. Section IIII Al is devoted to a formulation of cosmological reconstruction (in terms 
of cosmological time or in terms of e- folding N) for scalar-tensor gravity. Several examples of a (super) accelerating 
epoch are used for the reconstruction of scalar potentials. This formulation is extended for the theory with two scalars 
where the reconstruction of the transition from matter dominance to late-time acceleration as well as a reconstruction 
of the ACDM epoch is done. In section IIIIBI we study the reconstruction in the k-essence model, using both 
formulations. Cosmological reconstruction for the F(i?)-model, using its presentation with an auxiliary scalar, is 
developed in section HlI CI The differential equation that provides the form of function F(R) for a given spatially flat 
FRW cosmology is obtained. The ACDM epoch and dark era with a phantom divide crossing are reconstructed in 
F(i?)-gravity with/without matter in great detail. The general condition for the stability of the cosmological solution 
is derived. i ? (i?)-gravity with a Lagrange multiplier constraint is formulated, and its FRW cosmology is studied. In 
sections IIII Dl and IIII El we propose the reconstruction of modified Gauss-Bonnet and scalar-Einstein-Gauss-Bonnet 
gravities, respectively. In section IIII Ft we show that the reconstruction is easily generalized for the case of F(R) 
Hofava-Lifshitz gravity. The explicit reconstruction of the ACDM epoch and of phantom superacceleration is fulfilled. 
The cosmological reconstruction of the power-law FRW solution for non-minimal Yang-Mills theory is given in section 
InTGl 

The fourth chapter is devoted to an investigation of late-time dynamics for effective quintessence/phantom dark 
energy. It is demonstrated that, regardless of whether the dark energy model utilizes fluid, scalar or alternative 
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gravity, the future universe may evolve to a finite-time singularity of one of four known types. The prescription to 
resolve such a future singularity via the addition of a higher-derivative gravitational term is proposed. In section HV Al 
we consider dark fluid coupled with dark matter. For a power-law equation of state, it is shown that the future dark 
era may be singular, depending on the parameter choice. The coupling with dark matter may prevent some of the 
soft future singularities. However, such coupling does not resolve the future singularity in a general case. Only the 
addition of the i? 2 -term may cure all types of future singularities in the model under consideration. In section IIVB[ 
using the reconstruction method, we show that F(R) gravity may have late-time accelerating cosmological solutions 
with all four types of finite-time future singularities. To avoid this, one should additionally modify the theory by 
the R 2 term. It is interesting that the use of such a term may be motivated by conformal anomaly considerations. 
The same problem with qualitatively similar results is studied for modified Gauss-Bonnet and F(R) Hofava-Lifshitz 
gravities in section IIV CI and IIVDI It is remarkable that, in the last case, one may need other higher-derivative 
gravitational terms to resolve all types of future singularities. 
A summary and outlook are given in the Discussion section. 



II. MODIFIED GRAVITY UNIFYING THE EARLY-TIME INFLATION WITH LATE-TIME 

ACCELERATION 

This chapter is devoted to the general review of many popular models of modified gravity such as F(R) gravity, 
modified Gauss-Bonnet gravity, string-inspired theory, non-local gravity and non-minimal theories. Their general 
properties, different representations and background spatially flat FRW cosmologies are investigated. Special attention 
is paid to possible unification of early-time inflation with late-time acceleration within such theories. The extension 
of such a study to power-counting renormalizable gravity with explicit or apparent Lorentz symmetry breaking is also 
made. 



A. F(R) gravity 

In this section, we give a review of general properties of F(R) gravity. Its scalar-tensor description is presented, and 
some constant curvature solutions are discussed. A number of popular models are introduced, and their applications 
for a unified description of inflation with dark energy are studied. 



1. General properties 

Let us start with the general introduction of the model. In F(R) gravity Q, the scalar curvature R in the Einstein- 
Hilbert action 1 



SeH = j d^Xy/^g ^— + £mattcrj , (II. 1) 

is replaced by an appropriate function of the scalar curvature: 

S F (R) = J d^x^—g + £ roatter J . (II.2) 

Let us review the general properties of F(R) gravity. For F(R) theory, as for any other modified gravity theory, 
one can define an effective EoS parameter using its fluid representation. The FRW equations in the Einstein gravity 
coupled with perfect fluid are: 



3 1 

Pmatter = ' Pmattor = + 2Zf J . (H-3) 



We use the following convention for the curvatures: 

= ' Rpv , Rfiv = R . -R A MP i/ = — Fflp,v + F/jjap ~~ rJJ p r^ + V^vTpy , = — (p" (cjp„,A + 9Ai/,p — 9ft\,v) ■ 
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For modified gravity, one may define an effective equation of state (EoS) parameter as foliows: 

w ^-- 1 ~i§2- ( IL4 ) 

This relation is very useful for a generalized fluid description of modified gravity. 
The equation of motion for modified gravity is given by 

1 K 2 

-g^F(R) ~ R^F'(R) - g^UF'{R) + V„V V F'(R) = -—T mattcI ^ . (II.5) 
Assuming a spatially flat FRW universe, 

ds 2 = -dt 2 + a{t) 2 {dx 1 ) 2 , ( IL6 ) 

1=1,2,3 

the FRW like equations are as follows: 
F(R) 



= — 



3 (H 2 + ff) F' (R) - 18 (4H 2 H + HH} F" (R) + K 2 /Wter , (II.7) 



F(R) 



= — [H + 3H 2 ) F 



?'(R) + 6 (8H 2 H + AH 2 + 6HH + 'h) F"{R) + 36 (iHH + H^j 2 F"'(R) 



+ K Pmattor ■ (H.8) 

Here, the Hubble rate H is defined by H = a/ a and the scalar curvature R is given by R = 12H 2 + 6H. 

Several (often exact) solutions of pi.7[> maybe found. However, due to the complicated nature of field equations, 
the number of exact solutions is much less than, for instance, that in general relativity. Without any matter, assuming 
that the Ricci tensor is covariantly constant, that is, i? M „ cx g^, Eq. (III. 51) reduces to the algebraic equation: 

= 2F(R) - RF'(R) . (11.9) 

If Eq. pi.9p has a solution, the (anti-)de Sitter and/or Schwarzschild- (anti-)de Sitter space 

2 / 2MG r 2 \ , / 2MG r 2 \ _1 2 , 2 , 
ds 2 = -[l-—- T -yt 2 +[l-—- T -} dr 2 +r 2 dn 2 , (11.10) 

or Kerr - (anti-)de Sitter space is an exact vacuum solution. In pi.lO[) . the minus and plus signs in ± correspond to 

2 

the de Sitter and anti-de Sitter space, respectively. M is the mass of the black hole, G = and L is the length 
parameter of (anti-)de Sitter space, which is related to the curvature R = ±j§ (the plus sign for the de Sitter space 
and the minus sign for the anti-de Sitter space). For example, if we consider the model 

,,2m+2 

HR) = R-^ W r: (n.ii) 

with a constant /i with a mass dimension, the solution of (|II.9|) is given by 

R= (2 + m)^ rr /i 2 . (11.12) 

As another example, for a logarithmic model 

F(i?)= M 2 lnA (n.13) 

Ko 

with a constant i?o> one gets 

R = R e? . (11.14) 

The entropy of the Schwarzschild - (anti-)de Sitter space can be evaluated by the WKB approximation for the 
Euclidean partition function Z [lOj : 

Z ~ e SsdS . (11.15) 
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Here, 5sds is the value of the classical action calculated in the Schwarzschild - (anti-)de Sitter space. Because the 
value of 5*sds diverges, we regularize the action by subtracting the value when the solution of the (anti-)de Sitter 
space without a black hole is substituted [To| . The free energy T is given by 

T=~\a.Z. (11.16) 



And, using the standard thermodynamics relations, 



one can evaluate the entropy S. Here, Th is the Hawking temperature. Then, it follows that 



* = (".IS) 



Here, Ah expresses the area of the horizon. 

Let us now assume that F(R) behaves as F(R) oc foR m - Eq. ()II.7|) gives 

= f |-i (6H + 12H 2 Y + 3m (h + H 2 ^j (qH+WH 2 }™ * - 3mff^ | (6H + UH 2 }™ * j j + K 2 p Q a,-* {1+w) . 

(11.19) 

Eq. (111.81) is irrelevant because it can be derived from 19[) . When the contribution from the matter can be neglected 
(po = 0), the following solution appears: 



(m-l)(2m-l) 

H a-* , (H.20) 



with a corresponding EoS parameter (|II.4[) : 



6m 2 — 7m — 1 . „_ . 

WcS = ~ T( "TT^ -T • H.21 

3(m — l)(2m — 1) 

On the other hand, when the matter contribution is include with a constant EoS parameter w, an exact solution of 
(III.19P is given by 

h 2m 
a = atf , ho — 



3(1 + w) 



«o 



and the effective EoS parameter (|II.4[) is 



5^ (-6/io + I2h 2 )" 1 1 {(1 - 2m) (1 - to) - (2 - m)h } 



(11.22) 



w eS = -l + ^^. (11.23) 

TO 

These solutions (|II.20I) and (|II.22[) show that modified gravity may describe early/late-time universe acceleration. 
Moreover, it is very natural to propose that a more complicated modified gravity from the above class may give the 
unified description for inflation with late-time acceleration. 

2. Scalar-tensor description 

Let us now introduce perfect fluid and scalar-tensor representations for the modified gravity under consideration. 
These mathematically equivalent representations may simplify the investigation of a number of problems. 

Eqs. pi.7l) and ([II. 8[) show that the effective energy density p e g and the effective pressure p e g including the contri- 
bution from f(R) are given by (see, for instance, [111 ]) 



1 / 1 

Peff 



h— 



f(R) + 3 (H 2 + H) f'(R) - 18 (AH 2 H + HH^j f"(R)\ + p mattcr , (11.24) 



Peff = 4 (lf(R) - ( + H) f'(R) + 6 (8H 2 H + AH 2 + 6HH + h) f"(R) + 36 UhH + H) ' f"{R) 



k 2 \2' 

+Pmattcr • (11.25) 
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The equations pi.7[) and pi.8[) can be rewritten as in the Einstein gravity case (|II.3[) 

3 



„ = ~H 2 , p e£ f = 2 ( 3i?2 + 2 #) ' ( IL26 ) 



K 2 ' 



Such a mathematically equivalent fluid representation for the FRW equations in modified gravity may lead to a 
number of assumptions that are not well justified. Indeed, the generalized gravitational fluid contains higher-derivative 
curvature invariants. This fact is often ignored in the study of cosmological perturbations in modified gravity when such 
a theory is developed in strict similarity with general relativity, such that corresponding differential equations become 
of second order (see, for instance, [IH). Indeed, first of all, the field equations of the theory under investigation are of 
fourth order. It is evident that cosmological perturbation equations should be of fourth order as well. Moreover, such 
equations should be covariant. The corresponding covariant and higher-derivative cosmological perturbation theory 
(for review, see [13, HH) is developed, but it is extremely complicated. 

One can also rewrite F(R) gravity in the scalar-tensor form. By introducing the auxiliary field A, the action (|II.2|) 
of the F(R) gravity is rewritten in the following form: 

S = / d 4 xy^{F'(A) (R-A)+ F(A)} . (11.27) 



2k 2 d 

By the variation of A, one obtains A — R. Substituting A — R into the action (|II.27[) . one can reproduce the action 
in (|II.2|) . Furthermore, we rescale the metric in the following way (conformal transformation): 

9w -> e"ff M „ , a = - In F'(A) . (11.28) 
Thus, the Einstein frame action is obtained: 



S E = ^2 / d A x^—g (r - \gt>°d p ad a o - V(a 



V(a) = c°g (a"') - e 2 <7 (g ( e -)) = _ . (II .29) 

Here g (e _cr ) is given by solving the equation a — — In (1 + f'(A)) = — InF'(A) as A = g (e _<T ). Due to the scale 
transformation pi.28|) . a coupling of the scalar field a with usual matter arises. . The mass of a is given by 

2 _3^)_3f A 4F(A) 1 1 

a ~ 2 da 2 "2 \F'(A) (F'(A)) 2 F"(A) J ' 1 ' 

Unless m a is very large, the large correction to Newton's law appears. Naively, one expects the order of the mass m a 
to be that of the Hubble rate, that is, m a ~ H ~ I0 _33 eV, which is very light and could make the correction very 
large. 

In [lq . a "realistic" F(R) model was proposed. It has been found, however, that the model has an instability 
where the large curvature can be easily produced (manifestation of a possible future singularity). In the model of 
[l5j . a parameter m ~ 10~ 33 eV with a mass dimension is included. The parameter m plays a role of the effective 
cosmological constant. When the curvature R is large enough compared with m 2 , R ^> m 2 , F(R) [l^] looks as follows: 

F(R) = R-c 1 m 2 + ^^ + 0(R- 2n ) . (11.31) 
Here ci, C2, and n are positive dimensionless constants. The potential V(a) (jll.29j) has the following asymptotic form: 

V{*)~™£-. (11.32) 

Then, the infinite curvature R = A — > oo corresponds to a small value of the potential and, therefore, the large 
curvature can be easily produced. 

Let us assume that when R is large, F(R) behaves as 

F(R) - F R e . (11.33) 

Here, Fo and e are positive constants. We also assume e > 1 so that this term dominates compared with general 
relativity. Then, the potential y(cr) pi.29[l behaves as 
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Therefore, if 1 < e < 2, the potential V(a) diverges when R — > oo and, therefore, the large curvature is not realized so 
easily. When e = 2, V(a) takes a finite value 1/Fq when R — > oo. As long as 1/Fq is large enough, the large curvature 
can be prevented. 

Note that the anti-gravity regime appears when F'(R) is negative, which follows from Eq. (|II.27[) [2|. Then, we 
need to require 



We should also note that 



Therefore, if 



F'(R)>0. (11.35) 

dV(a) _ F"(A) 



dA F'{A)- 



(-AF'(A) +2F(A)) . (11.36) 



= -AF'(A) + 2F(A) , (11.37) 

the scalar field er is on the local maximum or local minimum of the potential and, therefore, a can be a constant. 
Note that the condition pi. 371) is nothing but the condition pi.9[) for the existence of the de Sitter solution. When 
the condition pi.37[) is satisfied, the mass pi. 301) can be rewritten as 



2F'(A) V F"(A) 

Then, when the condition (|II.35|) for the non-existence of the anti-gravity is satisfied, the mass squared m 2 is positive 
and, therefore the, scalar field is on the local minimum if 

- A + ¥$j >0 - (IL39) 

On the other hand, if 

~ A+ I^y <0 ' (IL40) 

the scalar field is on the local maximum of the potential and the mass squared m^. is negative. As we will see later, 
the condition (|II.39|) is nothing but the condition for stability of the de Sitter space. 

We have rewritten the action (|II.2[) of F(R) gravity in a scalar-tensor form (|II.29p . Inversely, it is always possible 
to rewrite the action of the scalar-tensor theory 

S = J d'x^—g | ±R - - V{tp) | , (IL41) 

as the action of F(R) gravity [l6j]. By comparing (|II.41|) with (jll.29[) . we find 

ip =^a, V{y) = v(^y\. (11.42) 



One now uses the conformal transformation 

<v-> e±K ^s^> ( IL43 ) 

to make the kinetic term of <p vanish. Then, we obtain 

S = / d 4 x^j { ° R - e ±2K *Vlv(<p) } . (11.44) 



2k 

The action pi.44[) is often called the "Jordan frame action" , whereas the action (|II.41|) is the "Einstein frame action" 
due to either the non-minimal coupling or the standard coupling in front of the scalar curvature. Because tp now 
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becomes an auxiliary field, which does not have a kinetic term, one may delete tp by using the corresponding equation 
of motion: 

R = e ±K *Vi (^K 2 V(<p) ± 2k^V'(<p)J , (11.45) 

which can be solved with respect to (p as a function of R as ip — <p(R). Thus, we can rewrite the action (|II.44[) in the 
form of F(R) gravity: 

S = [ d 4 x^=g~F(R) , F(R) = C±KV(R ^ R _ e ±2^(R)^y / /mx _ (11.46) 
J 

Note that one can rewrite the scalar-tensor theory f|II.41[) only when the sign in front of the kinetic term is negative, 
that is, when the scalar field is canonical. One cannot rewrite the action if the sign is positive, as 

S = J d i x^—g^R+ \d^<p-V(<p)\ , (11.47) 

which corresponds to the phantom scalar [l7j . 

Note that, for the metric transformed as in (|II.43[) . even if the Einstein frame universe is in a non-phantom phase, 
where the effective EoS w c e in (|II.4[) is larger than —1, the Jordan frame universe can be, in general, in a phantom 
phase. In other words, despite the mathematical equivalence between two frames there occurs some kind of physical 
non-equivalence [l8[ . More precisely, if physics in one frame is matched with observations and describes the observable 
accelerating universe, then another frame of physics becomes extremely unconventional. For instance, instead of the 
acceleration, the observer sees the deceleration, and matter couples with scalar, etc. 



3. Matter instability 

Let us discuss the matter instability issue in modified gravity. It is related to the fact that the spherical body 
solution in general relativity may not be the solution in modified gravity theory. The matter instability may appear 
when the energy density or the curvature is large compared with the average density or curvature in the universe, as 
is the case inside of a planet. Multiplying g^ v with Eq. (|II.5[) . one obtains 

F^(R) F'(R)R 2F(R) k 2 

DR + F~f 2 ~y(R~) R + 3FM(R) ~ 3F&(R) ~ 6F^{R) Tmattor ' (IL48) 

Here, T mat tcr is the trace of the matter energy-momentum tensor: T mat tcr = 7m a tt C r p p - We also denote d n F(R)/dR n 
by F<> n \R). Let us now consider the perturbation from the solution of general relativity. We denote the scalar 
curvature solution given by the matter density in the Einstein gravity by Rb ~ (K 2 /2)p mat t G r > and separate the 
scalar curvature R into the sum of Rb and the perturbed part R p as R = Rb + R P (\Rp\ -C \Rb\)- Then, Eq. (1II.48|) 
leads to the following perturbed equation: 

° - DRh+ F^(Rb) VpRbV Rb+ 3F(V(R b ) 

2t ' Uh 1 ,h ' OR p + 2^^-W p R b WPR p + U(R b )R P . (11.49) 



3F( 2 )(i? fc ) 3F( 2 )(i? b ) p F( 2 )(i? b ) 



Here, U(Rb) is given by 



"<*>-(a-»)w-* + ! 

F( 1 \Rb)F^(Rb)Rb F^(Rb) 2F(R b )F^(R b ) F^(R b )R b 
3F( 2 )(i? fc ) 2 3F(2)(i? b ) + 3F(2)(i? b )2 3F(2)(iT- b )2 



(11.50) 



It is convenient to consider the case that Rb and R p are uniform; that is, they do not depend on the spatial coordinates. 
Thus, the d'Alembertian can be replaced with the second derivative with respect to the time coordinate: \3R P — > 
—dfR p and Eq. (jll.50[) has the following structure: 

= -d 2 R p + U(R b )R p + const. . (11.51) 
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Then, if U(Rb) > 0, R p becomes exponentially large with time t: R p ~ eV^ 1 ^' and the system is unstable. In the 
1 /R model [2CH22| , because the order of the mass parameter m p is 

yT 1 ~ 10 18 sec ~ (lO-^eV)^ 1 , (11.52) 

one finds 

U(R b ) = -R b + ^| ~ 4 ~ (l0- 26 sec)^ 2 (Pe***) 3 , R b „ (lO^ec)" 2 ( Ei^Sl) . (11.53) 

Thus, the model is unstable and it presumably decays in 10 -26 sec (for a planet-sized object). On the other hand, in 
a 1/R + R 2 model @, one gets 

U(Rq)~^>0. (11.54) 

Then, the system could be unstable again, but the decay time is ~ 1,000 sec, that is, macroscopic. For the model 
(l5j . U(Rb) is negative 

UW~-£±*j^<0. (11.55) 

Therefore, there is no matter instability. Note that the models considered in the following sections successfully pass 
the matter instability test. 



4- Unified theories without a flat space solution 

We now review the F(R) models which do not have a flat space solution because of the non-analytical behavior at 
zero curvature. Some of these such models may induce large corrections to Newton's law. 

Modified gravity with negative and positive powers of the curvature 
In 0, the following action was proposed: 

F(R) = R - _° n +b(R- A 2 ) m . (11.56) 

Here n, to, c, and b are constants. In general, n and to may be fractional. Let us show that the above model leads to 
an acceptable cosmic speed-up and is consistent with some of the solar system tests. 
For the action pi. 561) . Eq. pi. 9ft has the following form: 

c{(n + 2)R - 2Ai} , , . , . 

= R - ; D I + b {(2 - m)R - 2Aa} (i? ~ A2) ■ (IL57) 

[R — AiJ 

Especially when n = 1, m = 2, and Ai = 0, one gets 



R = R ± = . (11.58) 

If bA 2 >, c < 0, and b 2 A 2 + 3c (1 + 2£A2) > 0, both of the solutions express the de Sitter space. Thus, the natural 
possibility for the unification of early-time inflation with late-time acceleration appears. Here, higher-derivative terms 
act in favor of early-time inflation, and the 1/R term supports a cosmic speed-up. 

We now neglect the contribution from matter. When n = 1, m = 2, and Ai = A2 = in (III.56|) and the curvature 
is small, we obtain the solution of pi.7|) : a a t 2 , which is consistent with the result in (20l - l22| . 

If the present universe expands with the above power-law, the curvature of the present universe should be small 
compared with that of the de Sitter universe solution in pi.58[) . As the Hubble rate in the present universe is 10 -33 eV, 

the magnitude of the parameter c, which corresponds to /1 4 in [2p| - [2^ j . is on the order of (lO _33 eV) . 

Let us consider the more general case that F(R) is given by (|II.56[) when the curvature is small and Ai = 0. 

(ti + 1)(2ti + 1) 

Neglecting the contribution from the matter again, solving (|II.7|) . we obtain a cut ™+ 2 (see pi.20p ). It is quite 
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remarkable that actually any negative power of the curvature supports the cosmic acceleration. This gives the freedom 
in modification of the model to achieve the better consistency with experimental tests of Newtonian gravity. 

(m-l)(2m-l) 

On the other hand, when the scalar curvature R is large, one obtains a a i m - 2 (see (III.20p ). When 

^"^m-™ -1 ^ ^ 0> ^ ne un i verse is shrinking but with a change in the arrow of time by t — >• t s — t, the inflation occurs 
with the inverse power-law. At t = t s , the size of the universe diverges. It is remarkable that when m is fractional 
(or irrational) and 1 < m < 2, the expression of a is still valid and the power becomes positive, the universe evolves 
with the (fractional) power-law expansion. 

It is interesting that the model in (|II.56j) is more consistent that refs. (20l - [22l | (for a discussion of the Newtonian 
limit in 1/R theory, see [23j]) as it may pass (at least, some) solar system tests. An example of the tests is matter 
instability as discusses in the previous subsection. Indeed, in 19], small gravitational object like the earth or the 
sun in the model [2fj| - l22j | is considered. It has been shown that the system becomes unstable. This may cause an 
unacceptable force between distant galaxies [IHIIil- As shown in the previous subsection (see ref. Q), however, by 
adding the positive power (higher than 1) of the scalar curvature to the action, the instability could be significantly 
improved. Furthermore, the account of quantum effects in modified gravity also acts against instability [26| . 

It has also been mentioned in [27| that the 1/R model [2014221 ] . which is equivalent to scalar-tensor gravity, is 
ruled out as a realistic theory due to the constraints to such Brans-Dicke-type theories. In it was shown that 
by adding the scalar curvature squared term to the action, the mass of the scalar field can be adjusted to be very 
large and the scalar field can decouple. Thus, the modified gravity theory (|II.56|) (after some fine-tuning Q) passes 
most of the solar system tests. In addition, at precisely the same parameter values, the above modified gravity 
has Newtonian limit that does not deviate significantly from that in general relativity. Nevertheless, as has been 
mentioned at the very beginning, such theory has non-analytical behavior at zero (or some constant) curvature which 
may be phenomenologically unacceptable. 

lni? gravity 

Other gravitational alternatives for dark energy may be suggested along this same line. As an extension of the 
theory (|II.56|) , one may consider the model [28j containing the logarithm of the scalar curvature R: 

R 

F(R) = R + a'ln— +/3iT\ (11.59) 
M 

We should note that the m — 2 choice simplifies the model. Assuming the scalar curvature is constant and the Ricci 
tensor is also covariantly constant, the equations (|II.9[) are as follows: 

= 2F(R) - RF'(R) = F(R) = R + 2a' In - a' . (11.60) 

If a' > 0, F(R) is a monotonically increasing function and lim^^o F(R) — — oo and lim^_ ! . +00 F(R) = +oo. There 
is one and only one solution of (jll.601) . This solution may correspond to the inflation. On the other hand, if a' < 0, 
limji^o F (R) = oo and lirnR_ s .+ 00 = +oo. Because F'(R) = 1 + %■, the minimum of F(R), where F'(R) = 0, is 
given by R = —2a'. If /(—2a') > 0, there is no solution of (|II.60[) . If /(—2a') = 0, there is only one solution and if 
f(—2a') < 0, there are two solutions. Because /(—2a') = —2a' ^1 — In , there are two solutions if —^t > e. 

Because the square root of the curvature R corresponds to the rate of the expansion of the universe, the larger of the 
two solutions might correspond to the inflation in the early universe and the smaller one to the present accelerating 
universe. 

We can consider the late FRW cosmology when the scalar curvature R is small. Solving (|II.7j) . it follows that the 
power-law inflation could occur: a oc 1 5 . 

One may discuss further generalizations like [28[ (see also [29j]) 

F(R) =R + 7 iT" (in - . (11.61) 

Here, n is restricted by n > —1 (m is an arbitrary integer) in order that the second term could be more dominant 
than the Einstein term when R is small. For this model, we find 

(r»+l)(2r» + l) 

a - t — — . (11.62) 
This does not depend on m. The logarithmic factor is almost irrelevant. The effective w c g (|II.4p is given by 

6n 2 + 7n-l 

w cft = -^7 — '. TTTo — T~n ■ 11.63) 

3(n + l)(2n + 1) 
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Then, w e g can be negative if 



-| n I / r 7 f i 

Kn< — oru> — - — = 0.1287- ••. (11.64) 

2 12 v ' 



From (|II.62|) . the condition that the universe could accelerate is given by i n + 1 )i^+ 1 ) -> ^ j s: 

n > ~ l+ 2 ^ = 0.366 • • • . (11.65) 

Clearly, the effective EoS parameter w may be within the existing bounds. Thus, the logarithmic model may also 
propose the unified description of inflation with dark energy era. 



5. Viable modified gravities 

Using the previous arguments, let us review the viable models of gravity that could unify the accelerating expansion 
of the present universe and the inflation of the early universe. 

The realistic F(R) models unifying inflation with dark energy 

In refs. [30l - l32l ]. viable models of F(R) gravity unifying the late-time acceleration and the inflation were proposed. 
To construct these models, we have required several conditions: 



1. To generate the inflation, one requires 



lim f{R) = -K . (11.66) 

R— >oo 



Here, Aj is an effective cosmological constant at the early universe; therefore, it is natural to assume A, 3> 
(lCT 33 eV) 2 . For instance, it is natural to have A l ~ io 20 ~ 38 (eV) 2 . 

2. To enable the current cosmic acceleration to be generated, the current f(R) gravity is considered to be a small 
constant, that is, 

f(Ro) = -2R , f'(Ro)~0. (11.67) 

Here, Rq is the current curvature Rq ~ (lO~ 33 eV) 2 . Note that Rq > Rq due to the contribution from matter. 
In fact, if one can regard f(Ro) as an effective cosmological constant, the effective Einstein equation gives 
Ro = Ro — K 2 T'mattor- Here, T mat t C r is the trace of the matter energy-momentum tensor. Note that f'(Ro) need 
not vanish completely. Because we consider the time scale from one billion years to ten billion years, we only 
require |/'(-R )| < (lCT 33 eV) 4 . 

Instead of the model corresponding to (|II.66[) , one may consider a model that satisfies 

lim f(R) = aR m , (11.68) 

R— >oo 

with a positive integer m > 1 and a constant a. To avoid the anti-gravity f'(R) > — 1, a > and, therefore, 
f(R) should be positive at the early universe. On the other hand, Eq. pi.67[) shows that f(R) is negative at 
the present universe. Therefore, f(R) should cross zero in the past. 



3. The last condition is 



Urn/(#) = 0, (11.69) 



which means that there is a flat space-time solution. 
One typical model proposed in ref. [3(| that satisfies the above conditions is 
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Here, a, (3, and 7 are positive constants and n is a positive integer. Eq. (III.70|) gives [3C 

l/r. 

I / I \ / >-. \ I 

i?0 = 



1 + \ 1 



(11.71) 



and, therefore, 



f(R ) ~ -2i? = -» 1 



1- & ) 22 



/3 7 



Then, it follows that 



O E? T}~' Zn R A T?'^ T}~' Zn T? n ~ 1 ~, Of? E?~ 2n r>n — 1 



(11.72) 



(11.73) 



In the model (|II.70p , the correction to Newton's law is small because the mass m a (|II.30[) is large and is given by 
to 2 ~ ^Q-i6o+i09n e y2 j n ^ e so l ar system and m 2 ~ io~ 144 + 98 " cV 2 in the air on the earth [30j |. In both cases, the 
mass m a is very large if n > 2. 

Another model corresponding to (|II.66I) was proposed in 32 1 



f(R) = — cxq ( tanh 
One now assumes 

and 



b (R-Ro) 



tanh 



b Rc 



— ai ( tanh 



bi (R-Ri) 



tanh 



bjRj 



Ri >• R , ai > a , 6/ <C 6 . 



6/i?/ > 1 . 



When R — !> or R < i? , Ri, f(R) behaves as 



a b 



OLibi 



2 cosh 2 (^f*) 2 cosh 2 (^) 
and /(0) = again. When R Rj, it follows that 



R. 



f(R) -2A/ ee -a 1 + tanh 



boRo 



— aj I 1 + tanh 



—a/ 1 + tanh 



6/i?/ 



On the other hand, when Rq <C R <Si Ri, one gets 



f(R) -> -a 



1 + tanh 



b Ro 



aibjR 



2 cosh 2 



Here, we have assumed the condition (|II.76j) . One also finds 

a b 



f(R) = - 



-2A ee -a 



1 + tanh 



&oi?o 



2 cosh 



2 (M^M) 2cosh 2 (M^«l)) 



which has two valleys when R ~ i?o or R ~ Rj. When R = Ro, 

aibi 



f(R ) = -a b 

On the other hand, when R = Rj, it follows that 

/'(#/) - -a/6, 



2,osh 2 ( MJ Y J?/) ) 



a b 



2 cosh 



2 f boiRo-Rt] 



> -a/6/ - a 6 . 



> -a/6/ - a 6 . 



(11.74) 

(11.75) 
(11.76) 

(11.77) 

(11.78) 
(11.79) 
(11.80) 

(11.81) 

(11.82) 
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Due to the condition (|II.35|) required to avoid the anti-gravity regime, one obtains 

a I b I + a Q b <2. (11.83) 
In the solar system domain, on or inside the earth, where R 3> Rq, f(R) can be approximated by 

f(R) ~ -2A off + 2ae- b ( R - Ro '> . (11.84) 
On the other hand, because R <C R <C Ri, by assuming Eq. (|II.76|) . f{R) (|II.74[) could also be approximated by 

f(R) ~ -2A + 2ac- b ^ R - R ^ , (11.85) 



which has the same expression, after having identified Ao = A c ff and &o = b. Thus, one may check the case (jlL 
only. The effective mass has the following form 



p b{R-R ) 
2 c 

m„ - 



Aab 2 



(11.86) 



which could be very large, that is, m 2 ~ 10 1,000 eV 2 in the solar system and m 2 ~ iqio, 000,000,000 e y2 - n a j r 
surrounding the earth. Thus, the correction to Newton's law becomes negligible. 
One may consider another model (30j : 

f(R) = fc^i ±^2 (IL87) 

fo + h{(R-Rof k+1 +Rl k+1 } 

It has been shown [30] that for k > 10, such modified gravity passes the local tests. It also unifies the early-time 
inflation with the dark energy epoch. In (III.87[) . Rq is the current curvature Rq ~ (lO~ 33 eV) . We also require 

h~\, h = Ti (11.88) 
Here Aj is the effective cosmological constant in the inflation epoch. When R ^> A i; f(R) (|II.8T[) behaves as 

f( R )--J l +Jf^I- (11-89) 

The trace equation, which is the trace part of (|II.5[) . is as follows: 

dDf(R) =R+ 2f(R) - Rf'(R) - K 2 T matter . (11.90) 

Here, T mat tcr is the trace part of the matter energy- momentum tensor. For the FRW metric (111.61) . one finds 

R~ (t s _t)- 2 /( 2 ™+ 3 ) ) (H.91) 

which diverges at a finite future time t = t s . By a similar analysis, we can show that if f(R) behaves as f(R) ~ R e for 
large R with a constant e, a future singularity appears if e > 2 or e < 0, which is consistent with the analysis (|II.34|) 
via the scalar-tensor form of the action. Conversely if 2 > e > 0, the singularity does not appear. Thus, adding the 
term R 2 f(R), where lim^^o f(R) = c\, hnifl->oo f{R) = C2, to f(R) (jll.87[) . the future singularity (III.91|) disappears. 
More details about future singularities in modified gravity will be given in the fourth chapter. 
Let us investigate the above situation in more detail. Assume 

f(R) ~F + W , (11.92) 

when R is large. Here Fg and F\ are constants where Fq may vanish but F\ 7^ 0. In the case of pi. 891) 



F = ~, *i = 4 , e = - (2n + 1) . (11.93) 
h h 



Under the assumption pi.92[) . the trace equation pi.90|) gives 



aFiDfl- 1 = {* , „ R£ w i en 6 < ? or e = ? . (IL94) 

1 (2 — e) Fiic when e > 1 or e f= 2 v ' 
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In the FRW background pi.6|) . when the Hubble rate has a singularity as 



with constants ho and f3, the scalar curvature R = QH + 12H 2 behaves as 

when/3>l 

R ~ } when 13 = 1 . (11.96) 

[ J^kr when (3 < 1 

In (|IL95|) or pL96)l . the /3 > 1 case corresponds to Type I (Big Rip) singularity in P, HI], 1 > (3 > to Type III, 
> (3 > -1 to Type II, and f3 < -1 (but (3 ^ integer) to Type IV. 

The classification of the finite-time future singularities used above is given in ref . [34[ : 

• Type I ("Big Rip") : For t — > t s , a — > oo, p c g — > oo and |p c ff| — > oo. This also includes the case of p e ff, Pes 
being finite at t s . This singularity has been discovered in ref. Its manifestations in different models have 
been studied in refs. [34j |. 

• Type II ("sudden") (35|: For t —> t s , a — > a s , p e g — > p s and |p e gf| — > oo. 

• Type III : For t — >■ t a , a —> a s , p e g °o and \p c g\ — > oo. 

• Type IV : For f — > f s , a — > a s , /9 e g — > 0, |p c ff| — > and higher derivatives of H diverge. This also includes 
the case in which p e g (fi e &) or both of p e s and p c g tend to some finite values, whereas higher derivatives of H 
diverge. 

Here, p e g and p e g are defined by 

p off EE \H 2 , p cff EE - -L (2ff + 3H 2 ) . (11.97) 

Substituting (jll.96[) into pi.94[) . one finds that there are two classes of consistent solutions. The first solution is 
specified by (3 = 1 and e > 1 (but e ^ 2) case, which corresponds to the Big Rip (ho > and t < t s ) or Big Bang 
(ho < and t > t s ) singularity at t = t s . Another solution is e < 1, and j3 = — e/ (e — 2) (—1 < (3 < 1) case, 
which corresponds to (|II.91|) and to the Type II future singularity. In fact, we find that e = — 2n — 1 and, therefore, 
(3+1 = -2/(2n + 3). We should note that when e = 2, that is, f(R) - R 2 , there is no singular solution. Therefore, 
if the above term R 2 f(R), where lim/?^o f(R) — c i, li m fl-4-oo f(R) = c 2, is added to f(R) in (|II.87|1 . the added term 
dominates, and the modified f(R) behaves as f(R) ~ R 2 - The future singularity (|II.91[) disappears. We also note 
that if the i? n -term with n = 3,4, 5, • • ■ is added, the singularity becomes (in some sense) worse because this case 
corresponds to e = n > 1 , that is the Big Rip case. Using the potential that appears when we transform F(R) gravity 
to scalar-tensor theory [3y|, it has been found that the future singularity may not appear in the case where < e < 2. 
Thus, to avoid the finite-time future singularity, one has to add the i? 2 -term to the above viable unification models. 

/(*) - ^"r^ +cR 2 > (n.98) 



1 + 7i? r ' 



f(R) = -oo ( tanh ( Mg^gg) ) + tanh 



2 7 V 2 



-a/ 



( tanh ( b i( R - R i) ) + tanh ( b J^L) \ +c R 2 i (H.99) 



2 \ 2 



m = ( r" } 1 ° ^ + ^ 2 ■ (n-ioo) 



/o + /i { (H - i2o) 2fc+1 + J2g* +1 } 



The addition of this term was proposed first in ref. [37|, where it was shown that, in this case, the Big Rip singularity 
disappears. Moreover, such a term, which effectively eliminates future singularity also supports the early-time inflation. 
In other words, adding such a R 2 term to the gravitational dark energy model (for instance, to viable dark energy 
models [lH [Hj]) may lead to the emergence of an inflationary phase in the model as was first observed in ref. Q. 
In the case when a model already contains the inflationary era, its dynamics will be changed by the i? 2 -term. The 
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investigation showing that the i? 2 -term cures all types of future singularities was carried out in refs. [36l. I39M41] . In 
fact, it has been recently realized that some phenomenological problems [421 - I46T ] of F(R) dark energy (like achieving 
a consistent description of neutron stars) may be resolved in the presence of the R 2 term. Moreover, the traditional 
phantom/quintessence (fluid/scalar) dark energy models often bring the universe to a finite-time future singularity. 
It was demonstrated in ref. [47] that the natural method for eliminating the singularity in these models is again the 
addition of the R 2 term. In other words, resolving future singularities in specific dark energy models requires (at least 
partly) modified gravity. We discuss these questions in more detail in the fourth chapter. Note also that in the above 
models, there is no matter instability (Eq. pi.51j l). 

It should be noted that the exponential-type model may represent a very interesting proposal for the resolution of 
the dark energy problem. Some versions of such a theory are 

F(R) =R + a {e- bR - l) or F(R) = R-a - ) , (11.101) 

with constants a, b, and Rq. This model, which was introduced in ref. [HJ provides the accelerating cosmological 
solutions without a future singularity. As a rule, the accelerating late-time cosmology leads asymptotically to a de 
Sitter universe. It was shown [48| that such a model describes a realistic dark energy epoch that is compatible with 
local and observational tests. Moreover, the number of free parameters is not big. Of course, one can consider different 
modifications of the above model. 

A more complicated but viable F(R) model 

The above analysis shows that, to obtain a realistic and viable model, F(R) gravity should satisfy the following 
conditions: 

1. When R 0, the Einstein gravity is recovered, that is, 

F(R) -» R that is, - . (11.102) 

R z R 

This also means that there is a flat space solution as in (|II.69[) . 

2. As will be discussed after Eq. (IIII.154[) . there appears a stable de Sitter solution, which corresponds to the 
late-time acceleration and, therefore, the curvature is small R ~ Rl ~ (l0 _33 eV) . This requires, when 



R~Rl, 



^ = f 0L - f 1L (R - R L ) 2n + 2 +o((R- R L f n+2 ) . (11.103) 



Here, /ql and hh are positive constants and n is a positive integer. Of course, in some cases this condition may 
not be strictly necessary. One can use it, for instance, to avoid a future singularity. 

3. As will also be discussed after Eq. 1 54[) . there appears a quasi-stable de Sitter solution that corresponds 
to the inflation of the early universe and, therefore, the curvature is large R ~ Ri ~ (io 16 ~ 19 GeV) 2 . The de 
Sitter space should not be exactly stable so that the curvature decreases very slowly. It requires 

1 ' U) - hi - hi (R - Rif m+1 +o((R- Rj) 2 ^ 1 ) . (11.104) 



R 2 

Here, hi an d hi are positive constants and m is a positive integer. 
4. Following the discussion after (jll.33[) . when R oo, to avoid the curvature singularity, it is proposed that 

F(R)^f oc R 2 that is S^->/oo. (11.105) 

R z 

Here, is a positive and sufficiently small constant. Instead of (|II.105|i . we may take 

F(R) -> f^R 2 -' that is (11.106) 

rt z R e 

Here, is a positive constant and < e < 1. The above condition (|II.105|) or (III. 106[) prevents both the future 
singularity and the singularity due to large density of matter. 
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*- R 



FIG. 1: The qualitative behavior of F ^P versus R in a viable model 



5. As in (|II.35j) . to avoid the anti-gravity, we require 



F'(R) >0. 



which is rewritten as 



2 
R 



6. Combining conditions pi. 102[) and pi.l07[) . one finds 

F(R) > 0. 

7. To avoid the matter instability in (III.51[) [X9j . we require 

TKR \ = ?±- F (1) (Rb)F^(R b )R b _ F^(R b ) 2F(R b )F^(R b ) _ F^(R b )R b 
[ b) 3 3F( 2 )(R b ) 2 3FW(R b ) 3F( 2 )(i? 6 ) 2 3F( 2 ) (R b ) 2 



(11.107) 
(11.108) 

(11.109) 

(11.110) 



The conditions Q] and [2] show that an extra, unstable de Sitter solution must appear at R = R e (0 < R e < Rl) (see 
FigHJ. The universe evolution will stop at R — Rl because the de Sitter solution R = Rl is stable; the curvature 

F(R) 

never becomes smaller than Rl and, therefore, the extra de Sitter solution is not realized. The behavior of ^ 2 , 
which satisfies the conditions (jII.102j) . (|II.103I) . (IIL104|) . pi.!05|l . and (|II.109|) is given in FigHJ 
An example of such F(R) gravity is given in [49j| 



F(R) 
R 2 



X m (Ri;R) — 



|(A m (Ri; R) — X m [Rj] Ri)) (X rn (Ri; R) — X m [Rj] Rl)) 



2n+2 



+X m (Ri; i?i) X m (Rj; Rl) 

(2m + 1) Rj m 
(R - R!) 2m+1 + Rj m+1 ' 



2n+2 



(11.111) 
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Here, n and m are integers greater or equal to unity, and n, m > 1 and i?i is a parameter related with i? e by 

Ym p \ (gn + 2) X (fl f ; gx) X gi) + X (ifr ; 

X i? e ) = — . (11.112) 

It is assumed that 

< Ri < R L < i?/ . (11.113) 
X (i?j; i?) is a monotonically decreasing function of R and, in the limit R — > 0, X (Ri; R) behaves as 

X(Rr,R)^^, (11.114) 

which shows that in the limit J2 — 5> 0, F(R) (|II.111[) reproduces Eq. (|II.102j) and, therefore, the condition Q] is satisfied. 
On the other hand, when R — > oo, 

X(Rr,R)^ {2m +l\T r ' ->(), (11.115) 

and, therefore, -F(i?) behaves as 105[) and the condition 2] is satisfied. 
When i? ~ Rl, we find the behavior of ()II.103|) . where 

i 

foL — (x m (Ri; Ri) X m (Ri; R L ) 2n+2 + /^™ +3 | , 

/il = ^ — g |X m (J?/; i?i) X m (Ri;Rl) n + f^ +3 j (X m (Ri; Ri) — X m (Rj; Rl)) 

(2m+l) 4(m+1) {Rj (R L -Ri)} im{m+1} 

x ~~ ; m ; i^+n • ( IL116 ) 

{(ij i -ij J ) am+1 + fl? m+1 } 

Then, the condition [2] is satisfied. 

On the other hand, when R ~ i?/, we also find the behavior of (|II.104[) . where 

foi = ^(X m (Ri: Ri) — X m (Rr, Ri)) (X m (Ri; Ri) — X m (Ri; Rl)) 2u+2 

i 

+X m (Ri;Ri) X m (Ri; Rl) 2u+2 + /^I l+3 | 

2?7i -(- 1 r 2 _|_2 

■f 11 = n 9 m+2 i ("^ m (~^ 7 ' ~ (~^ J ' ~^)) (X m (Ri;Ri) — X m (Ri; Rl)) 

Rj - y 

+X m (Ri;Ri) X 711 (Ri;Rl) n + f^ +3 > Ux m (Ri;Ri) — X„ 7 (Ri;Rl)) n 

+ (2n + 2) (X ro (i?/; Jfr) - X m (i?/; i?i)) (X m (Ri;Ri) - X m (R r , R L )) 2n+1 } . (11.117) 

Thus, the condition [3] is satisfied. 

Let us now investigate the mass of the scalar field a to determine whether or not the Chameleon mechanism [50l [5lj 
works or not. For this purpose, one investigates the region i?i < Rl <U< Ri- In this region, X m (Ri;R) can be 
approximated as 

X m (Ri; R) ~ — , X m (i?/; i? x ) ~ — , X m (Ri;R L ) ~ — . (11.118) 
ri /ti it/, 

Thcnj F(R) can be approximated as 

M f > , _ i?i + (2n + 2)fl L / 0O 2 "- 2 Z^ 2 - 2 
i? 2 ~ /o ° + i?' /n " (2n + 3)i? 1 i?i"+ 2 ~ i? 2 ^ 2 • [ ' 

In the last equation, it was assumed that R\ ~ If fooR "C / n or 



1 



» JtfZf^ , (11.120) 
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one gets 

m -~vz- (IL121) 

If we consider the region inside the earth, because 1 g ~ 6 x 10 32 eV and 1 cm ~ (2 x 1CP 5 eV) , the density is about 

Pmattcr ~ lg/cm 3 ~5x 10 18 eV 4 . This shows that the magnitude of the curvature is R ~ /t 2 pmatter ~ (lCP 19 eV) 2 . 

In the air on the earth, one finds p ma ttcr ~ 10~ 6 g/cm 3 ~ 10 12 eV 4 , which gives i?o ~ K 2 p m attcr ~ (lO~ 25 eV)~\ In the 
solar system, there could be interstellar gas. Typically, in the interstellar gas, there is one proton (or hydrogen atom) 
per 1cm 3 , which shows p m attcr ~ 10~ 5 eV 4 , Rq ~ 10~ 61 eV 2 . Then, the condition 120[) can be easily satisfied; 
for example, we can choose -j— ~ MeV 2 . Then, the Compton length of the scalar field becomes very small, and the 

Joe 

correction to Newton's law is negligible. 

One can now check the matter (in)stability issue [llj]. In the earth or sun, the condition i?i < Rl <C R <C Ri is 
satisfied and, therefore, F(R) behaves as (|II.119|) . Then, U(Rb) in (jll.1101) is given by 

U(Rb)~ -J?- <0. (11.122) 

Because U(Rb) is negative, the matter instability does not occur in the model under consideration. 

The condition [5] remain to be checked. Because F(R) (|II.111|) satisfies Eqs. (|II.102|) and (|II.105p . it is clear that 
Eq. (|II.107p or ([II.108P and, therefore, the condition [5] are satisfied when R — > or R — > oo. Because F & is a 
monotonically increasing function of R in the region R e < R < Rr,, Eq. (|II.107[> or (III. 108|) is trivially satisfied in 
the region. In the region Ri < R L <C R < Ri, because behaves as qiI.119|> . Eq. (|II.107|> or (TH.108P is again 
satisfied. Then, the condition [5] seems to be satisfied throughout the whole region. 

After the inflation at R = Rj, the radiation and the matter are generated. If the matter and radiation energy 
densities dominate compared with the contribution from f(R) = F(R) — R, that is, only the first term in (|II.24p 
dominates, a radiation/matter dominated-universe will be realized. The radiation and matter densities decrease 
rapidly compared with the contribution from the f(R) term at late-time universe. When curvature arrives at R ~ Rl, 
the late-time acceleration occurs. 

The condition pi.1061) indicates that the R 2 term appears in the very high curvature region. The R 2 term will 
generate an inflation besides the inflation at R = Rj if the universe started with a very high curvature R 3> Rj. 
Because the inflation due to the i? 2 -term is unstable, the inflation at the very early stage would stop; however when 
the curvature decreases and reaches R ~ Ri, the inflationary phase will occur again. 

Thus, we demonstrated that a number of viable F(R) gravity models may explain the early-time inflation in 
addition to the dark energy epoch in a unified way. It will be shown in the third chapter that when such a unification 
is problematic, one can use cosmological (partial/complete) reconstruction to ensure it. 



B - f(S) gravity 
1. General properties 

Let us consider another class of modified gravity where the arbitrary function added to the action of the general 
relativity depends on topological Gauss-Bonnet invariant. It has been shown recently that such a class of modified 
gravity is closely related with (super)string theory. The starting action is given by |52h54| : 

S = J d^xj—g [j^R + f{G) + ^matter) . (11.123) 

Here, £ m attcr is the Lagrangian density of matter and Q is the Gauss-Bonnet invariant: 

G = R 2 - AR^R» U + R^frRP*" . (11.124) 
By the variation of the metric g^ Vl we obtain 

o - ^ {-R^ + \g» v R>) + TZu* + \g» v f{Q) - V'{G)RR>* V 

+4f'(g)R%R» p - 2f'(g)R^ T R" paT - Af'(Q)R t2pau R pa + 2 (V^f'iG)) R 
-2g^ {V 2 f\g)) R-A (V„V7'(a)) R vp - 4 (V P V7'(£)) R^ 

+4 (V 2 f'(G)) R» v + Ag^ (V p V a f'(G)) R pa - 4 (V p V CT /'(£)) R wv ° ■ (H.125) 
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It is interesting to note that some memory of the fact that the Gauss-Bonnet term is topological invariant is contained 
in the above equation of motion. Indeed, it will be seen that such an equation of motion does not contain the 
higher-derivative terms as, for instance, the theory in the previous section. 

Choosing the spatially flat FRW universe metric (|II.6[) . the equation corresponding to the first FRW equation has 
the following form: 

= -^H 2 - f(G) + gf(G) - 24Gf"(g)H 3 + Pmattcr . (11.126) 

K 

In the FRW universe (|II.6j) . Q looks like: 

Q = 24 (h 2 H + i? 4 ) . (11.127) 
Then, from Eq. (|II.126|1 . the FRW like equations, as in the Einstein gravity case (|II.3|) . we find 

Pcff = \h 2 , p% = -1 h,H 2 + 2H) . (11.128) 



Here, 



p° oS = -f(G) + gf(g)-mf"(G)H 3 + Pnillttc 



Pes = f(Q) - Qf'(G) + i^f(G) + 8H 2 Gf"(G) + 8H 2 g 2 f"(G) + p matter . (11.129) 

When Pmattcr = 0, Eq. (jll.1261) has a de Sitter universe solution where H, and therefore G, are constant. For H = H , 
with a constant H , Eq. (jll.1261) turns into 

= - \h 2 + 24tf 4 /' (24i7 4 ) - / (24if 4 ) . (11.130) 
As an example, we consider the model 

f(G) - /o \Gf , (11.131) 
with constants /o and (3. Then, the solution of Eq. (|II.130|) is given by 



H* = . (11.132) 

24 (8 (n-l)« 2 /o)^ 

For a large number of choices of the function f(G), Eq. (|II.130[) has a non-trivial (H ^ 0) real solution for H 
(de Sitter universe). The late-time cosmology for above theory without matter has been discussed for a number of 
examples in refs. [52h54| . We will refer to these studies in the following. 

The de Sitter universe can be expressed as an accelerating, expanding universe, 

ds 2 - -di 2 + e 2Hot i dxl f ■ ( IL133 ) 

i=l,2,3 

Changing the coordinates, the de Sitter universe can be described by the static metric, which can be obtained by 
putting M = in the metric (|II.10[) (and choosing the minus sign in =p): 



r 



'Is 1 ( I - —J dt 2 + f 1 - ^) dr 2 + r 2 dfl 2 . (11.134) 



The length parameter L can be identified with 1/Hq. Note that there is a cosmological horizon at r = L. When 
r ~ L, by the Wick-rotating time coordinate t = it, the metric has the following form: 

ds 2 ~ 2 fl - dr 2 + + l 2 dfl 2 . (11.135) 

Changing the radial coordinate r to p as 

= L (l - ^\ , (11.136) 
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the metric (|II.135|> can be rewritten as 



n 2 

d s 2 ~ ^-dr 2 + dp 2 + l 2 dn 2 . (11.137) 



To avoid the conical singularity at p = 0, the periodicity condition for the coordinate r should be imposed as 
r ~ t + 2irL. Thus, the Hawking temperature Tjj is defined as 

T H = ^. (11.138) 

Let us now calculate the entropy like in Eq. pi,18|) . The partition function Z 16j> can be evaluated from the action 
123|) by neglecting the matter Lagrangian £ ma tter an d substituting the metric of the de Sitter universe (III. 134|) : 

-lnZ~S = 4«j Q d tj o drr 2 ^ + f^ 

8tt 2 L 4 / 6 / 24 

^Z 2 +J [j2 



1 /24tt 2 1| 



6tt 2 T 4 V 



/ (3847r 4 T H ) . (11.139) 



Note that R = $ and Q = f£. Using (|rEL6]) and piI.148|l . the free energy T and the entropy S are 

1 /247r 2 T 2 



T = 6^7f ( v ^ Ji + ^ 384 " 4r H)) > (TLI4U) 



4_ _ /(384^ 4 T 4 ) 

H 



^ = --o^ - \ M 2567T 2 /- (384^ 4 T H 



8^ 2 L 4 /(— j+ 2567r 2 /'f-^l . ! 1 1.141) 



Because H = \, Eq. (|IL130j) can be rewritten as 



k 2 L 2 L 4J V^V V^ 4 
By using (III. 142|) . the entropy (jll. 141|) has the following expression 



■5 = ^+64^/'^ • (H.143) 



4G \L 4 J 

Here, Ah is the area of the cosmological horizon: An — 4nL 2 and G — When f(G) is a constant, which corresponds 
to the cosmological constant, Eq. (|II.143[> reproduces the standard result. 

One now considers the case p ma tter 7^ 0. Assuming that the EoS parameter w = Pmatter/ Pmatter is a constant, using 
the conservation of energy: 

Pmatter + 3-ff (^matter + Pmatter) = , (11.144) 

we find p — pqcl~ 3 ( 1+w \ When the function f{Q) is chosen as in (|II.131j) . again, if < 1/2, the f(G) term becomes 
dominant compared with the Einstein term when the curvature is small. If we neglect the contribution from the 
Einstein term in (|II.126|) . the following solution may be found 

a = aot h ° when ho > , a = ao (t s — t) h ° when ho < , 



4/3 

a 



{24\h 3 (-l + ho)\} P (h - 1 + 4(3) i{1+w > . (11.145) 



3(1 + w)' L (h -l)p 

Thus, the effective EoS parameter w e g (|II.4[) is less than — 1 if j3 < 0, and for w > —1, it is 



w °*=- 1+ Wo = - 1+1 -w> (IL146) 
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which is again less than — 1 for (3 < 0. Thus, if (3 < 0, we obtain an effective phantom with negative ho even when 
w > — 1. In the phantom phase, the Big Rip might seem to occur at t — t s . Near this Big Rip, however, the curvature 
becomes dominant and then the Einstein term dominates, so the f(G) term can be neglected. Therefore, the universe 
behaves as a — ao£ 2 / 3< - IU+1 - ) and as a consequence, the Big Rip does not occur. The phantom era is thus transient. 

The case when < (3 < 1/2 may also be considered. As (3 is positive, the universe does not reach the phantom 
phase. When the curvature is large, the f(G) term in the action (|II.123[) can be neglected and one can work with the 
Einstein gravity. Then, if w is positive, the matter energy density p ma ttor should behave as p ma ttcr ~ t , but f(G) 
goes as f(G) ~ t~ 4 ^ . For late times (large t), the f(G) term may become dominant as compared with the Lagrangian 
density of matter. Neglecting the contribution from matter, Eq. pi.l26[) has a de Sitter universe solution where H, 
and, therefore, G, are constant. If H = H with a constant H , Eq. (III. 126|) looks as (|II.130p . As a consequence, even if 
we start from the deceleration phase with w > —1/3, an asymptotically de Sitter universe emerges. Correspondingly, 
a transition also occurs here from deceleration to acceleration of the universe. 

Now, we consider the case in which the contributions from the Einstein and matter terms can be neglected. 
Eq. (|II.126|) reduces to 



If f(G) behaves as 131 [) . assuming 



= Gf'{G) - f(G) - 2AGf"(G)H s . (11.147) 

aot h ° when ho > (quintessence) , . 

a (t s — t) h ° when ho < (phantom) 



one obtains 



= (0 - 1) fcg (ho - 1) (ho - 1 + 4/3) . (11.149) 
As ho = 1 implies G = 0, one may choose 

h = l-4l3, (11.150) 

and Eq. (|II.4[) gives 

Therefore, if j3 > 0, the universe is accelerating (w e g < —1/3), and if (3 > 1/4, the universe is in a phantom phase 
(u> c ff < —1). Thus, we are led to consider the following model: 

HG)^fi\g\ Pi + fi\gf l , (n.152) 

where it is assumed that 

h>\, \>h>\- ( IL153 ) 

Then, when the curvature is large, as in the primordial universe, the first term dominates, compared with the second 
term and the Einstein term, and it gives 

-l>- cff = -l + ^^>4 (11.154) 

On the other hand, when the curvature is small, as is the case in the present universe, the second term in (|II.152|) 
dominates compared with the first term and the Einstein term and yields 

- off = -l + — (11.155) 

Therefore, theory (|II.152j) can produce a model that is able to describe inflation and the late-time acceleration of the 
universe in a unified manner. 

Instead of (|II.153I) . one may also choose /3i as 

J>A>0, (11.156) 
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which gives 

-^>w eB >-l. (11.157) 

Then, we obtain an effective quintessence epoch. Moreover, by properly adjusting the couplings fi and /; in 152[) . 
one can obtain a period where the Einstein term dominates and the universe is in a deceleration phase. After that, 
a transition occurs from deceleration to acceleration when the Gauss-Bonnet term becomes the dominant one. More 
choices of f(G) may be studied for the purpose of the construction of the current accelerating universe. Nevertheless, 
many non-linear choices for this function may be approximated by the above model. For instance, one can mention 
some realistic examples of f(G) gravity: 

f (G\ aiQn + hl f(c\ aiQn+N ± bl miw 

Let us address the issue of the correction to Newton's law. Let gm) be a solution of (jll. 125|) and represent the 
perturbation of the metric as g^ v = g(o)tiu +h flI/ . First, we consider the perturbation around the de Sitter background. 
The de Sitter space metric is taken as 5(o)/ji/j which gives the following Riemann tensor: 

R(o)^ P a = H (5(0)^,5(0)^ - 5(o)a«t5(oVp) ■ (11.159) 

The flat background corresponds to the limit of H — > 0. For simplicity, the following gauge condition is chosen: 
9$ V = V fo) V = °- Then E q- AIT251) gives 

= -j— 5 (V 2 V - 2ff 2 V) + T matterM „ . (11.160) 

The Gauss-Bonnet term contribution does not appear except in the length parameter 1/Hq of the de Sitter space, 
which is determined by taking into account the Gauss-Bonnet term. This may occur due to the special structure of the 
Gauss-Bonnet invariant. Eq. 160[) shows that there is no correction to Newton's law in de Sitter space and even in 
the flat background corresponding to Ho — > 0, regardless of the form of / (at least, with the above gauge condition). 
In contrast, the study of the Newtonian limit in l/R gravity indicates that significant corrections to Newton's law 
may appear as is shown in refs. [231 ] . For most l/R models, the corrections to Newton's law do not comply with solar 
system tests. 

By introducing two auxiliary fields, A and B, one can rewrite the action (jll.1231) as 



S = I <f aV^l + B{G-A) 



+f(A) + £ mattC r) • (11.161) 

By the variation of the action (|II.161|) with respect to B, it follows that A — Q. Using this in 161|) . the action 
(jll. 123|) is recovered. Alternatively, varying the action 161|) with respect to A in (|II.161[) . one gets B = f'(A), and, 
thus, 

S = J d'x^—g^R + f\A)Q - Aj'(A) + f(A)) . (11.162) 

By varying the action fjll. 162[) with respect to A, the relation A = Q is obtained again. The scalar is not dynamic as 
it has no kinetic term. One may add, however, a kinetic term to the action by hand 

S = J cfxV^g-^R - € -d^A + f'(A)g - Af'(A) + f(A)) . (11.163) 

Then, one obtains a dynamical scalar theory coupled with the Gauss-Bonnet invariant and with a potential. It is 
known that, in general, a theory of this kind has no ghosts and is stable. Actually, it is related to string- inspired 
dilaton gravity, which was proposed as an alternative for dark energy [55j-|57|. Then, when the limit e — > can be 
obtained smoothly, the corresponding f(Q) theory has no ghost and could actually be stable. It may be of interest to 
study the cosmology of such a theory in the limit e — > 0. 

Thus, we investigated the modified Gauss-Bonnet gravity and demonstrated that it may naturally lead to the unified 
cosmic history. Further study of cosmological properties of such a theory as well as discussion of particular models of 
the modified Gauss-Bonnet gravity may be found in refs. [58j . 
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2. F(G,R) gravity 

More general models include the dependence from curvature as well as from the Gauss-Bonnet term [54| : 

S = J d*Xy/=g (F(Q, R) + Anattcr) ■ (H.164) 



An example is given by 



which gives a solution describing the FRW universe: 

^0 , ^~ 2 fo ± V 8 /o(/o- 8^) 

H=—, h = 6 - , ■ (11.166) 

Then, for example, if n 2 fo < —3, there is a solution describing a phantom with ft,Q < — 1 — \/2 and a solution describing 
the effective matter with ho > — 1 + \/2. Late-time cosmology in other versions of such theory can be constructed. 
The unification with inflation in such model may be achieved via cosmological reconstruction, which is developed in 
the third chapter. 

C. String-inspired model and scalar-Einstein-Gauss-Bonnet gravity 

In string theories, the compactification from higher dimensions to four dimensions induces many scalar fields, such 
as moduli or dilaton fields. These scalars couple with curvature invariants. Neglecting the moduli fields associated 
with the radii of the internal space, we may consider the following action of the low-energy effective string theories 
MM: 



S = I d A XyJ—g 



R 

- + C, P + C C + 



(11.167) 



where </> is the dilaton field, which is related to the string coupling, C^ is the Lagrangian of <$>, and C c expresses the 
string curvature correction terms to the Einstein-Hilbert action, 

£0 - -d^4> ~ V(<f>) , Cc = c 1 a'e 2 ^C i c 1) + c 2 a ,2 e 4 ^C c 2) + c 3 c/ 3 e 6 ^ Cf ] , (11.168) 

where I /a' is the string tension, C c , C c 2 \ and C c ^ express the leading-order (Gauss-Bonnet term Q in (III.124[) ). the 
second-order, and the th 
have the following form 



second-order, and the third-order curvature corrections, respectively. The terms C C X \ C c 2 ^ and C c ^ in the Lagrangian 



£W=fi 2 , C^=2il 3 +R^Rf p R^ v , C^ =C 31 -5 H C 32 - 5 -^-C 33 . (11.169) 



Here, 5 b and 5h take the value of or 1 and 

n 2 = g, 

„ ^ c ^uparr, c 7 

^pv ^pa TTj 



c 32 = -(R^apR^) + - A R^RJ a Rp^RS-^Rp^RS R ^R P vlS -Rp^K p a v R P fR" a 



2 1 

~^ ±L pv * L *jS ±L pa * L af3 ML a/3 ±L o~f8 ±L p 1L pi> ML a/3 ±L pu * L ~f6 

C 33 = (R^R^ a0 f - lOR^ aP R^R aiSp R 07Sp - Rp, a pR^ZR PayS Rs 7 p ■ (H-170) 



The correction terms are different depending on the type of string theory; the dependence is encoded in the curvature 
invariants and in the coefficients (ci, c 2 , c 3 ) and 5h, 5b, as follows, 
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• For the Type II superstring theory: (ci, c%, c 3 ) = (0, 0, 1/8) and 5h = 5b = 0. 

• For the heterotic superstring theory: (01,02,03) = (1/8,0,1/8) and 8h = 1,#b = 0. 

• For the bosonic superstring theory: (01,02,03) = (1/4. 1/48, 1/8) and 5h =0,Sb = 1. 

Motivated by the string considerations, we consider the scalar-Einstein-Gauss-Bonnet gravity 2 based on [55l l60l]. 
It was first proposed in ref. [55[ to consider such a theory as a gravitational alternative for dark energy, so-called 
Gauss-Bonnet dark energy. 

The starting action is: 



S = d x^f—g 



R 1 

2^2 _ 2 



(11.171) 



Here, we do not restrict the forms of V(4>) and £(</>) which should be given by the non-perturbative string theory 
(III. 168|) . Note also that the action (|II.171|) is given by adding the kinetic term for the scalar field 4> to the action of 
the F (£/) gravity in the scalar-tensor form that appeared in the previous section. 

Making the variation of the action (|II.171[) over the metric g^, the field equations follow as: 



= -L (-sr + \g^B\ + \d»<t>d v <t> - w + \fT (-v&) + mo) 

-2£(<j>)RR' a ' - i^B^B"" - 2^{<P)R^ T R v prJT + ^{<P)R^ a R pa 

+2 (VV^)) R - 2g^ (V 2 £(</>)) R - 4 (VpV£(^)) B vp - 4 (V p V"f (<£)) R w 

+4 (V 2 £(0)) Rf v + 4g^ (V p V CT £(0)) R pa + 4 (V p V CT £(</>)) fl w "" r ■ (11.172) 

In Eq. (|II.172[) . the derivatives of curvature such as VR, do not appear. Therefore, the derivatives higher than two do 
not appear, which can be contrasted with a general aR 2 + fiR^R^ + ^R^poR^ 9 " gravity, where fourth derivatives 
of g^y appear. For the classical theory, if we specify the values of g^ v and g^ v on a spatial surface as an initial 
condition, the time development is uniquely determined. This situation is similar to the case in classical mechanics, 
in which one only needs to specify the values of position and velocity of particle as initial conditions. In general 
aR 2 + fiR^R^" + jR^paB^P' 7 gravity, we need to specify the values of g\ lL , and 5 pt , in addition to <? M1/ so that 
a unique time development will follow. In Einstein gravity, only the specific value of g^, gp U , should be given as an 
initial condition. Thus, the scalar-Gauss-Bonnet gravity is a natural extension of the Einstein gravity. 
In the FRW universe (|IL6|) . Eq. (|II.172|) becomes the following: 

= ~H 2 + + V{4>) + 24ff 3 ^ ( f W) , (11.173) 
k 2 dt 

= -1 (2H + 3H 2 ) + - VQ) - 8i7 2 ^» - IGHH^l - 16H^ im . (11.174) 
k V / 2 dt dt dt 

On the other hand, by the variation of the action (|II.171[) with respect to the scalar field, the scalar equation of motion 
follows as 

= 4> + 3Hcj) + V'((j>) +£(<f>)G. (11.175) 
In particular when we consider the following string- inspired model [55j . 

V = V e~^, £((/)) =£ Q e^, (11.176) 

the de Sitter space solution follows: 



2 f0 



e 



H 2 = H 2 ^-—^, = ^ o . (11.177) 



2 For pioneering work on the scalar-Einstein-Gauss-Bonnet gravity, see [59( . 
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Here, ipo is an arbitrary constant. If ipo is chosen to be larger, the Hubble rate H = Hq becomes smaller. Then, if 
£o ~ by choosing tpo/(f>Q ~ 140, the value of the Hubble rate H = H is consistent with the observations. The 

model (|II.176| also has another solution: 

H = !f , <t> = 00 In JL when h > , 
H = -^, = o ln^ when/i o <0. (11.178) 

Here, fto is obtained by solving the following algebraic equations: 

= _ ^ + | + V tf - ^ , = (1 - 3ft ) <Pl + 2V 4 + (ft - 1) • (H.179) 



Eqs. (III.179[) can be rewritten as 



1 \ou2„ , ^ , ^ 2 (l-5/io) 



^ = "^TTMr° (1 ^ o)+ ^ — ^ (TL1S,,) 

48 ^° - 6 - ^ ] . (IL181) 



tf K 2 (1 + fto) 

The arbitrary value of fto can be realized by properly choosing Vq and £o- With the appropriate choice of Vq and 
£oj we can obtain a negative fto and, therefore, the effective EoS parameter (|II.4[) is less than —1, w e g < — 1, which 
corresponds to the effective phantom. In usual (canonical) scalar-tensor theory without the Gauss-Bonnet term, the 
phantom cannot be realized by the canonical scalar. 

For example, if fto = —80/3 < —1 and, therefore, w = —1.025, which is consistent with the observed value, we hnd 

2 1 / 531200 403 2 ,\ /o 1 / 9 27 2 2 

Kntl = — ^ 1<T>[\K > , -T7 = 7t l<Pt\K 

1 k V 231 154 ^° / ' ^ \ 49280 7884800 



For other choices of scalar potentials one can realize other types of dark energy universes, for instance, the effective 
quintessence. Moreover, one can propose the potentials in such a way that the unification of the inflation with dark 
energy naturally occurs. Different cosmological aspects of scalar-Gauss-Bonnet dark energy have been studied in 
refs. [6l|. An extension of these studies taking into account the Euler invariant and other higher-derivative stringy 
terms has also been done [63 [. 



D. Non-local gravity 

1. Non-local F(R) gravity 



In this section, we review non-local gravity 63 65] and present its local (scalar-tensor) formulation. The explicit 
example of such a theory, which naturally leads to the unification of inflation with late-time cosmic acceleration, is 
worked out. 

The starting action of the simplest non-local gravity is given by 



S = J cfx^—g | ^ {R (1 + fiO^R)) - 2A} + £ mattcr (Q; «?) J 



(11.183) 



Here, / is some function, □ is the d'Alembertian for the scalar field, A is a cosmological constant and Q represents 
the matter fields. The above action can be rewritten by introducing two scalar fields, rj and £, in the following form: 



d Xy/—g 
d 4 X\Z—p 



1 

2^2 

1 

2k 



{R (1 + /(//)) + £ {Dq -R)- 2A} + C n 



2 {R(l + f(r))) - d^ri -£R- 2A} + £ D 



(11.184) 



By the variation of the action 184[) over £, we obtain Drj = R or 77 = □ 1 R. Substituting the above equation into 
(|II.184j) . one re- obtains (|II.183j) . 
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The variation of (jll.1841) with respect to the metric tensor gives 

o = \g»v {r (i + m - o - d P ^ v - 2A} - iv (i + m - o 

+ \ {d^dur] + d^duZ) - (g^D - V„V„) (/(r?) - £) + k 2 T> . (11.185) 

On the other hand, the variation with respect to r\ gives 

= n£ + f'(ri)R. (11.186) 

Now, we assume that the spatially fiat FRW metric and the scalar fields r\ and £ only depend on time. Then, 
Eq. (|II.185|) has the following form: 

= -3H 2 (1 + f( V ) -0 + \i"h - ZH (/' (77)17 - t) + A + « 2 P mattcr , (11.187) 

= (2iJ + 3ff 2 ) (1 + /(„) - + \iv + + 2 ^|) (/fa) - A + « 2 p . (11.188) 

On the other hand, the scalar equations are: 

= ij + 3Hf) + 6H + 12H 2 , (11.189) 

= £ + 3Hi - Uh + 12ff 2 ) f'(rj) . (11.190) 



For the de Sitter solution H = Hq, Eq. 189[) can be solved as 

77 = -4flo* - VoC- 3Hot + Vi , (H.191) 
with constants of integration 770 and 771. For simplicity, let 770 = 771 = and f(j]) be given by 

V iH O t 

f(v) = fo^ = foe~— - (11.192) 
Here, j3 is a constant. Then, Eq. 190[) can be solved as follows: 

+ (n , 93) 

Here, £0 and £1 are constants. For the de Sitter space a behaves as a — a$e Hot . For matter with a constant equation 
of state 7(;, we find 

Pmattcr = p e^ W + ^ Hot . (11.194) 

Then, after setting £ = 0, by substituting (|II.191|) . pi.l93|) . and (|II.194p into (|II.187p . one obtains 

= -3H 2 (1 + £i) + QH 2 a f Q - e" 4 ^ + A + n 2 p^ w+1 ^ t . (11.195) 



When po = 0, with the choice 

A 

the de Sitter space can be a solution. Even if /o ma tter ^ 0, if we choose 



= 2, £i=-l + _5, (IL196) 



4 K 2 p A 

P= T7TT x ■ /o= ^77271 "o \ , 6 = -1 + ^772, H.197 

3(1 + t«) 3ii5 (1 + Sw) 3Hq 



there is a de Sitter solution. 

Eq. (III.196|) or pi.!97|) shows that 



3A 

H 2 =—-. (11.198) 

1 +?i 
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This indicates that the cosmological constant A is effectively screened by £. We should also note that in the case 
without the cosmological constant A = 0, if we choose £1 = — 1, H can be arbitrary. Then, H can be determined by 
the initial condition. Because Hq can be small or large, the theory with the function (|II.192|) with /3 — 2 in (|II.19G|) 
can describe the early-time inflation or current cosmic acceleration. 

In the presence of matter with w ^ 0, the de Sitter solution H = Hq is realized even if f(rj) given by 



/W-/oe" /2 +/ie 3 ^ +1 )"/ 4 

The following solution exists: 

-AH t , £ = 1 + 3f e- 2Hot + A e -3(-+i)ffot 



V 



Pmattcr 



3(3^ + 1)^/! m+w)Hot 



(11.199) 



(11.200) 



The theory with multiple de Sitter solutions may describe the unification of inflation with dark energy [64|. Other 
cosmological aspects of such non-local gravity were studied in refs. [66j . 

One can construct more general non-local gravity along the above line. Let us start from a general non-local action: 



S = 



-g [F (R, OR, D 2 R, • • • , D m R, D^R, U~ 2 R, • • • , U~ n R) + £ mattor ] 



(11.201) 



with m and n being the positive integers. By introducing scalar fields A, B, Xki ?7fc, (k = 1, 2, • • • , to), and £j, <pi, 
(I = 1, 2, • • ■ , n), the action (|II.201I) can be rewritten in the following form: 



d X\J — g 



BR-BA + F(A,rn,ri2,-" , r? ro , <f>i, <fo, ■ ■ ■ > 4>n) + d^A + ^d^Xhd^Vh-i 



k=2 



1=1 k=l 1=2 

The variations over A, B, Xki an d £i lead to the following equations 

= R - A = m - DA = r lk - Ur} k -i =nfa-A= Ufa - fa-i {k = 2, 



, TO , I 



which give rjk = D k R and □ 4>l — R- Thus, the actions (|II.202[) and pi.201[) are equivalent. 
Let us now consider a scale transformation given by 



(11.202) 



(11.203) 



9p.ii ~ ' 2j^9p.u ■ 

Then, the action (1II.202[) is transformed into the action in the Einstein frame: 



S = I dfixyj—g 



/ m n 

2 ~ 2& d » Bd " B + 2B + 2 ^Xkd^m-t + £ W>< 

\ k=2 1=1 



1 

IP 



-BA + F(A,r]i,--- , 77^, 0i, • • • ,4> n ) + J2xiVi +Mi +J2&<t>i-i ) +£ 

\ k=l 



A 

matter 



1=2 



(11.204) 



(11.205) 



Here, £ mattcr can be obtained by scale transformation of the metric tensor in £ ma ttcr by pi. 2041) . The above action 
can be regarded as a non-linear a model coupled to gravity. The action pi.!83[) can be considered as a particular case 
of the general non-local action. In this case, however, the continuity equation pi,144[) is no longer valid; additional 
terms on the right-hand side of this equation are generated and are responsible for inducing the interaction of matter 
within the background. This aspect is crucial for the discussion of local gravity constraints for such models based on 
the modified theories of gravity. 



2. Non-local Gauss-Bonnet gravity 



It may be interesting to discuss the extension of the model presented in the above section by including the Gauss- 
Bonnet invariant [63, [68| . 
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The corresponding non-local action generalizing the modified Gauss-Bonnet gravity is 



s = J d^xyf—g f A _ ^gn-ig + c^uer) ■ (n.206) 

Here, £ ma ttcr is the matter Lagrangian, and g is the Gauss-Bonnet invariant. By introducing the scalar field <f>, one 
may rewrite the action (III.206|) in a local form: 

5 = / (2^ ~ W + <t>G + £roatter ) ' (IL207) 

In fact, from the <j> equation, it follows that <f> = — ! ^'Or 1 g. By substituting this expression into 207[) . we re-obtain 
(III.206|) . It is remarkable that the above action belongs to a specific class of string- inspired, scalar-Gauss-Bonnet 
gra vity that is known to be consistent with the local tests and may serve as a realistic candidate for dark energy 
|55j (see section C of this chapter). Unlike the situation in standard stringy compactifications, the scalar-Gauss- 
Bonnet coupling is not an exponential scalar function. The scalar potential is zero, as it should be in the low-energy, 
perturbative approach to the string effective action. 

Let us now study the de Sitter solutions that may describe late-time universe acceleration in such non-local gravity. 
Considering the FRW metric with a flat spatial part and a <j> that only depends on the cosmological time, the equations 
of motion are 

= - A" + 4? - 24</># 3 + p mattcr , = - 4 U + 3ff $) + 24 (HH 2 + H 4 ) . (11.208) 

When the matter contribution is neglected, assuming <f> and H are constant <f> = c^, H — Ho, one can solve Eqs. (|II.208|) 

as follows: = ^-£/q, Hq = — 16 3 n " 4 . Therefore, if a is allowed to be negative, there is a de Sitter universe solution, 
namely 

H = H «={-4tT ■ (il2o9) 

On the other hand, a negative a means that the scalar field has a 'wrong' kinetic term, like a ghost or phantom. As 
we will see later in pi.218[) . the de Sitter space solution corresponds to the case a — > oo (late-time acceleration) when 
matter is included. Even in the early universe, before matter was created, this solution can be used for the description 
of the inflationary epoch. 

Note that the second equation (jll.208j) yields 

+ 8R2 ) ° 3 = C ■ (11.210) 

Here, C is a constant. The action f|II.207|> is invariant under a constant shift of the scalar field (j>: cf> — > <\> + c (c is a 
constant) because g is a total derivative. The quantity C is the Noether charge corresponding to the invariance. By 
using (|II.210p . we may delete 4> in the first equation in (|II.208[) and obtain 

3H 2 160k 2 H 6 Wk 2 CH 3 k 2 C 2 

0= - + ; + fT+Pmatter. (11.211 

k z a aa* aa° 

When pmattcr is given by the sum of the contributions from all forms of matter witha constant EoS parameter Wi , one 
gets Pmatter(a) = Yli Poi a ~ d ' <yWi+1 ^ ■ Then, we can solve algebraically Eq. (jll.2111) with respect to H as a function of a: 
H = H(a). Using the definition of H, the equation has the form a = aH(a). Thus, t = J which gives, at least 

formally, the solution of pi.211[) and, therefore, the FRW universe. 

As a simple case, one may consider C = 0. Then, Eq. (|II.211[) reduces to the following form 

= Fix) = + -^- A x - " PmattCT 9 (a) . (11.212) 
V ' 160k 4 160k 2 V ' 

Here, x = H 2 and, therefore, x > 0. Because ^'(x) = 3 (x 2 + 16 q 4 ) , when a > 0, T{x) is a monotonically increasing 
function of x. On the other hand, when a < 0, we find 

^(x ± )=0, ^(^^(-^f 2 -^, I±s± (__y" J . ( „. 213) 
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Because 



when a > 0, because .F(O) < and F(x) is a monotonically increasing function, there is only one positive solution in 
(1TL2T21 : 



We now consider the case that p ma ttcr comes from a unique kind of matter with a constant EoS parameter w. Then, 
when a -> 0, one obtains ai -> ,S? 2 5 ra~ 3(1+,u) and a_ -> 0. Therefore, H cx a -(i+^')/2 an d a K ^/(i+u-) w hi cn 

' + 160^^ ' 

could be compared with the usual FRW universe in Einstein gravity, where a cx t 2 / 3 ( 1+w \ On the other hand, when 
a — > oo, Eq. (|II.215[) implies that H cx a~ 3 ( 1+w ^ 2 and, therefore, a cx t 2 / 3 ( 1 + u ') ; which reproduces the result in general 
relativity. This indicates that such a non-local term becomes less dominant as compared with the Einstein-Hilbert 
term in the late universe. 

We may consider the case when a < 0. Because J-(0) > in this case, if T (x+) > 0, there is no positive solution 
for x. On the other hand, if F (x+) < 0, that is, when 

-a>40 K 6 pLtter, (11.216) 

there are two positive solutions and one negative solution in (III.212[) , namely 



Eq. (|II.216|) indicates that there is a lower bound in a as a > a mm = (— 4Qk4 a J . Here, it is assumed that 

w > — 1. Note that even if a — a m - m , x and, therefore, H do not vanish, which means that the existence of a m i n does 
not imply the bounce solution. It is likely that a universe with a — a mm could appear, for example, as a quantum 
correction. 



When a is large, we find that J-(0) in (|II.214|) vanishes and a± tends to a constant: a\ — > zti ^^ y — fjj, which 
gives 



H 2 =x = 0, ±J ^-r. (11.218) 

' V 160k 4 v ' 

The positive solution of x in (|II.218|) corresponds to the de Sitter solution (|II.209|) . 
When the Einstein-Hilbert term in (|II.206j) or (|II.207P can be neglected, by assuming 

<t>=^, H=y (a = a t h °) , (11.219) 

Eqs. (|II.208j) reduce to the following algebraic equations 

= ^ (^ a <f> + 24h o) > = (ho - 1) (-^ + 4h£) . (11.220) 

The equations in (|II.220[) have a non-trivial solution (a^, ^ and ho ^ 0): ho — 1 and = — ^j^- 
Matter can also be accounted in the model; in that case, the first equation in pi.220|) is modified as 

= + ^j 2 - 2^H 3 + p mattor . (11.221) 

When one may neg lect the Einstein-Hilbert term in (III.206|) or (|II.207p . and if we assume pi. 2191) , there is a non-trivial 
solution when ho = g7~jrn ■ I n fact; Eqs. (|II.208j) reduce to the following algebraic expressions 

= J al + + 48«» (5^)' , = (^ - l) + 4 (j^)") . (H.222) 
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If the matter is radiation with w = 1/3, the second equation is satisfied and the first equation acquires the following 
form: = ^a| + p o a o " 4 + 48a . Then, if D = (24) 2 - > 0, there are non-trivial solutions: a = ^ (-24 ± y/o). 

On the other hand, when w ^ 1/3, we find = — ( 3(1^) ! > /°o a o~ 4 = 16 ° K " ( 3(1+^) ) • The secon d equation 
requires a to be positive. Thus, a phantom like universe solution from purely non-local gravity, without the R term, 
has been found. In the same way, in the case where several de Sitter solutions exist one can expect the unified 
description of inflation with dark energy. 



3. Other classical non-local Gauss-Bonnet models and their de Sitter solutions 



Another model giving rise to an interesting non-local action is the following 

R k 2 



S 



d Xy/—g 



2k 2 2a 



F{Q)U^ 1 F(Q) 



(11.223) 



where a is a dimensional constant and F(Q) is an adequate function of Q. By introducing three scalar fields, <j), £, 
and 77, one can rewrite the action (|II.223[) in the following form: 



S 



^ 2^ + 2^ 



+ <l>F{v) +Uv- G) 



We may further add a potential V(4>) to the action 



R a 



S = / <?Xy/=g — + —,d u 4>d^ + <t,F(r,) +t(v-G)- V(cf>) 



2k 2 2k 2 

In the FRW universe, this action leads to the following equations 



(11.224) 



(11.225) 



,^4> 2 + 2^H 3 + V(<f>) , 0=^U + 3H<j>)+F( V )-V^), 
= <j)F{ii) + £ , = 77 - 24 (h 4 + HH 2 ^ 



°=-- H -2 



(11.226) 



However, after adding the potential V (</)), it is difficult to obtain the corresponding non-local action explicitly. 

When V(4>) = 0, assuming that </> = cpt, r\ = tjq, £ — C£b, and H = Hq, with constant cp, 770, eg, and Hq, the 
equations in (jll.2261) reduce to the algebraic equations 



3 rJ 2 a „2 1 o/i„ zj3 - ^a 



° = -^2 H o - ^2 4 + 24c ?^o > = 72 Hocp + F( m ) , = c^F\n Q ) + c t , = m - 2AH 4 . (11.227) 



One can solve Eqs. (|II.22T[) with respect to c^, eg, and t]q a s follows: = — 

~^- F (24H 4 ) F' (2AH*), m = 24iJ 4 , and we find 



*F(24Hg 
3 a Ho 







Here, G - 24ffjf. 
For example, if 



Eq. pi.2286 gives 



8«r 18a 9a 



F(G) = f G 2 , 



°^ Go '-8^ 



1 , H« 2 /o 2 G 3 



18a 



which has a trivial solution Go = corresponding to the flat background and 

3 _ 9a 
° " U^f 2 ' 



c c = 



(11.228) 



(11.229) 



(11.230) 



(11.231) 
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which corresponds to the de Sitter universe. 
As another example, one can choose 

F {G? = -g {Q + gi){G-92) ■ (11.232) 

Here, go, gi, and gi are positive constants. For F (Q) to be real, the value of Q is restricted as —g\ < Q < gi- Then, 
Eq. (|II228I) yields 



= G o + 7^2 G ° + -9152, (H.233) 
4g n z 

which can be solved as 



= <"*'" 4 : * ^ — ] 4,/ " /2 l ' (IL - :i4 ' 




The above solution is positive, provided that a < 0. For —g\ < Gq± < g2, we require 

9a 9a 
4^' 92 > "8^o 



.9i>-7U^' ff2>-^-. (IL235) 



Then, as long as Eq. (|II.235[> is satisfied, there are two solutions describing the de Sitter universe. By choosing a more 
general F(Q), we may find multiple de Sitter universe solutions, among which the solution with the largest value of 
Hq could perfectly correspond to the inflation epoch and the one with the smallest value of Hq to the dark energy 
era. Note that, in the model (|II.232[) . as long as Eq. pL235|) is satisfied, there is no singularity, for example at Q = 0. 
Thus, a smooth transition could occur and there exists a true possibility to realize matter dominance before late-time 
acceleration. The important lesson from the above investigation is again the fact that the inflation and dark energy 
epochs may be natural solutions of the field equations. 

E. Gravity non-minimally coupled with the matter Lagrangian 

1. Non-minimally coupled scalar theory 

It could be interesting to investigate the role of the non-standard non-minimal coupling of modified gravity with 
the matter Lag rangian in the cosmic acceleration epoch. The corresponding, quite general theory has been suggested 
in refs. (69l. l70l| (see also [7l|). As an example of such a theory, the following action is considered: 

,4 . f R 



S = j d^x^g | — + / (R) L d | . (11.236) 

Here, Ld is a standard matter like action that produces dark energy due to non-trivial coupling with the curvature. 
In principle, the first term may be more general, for instance, of the f(R) form. However, for illustrative purposes it 
is enough to consider this simpler variant of the non-minimal theory. The second term in the above action describes 
the non-linear coupling of matter with gravity. Very often, such terms may be induced by quantum effects such as 
the non-linear effective action. Thus, it is natural to consider that the second term belongs to matter sector. 
By the variation over g^ v , the equation of motion follows as: 

1 5S 1 fl 

^gSg^u 2k 2 [2* 



= ^=-^— = —\- g t»'R-Ri iV \+T tu '. (11.237) 



Here, the effective energy-momentum tensor T M „ is defined by 



f ee -f'(R)R^L d + (V^ - g^V 2 ) (f'(R)L d ) + f(R)T^ , 



g bg^v 

As a matter we consider the massless scalar 



T"" ee —L__*_ ( / dSV^Ld) . (11.238) 



L d = ~g^d^<f>. (11.239) 
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Then, the equation given by the variation over <f> has the following form: 



= -L^ = -La M {f(R)V^g^d^) . (11.240) 



It is enough to consider the case that 



/(*)=( 4) ■ ( IL2 41) 



We again assume the spatially flat FRW universe (|II.6[) and that the scalar </> only depends on t ((f) = <f>(t)). The 
solution of the scalar field equation pi.240p is given by 

4> = qa- 3 R- a . (11.242) 

2 

Here, q is a constant of the integration. Thus, R a Ld = 2 J>r^ , which becomes dominant when R is small (large) 
compared with the Einstein term if a > —1 (a < —1). Thus, one arrives at the remarkable possibility (69l. |70T| 
that dark energy grows to asymptotic dominance over the usual matter with a decrease of the curvature. In the 
current universe, this solves the coincidence problem (the equality of the energy density for dark energy and for 
matter) simply by the fact of universe expansion. 

Substituting (jll.2421) into 23T[) . the (/i, v) — (t, t) component of the equation of motion has the following form: 

3 u2 36g 2 (a(a + l)~ a + l- 2 



l + —a + a 2 )HH 2 +(l + -a)H i \. (11.243) 



4 / V 2 

The accelerated FRW solution of (III.243|) exists as [gl, |7(J 



a + l\ 6 _ 2n 2 q 2 (2a-l)(a-l) 



St 



\H = —\, og= — 1 — T2 - (IL244) 



M 2Q 3(a + l)" +i (f (2a-l)) 



Eq. (|II.244[) shows that the universe accelerates, that is, a > if a > 2. If a < — 1, the solution (1II.244[) describes 
a shrinking universe if t > 0. If the time is shifted as t — > t — t s (t s is a constant), the accelerating and expanding 
universe occurs when t < t s . In the solution with a < -1 a Big Rip singularity appears at t = t s . For a oc t h °, the 
effective EoS parameter w c s in pi.21j) is given by w c g = — 1 + and the accelerating universe expansion (ho > 1) 

occurs if — 1 < w e s < — |. For the case of (|II.244|) . one finds 

w cS = . (11.245) 

1 + a 

Then, if a < —1, we have w e s < — 1, which is an effective phantom regime. For general matter with the constant EoS 
parameter w, the energy E and the energy density p m attcr behave as E ~ a~ 3w and p ma ttcr ~ a~ 3 ( w+1 \ Thus, for the 
standard phantom with w < —1, the density becomes large with time and might generate the Big Rip. Thus, it is 
demonstrated that the non-linear coupling of gravity with matter may produce an accelerated cosmological expansion 
with the effective EoS parameter close (above, equal or below) to —1. More complicated couplings [69|, [70, l72( may 
be considered in the same fashion. Note that such theories show some very interesting properties, like the induction 
of extra gravitational force (73l |. Other cosmological aspects of such non-minimally coupled non- linear theories are 
discussed in refs. 17411. 



2. Non-minimally coupled vector model 

The non-minimal coupling of gravity with the matter Lagrangian of the previous section may be considered for 
more complicated models. Let us investigate the abelian vector field Lagrangian non-minimally coupled with gravity 
75]: 

S=f d'x^l^ ^F^FA . (11.246) 



3G 



Here, I(R) is an appropriate function of the scalar curvature R, and is the curvature of the f (1) gauge field. 
Making the variation of the action Eq. (|II.246[) with respect to the metric g^ v and the U(l) gauge field A a, we obtain 
the following equations: 



R - -a R - k- 2 t( em ) 

T^ M) - (g aP F^F va - \g,„F aP F^ 



~ {f(R)F a0 F^R^+g^n [f'(R)F a pF a P] - V M V„ [f(R)F a pF^] } , 



and 



--()„ 



^I(R)F^) =0. 



(11.247) 

(11.248) 
(11.249) 



In the FRW universe (|II.6[) . the (jj,,v) = (0,0) component and (//, v) = component of Eq. (III.247[) gives 

.2 r / i 



H 



^(i?) ^^Jb^Oa + ^i^^ 

3 
2 



-/'(#) (tf + # 2 ) + 6f"(R)H (h + 4HH) 



F a f}F° 



(11.250) 



2H + 3H 2 = 



k ri 

6 



I(R)F a0 F 



a/3 



f'(R) [H + 3H' 



6f"{R) H + 7HH + AH 2 + 12# 2 # + 36 f"'{R) H + AHH 



FapF 



a0 



f'(R)H + 4/"(i?) (h + <LHH\ 



lf'(R)^A(F aP F^) 



(11.251) 



By properly choosing I(R), a model unifying the late acceleration and the inflation emerges. However, the properties 
of the model may be investigated only qualitatively 75] or numerically because the equations of motion are very 
complicated. Alternatively, as will be shown in the next chapter the cosmological reconstruction of the theory under 
investigation may be conducted so that unified cosmic history is realized. Different properties of the non-minimally 
coupled vector theory have been studied in refs. [76| . 



3. Non-minimal Yang-Mills theory 

Here, we consider the cosmology in the non-abelian non-minimal vector- -F^i?) gravity based on ref. [77]. Different 
cosmological proposals regarding the non- minimal Yang-Mills theory may be found in refs. [78j |. 
The proposed action follows: 



Smg ■ 



where F^ v is given by 



d 4 Xy/^g [Cmg + £v] , 



±[R + F(R)} , £ V = I(R) {-\F; v F a ^ - V[A 2 ] 



f± = a,>Ai - d v Ai + r hc K A i 



A 2 = A a1 



(11.252) 



(11.253) 



and A 2 = g^A^A®. We should note that the last term V[A 2 ] in the action (|II.252[) is not gauge invariant, but it can 
be rewritten in a gauge-invariant way. For example, if the gauge group is a unitary group, one may introduce a a 
model like field U, which satisfies Wll = 1. Then, the last term could be rewritten in the gauge-invariant form: 



V[A 2 } -> V [ctr (tfA'U) (U^A^U)] 



(11.254) 
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Here, c is the normalization constant. Choosing the unitary gauge U = 1, the term in (jll.2541) reduces to the original 
one: V^[A 2 ]. This may indicate that the action (|II.252|) describes the theory in which the gauge group is spontaneously 
broken. 

The field equations can be derived by taking variations of the action pi.252[) with respect to the metric g^ v and 



the vector field Af : as follows: 



with 



and 



[1 + F'(R)} - X - 9llv [R + F(R)} + 9llu UF'{R) ~ V M V„F'(i?) = n 2 T^ , (11.255) 



<a Aa dV{A 2 } 1 



+ i {f'(R)FR^ + g^D [f'(R)f] - V M V„ [f'{R)P\ } , (11.256) 
T = F"F a ^ + W[A 2 }, (11.257) 



=fy [V^gI(R)F a ^] - I(R) \ ^A\F C ^ V + 2 dV }^ A av \ = , (11.258) 



y/=g " LV * K 1 J y ' y M dA 2 

where T^) is the contribution to the energy- momentum tensor from A^. As will be shown in the next chapter, the 
cosmological reconstruction of the above theory may result in the accelerating expanding universe. The direct solution 
of the equations of motion is extremely complicated. 



F. Modified F(R) Hofava-Lifshitz gravity 

1. General formulation 

This section is devoted to a review of the F(R) Hofava-Lifshitz gravity [7{|. The FRW equations for this theory 
are also formulated. The Hofava-Lifshitz version of general relativity was originally introduced as a candidate for 
renormalizable quantum gravity at the price of the violation of Lorentz invariance. The below model represents 
natural generalization of the standard F(R) gravity (III.ll) discussed at the beginning of this chapter. 

Let us present the general formulation of the theory. By using the ADM decomposition [8(J Ell (f° r reviews, see 
[HI HI]); one can write the metric in the following form: 

ds 2 = -N 2 dt 2 + g<lf (dx l + ISTdt) (dx j + N J dt) , I = 1, 2, 3 . (11.259) 

Here, N is called the lapse variable and N l, s are the shift variables. The scalar curvature R becomes: 

R = K l: >K tJ - K 2 + i? (3) + 2V M {n^ v n v - n v \s ' v n») , (11.260) 

and \f—g = \J g^N. Here, R^ is the three-dimensional scalar curvature defined by the metric g^\ and l£y is the 
extrinsic curvature defined by 

Kij = ~ " V ^ Nl ) ' K = K *i- ( IL261 ) 

(3) 

is a unit vector perpendicular to the three-dimensional hypersurface defined by t = constant and V i expresses 
the covariant derivative on the hypersurface T, t . 

The first attempt of the extension of F(R) gravity to a Hofava-Lifshitz-type theory [84| which is known to be a 
power-counting renormalizable model was proposed in ref. [85j j . The starting action is 

S Fhl(r) = J d^xy/g^NF^) , i? HL = K ll K l3 - XK 2 - E l ' 3 Q l]kl E kl . (11.262) 

Here, A is a real constant in the "generalized De Witt metric" or "super-metric", 

gijkl = 1 (g(3)ik g (W + 5 (3)« 5 (3 W ^ _ Xg m g (3)kl f (IL263) 
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defined on the three-dimensional hypersurface £ t , and E lJ can be defined by the so-called "detailed balance condition" 
by using an action W[<jrj. ; ] on the hypersurface E t 

V^W J ' = SW } 9kl ] , (11.264) 



and the inverse of Q^ kl is written as 

- 2 \ 9 ik 9ji + 9 a 9 jk ) ~ x 9ij 9 M , A - gA _ - . (11.265 J 

The action Wfg^ ] is assumed to be defined by the metric and the covariant derivatives on the hypersurface £ t . The 
original motivation for the detailed balance condition is its ability to simplify the quantum behavior and renormaliza- 
tion properties of theories that satisfy the detailed balance condition. Otherwise, there is no a priori physical reason 
to restrict E ij to be defined by pi.264[) . 

There exists an anisotropy between space and time in the Hof ava-Lifshitz approach. In the ultraviolet (high-energy) 
region, the time coordinate and the spatial coordinates are assumed to behave as 

x^bx, t^b z t, z = 2,3,--, (11.266) 
under the scale transformation. In ref. [84} . W A [gj^' > ] is explicitly given for the case z = 2, 

W = -3- [ d3x (R - 2A W ) , (11.267) 
K w J 

and for the case z = 3, 

W= \ [ uj 3 (T) . (11.268) 

w JSt 

Here, % is a coupling constant of dimension —1/2 and w 2 is the dimensionless coupling constant. 1^3 (F) is given by 

ws(r) = Tr (V a dr + |r a r a rj = e^ k (r%d 3 r{ m + |rsrj. m rj^) Sx . (n.269) 

A general E v consists of all contributions to W up to the chosen value z. 

In the Hofava-Lifshitz like F(R) gravity, it is assumed that N can only depend on the time coordinate t, which 
is called the "projectability condition". The reason for this is that Hofava-Lifshitz gravity does not have the full 
diffeomorphism invariance, but rather is invariant only under "foliation- preserving" diffeomorphisms, i.e., under the 
transformations 

8x l =C(t,x), St = f(t). (11.270) 

If N depends on the spatial coordinates, it could not be fixed to be unity (N = 1) by using the foliation-preserving 
diffeomorphisms. There exists a version of Hofava-Lifshitz gravity without the projectability condition, but it is 
suspected to possess a few additional consistency problems [86|, [8?J ■ Therefore, we prefer to assume that N depends 
only on the time coordinate t. 

Let us consider the FRW universe with a flat spatial part and the lapse function N: 

ds 2 = -N 2 dt 2 +a(t) 2 ( dxl ) 2 ■ (11.271) 

i=l,2,3 

It is clear from the explicit expressions (jll.2671) and (III.268[) that W[<?^] vanishes identically if Aw = 0, which is 
assumed because a non-vanishing Aw gives a cosmological constant. Then, one can obtain 

R _12H 2 6 d fH\ 6H 2 6 d f Ha^\ (3 - 9A) H 2 
R -— + Nd-t(N = ~~ + ^d-t(—h Rnh - N 2 • (IL2?2) 



In the case of general relativity, the second term in the last expression for R becomes a total derivative 

- 1 *• ^ < + a?l (£) W **{-« nv < 1 > " 
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Therefore, this term can be dropped in Einstein gravity. The total derivative term comes from the last term 
2V^ (n^V v n 1 ' — n"V v n^) in (|II.260p . which is dropped in the usual Hofava-Lifshitz gravity. In F(R) gravity, however, 
this term cannot be dropped due to its non-linearity. Then, if we consider FRW cosmology with the flat spatial part, 
there is almost no qualitative difference between Einstein gravity and Hofava-Lifshitz gravity, except that an effective 
dark matter could appear as a kind of integration constant in Hofava-Lifshitz gravity (88|. The effective dark mat- 
ter appears because the constraint given by the variation over N becomes global in the projectable Hofava-Lifshitz 
gravity. 

Now, we introduce very general Hofava-Lifshitz like F(R) gravity [79[ by 

S F(k) = J d 4 xy^NF(R) , R ee K ij K i:j - XK 2 + 2^V M (n"V ' v n v - n v V v n») - E lJ g ijkl E kl . (11.274) 

In the FRW universe with a fiat spatial part, R has the following form: 

- (3-9A)g 2 6 M d fHa 3 \ (3 - 9A + lgg) H 2 6/i d ( H\ 

R= ^^ + ^Ndt\—)= + Nd-t{Nj ■ (IL275) 

The case one obtains with the choice of parameters A = fj, = 1 corresponds to the usual F(R) gravity as long as 
we consider spatially flat FRW cosmology, because R reduces to R in (jll.272p . On the other hand, in the case of 
/i = 0, R reduces to i?HL pi.272[) and, therefore, the action (III.274[) becomes identical with the action (|II.262[) of the 
Hofava-Lifshitz like F(R) gravity in ref. [85]. Thus, the fi = version corresponds to some degenerate limit of the 
above general F(R) Hofava-Lifshitz gravity. We call this limit degenerate because it is very difficult to obtain the 
FRW equations when /x = is set from the very beginning. In our theory the FRW equations can be obtained quite 
easily, and then fi — is a simple limit. 

For the action (|II.274|) . the FRW equation given by the variation over has the following form after assuming 
the FRW space-time pi.271[) and setting N = 1: 



/ _\ f . \ / ~\ dF'(R) d 2 F'(R) 
= F i^Rj - 2 (1 - 3A + 3/i) [H + 3H 2 j F' i^RJ - 2 (I - 3A) H + 2^ -^-L 

On the other hand, the variation over N gives the global constraint: 



Pmatter ■ (11.276) 



= / d A x 



dF' I R 

F (jlj - 6 1 (1 - 3A + 3 M ) H 2 + F' (rJ + 6fiH ^J- - p matter 



(11.277) 



after setting N = 1. Because N only depends on t, and thus does not depend on the spatial coordinates, we only 
obtain the global constraint given by the integration. If the standard conservation law (jII.I44l) is used, Eq. (III.276[) 
can be integrated to give 



f \ r -\ r \ dF'(R) c 
= F \R) - 6 { (1 - 3A + 3^) H 2 + fiHj F' \R) + 6fiH - p mattcr - -3 . (11.278) 

Here, C is the integration constant. Using (III.277[) . one finds C = 0. In 88], however, it was claimed that C need not 
always vanish in a local region, because (|II.277|) needs to be satisfied throughout the entire universe. In the region 
C > 0, the Ca~ 3 term in (|II.278[) may be regarded as dark matter. 

Note that Eq. (|II.278[) corresponds to the first FRW equation and (|II.276|) to the second one. Specifically, if we 
choose A = M = 1 and C = 0, Eq. (|II.278[) reduces to 



dF' 



*(r) 

= F (Jlj — 6 [H 2 + H) F' + 6H — - p mattc 



= F (r) - 6 (H 2 + ij) F' (R) + 36 (4H 2 H + H^j F" (r) - p mattcr . (11.279) 

This coincides with the corresponding FRW equation for traditional F(R) gravity introduced in the first section of this 
chapter. Note that in the degenerate /i — case [85], the action (|II.274|) or f|II.262|) does not contain any term with 
second derivatives with respect to the coordinates. Such a term appears in the standard F(R) gravity. The existence 
of the second derivatives in the usual F(R) gravity induces the third and fourth derivatives in the FRW equation as in 



40 



pi,276[> . Due to such higher derivatives, an extra scalar mode appears, which is often called the scalaron in the usual 
F(R) gravity. This scalar mode often affects the correction to Newton's law as well as other solar tests. Therefore, 
such a scalar mode does not appear in the F(R) Hofava-Lifshitz gravity with fi = 0. Thus, we have formulated a 
general Hofava-Lifshitz F(R) gravity that describes the standard F(R) gravity or its non-degenerate Hofava-Lifshitz 
extension in a consistent way. 



2. FRW cosmology 

This section is devoted to the study of the FRW Eqs. (jll.2761) and (|II.277p . which admit a de Sitter universe solution. 
We follow ref. [79]. It is convenient to neglect the matter contribution: p m attcr = Pmatter = 0. Assuming H = Hq, 
both of Eqs. (|II.276|) and (|II.277|I lead to the same equation 

= F (3 (I - 3A + 6fj) Hi) — 6 (I - 3A + 3/x) H%F' (3 (I - 3A + 6/x) Hi) , (11.280) 

as long as the integration constant vanishes (C = 0) in Eq. (1II.278[) . 
First, we consider the popular case that 



F 



(it) (xR + (3R 2 . (11.281) 



Then, Eq. (|II.280p gives 



= # ( 2 {l-3A + 9/?(l-3A + 6/i)(l-3A + 2^).ff 2 } . (11.282) 



In the case of usual F(R) gravity, where A = fi = 1 and, therefore, 1 — 3A + 2fi = 0, there is formally only a trivial 
solution Hi = 0. For our general case, however, there exists the non-trivial solution 

H ° = ~ / 3(1-3A + 6 M )(1-3A + 2 M ) ' (IL283) 

as long as the right-hand side of (|II.283[) is positive. If the magnitude of this non-trivial solution is small enough, this 
solution might correspond to the late-time accelerating expansion. Thus, the R 2 term may generate the late-time 
acceleration. On the other hand, the above solution may serve as an inflationary solution for the early universe (with 
the corresponding choice of parameters). 

Instead of (|II.281|) one may consider the following model: 



F (#J oc R + PR 1 + jR 6 . (11.284) 

Eq. (|II.280|) becomes 

= Hi |l - 3A + 9/3 (1 - 3A + 6/i) (1 - 3A + 2/i) Hi + 9 7 (1 - 3A + 6fif (5 - 15A + 12 n) H^ , (11.285) 
which has the following two non-trivial solutions, 



Hi = -- (l-3* + W i± A _4(I~3A)(5-I5A+12,) 7 

2(1-3A + 6 M )(5-15A+12 M ) 7 I V 9 (1 - 3A + 2^) 2 /3 2 )' 1 ' 



as long as the r.h.s. is real and positive. If 

4(1-3A) (5 - 15A + 12/i) 7 



9 (1 - 3A + 2/i) /3 2 



< 1 , (11.287) 



one of the two solutions is much smaller than the other solution. Then, one may regard that the larger solution 
corresponds to the inflation and the smaller one to the late-time acceleration, similarly to the modified gravity model 
0, where such unification had been first proposed. Note that some of the above models may possess the future 
singularity in the same way as the usual F(R) gravity. However, it would be possible to demonstrate that adding 
terms with even higher derivatives might cure this singularity, as the addition of the R 2 term did in the standard 
F(R) gravity. Thus, we demonstrated the qualitative possibility to unify the early-time inflation with the late-time 
acceleration in the modified Hofava-Lifshitz F(R) gravity. 
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3. Unified inflation and dark energy in modified Hofava-Lifshitz gravity 

Let us consider here some viable F(R) gravities that admit the unification Q of inflation with late-time acceleration. 
In the traditional F(R) theory, a number of viable models (see [30l - l32j ). which pass all local tests and are able to 
unify the inflationary and the current cosmic accelerated epochs, have been proposed. Here, following ref. [89], it is 
shown that this class of models may be extended to the Hofava-Lifshitz gravity. The starting action is 

F(R) = R + f(R) , (11.288) 

where it is assumed that the term f(R) becomes important at cosmological scales, whereas for scales comparable to 
that of the solar system, the theory becomes linear on R. As an example, the following function [3ljj is considered 

f(R) = fr(«fr-0, (II .289) 

where a, /?, and 7 are constants and n > 1. This theory reproduces the inflationary and cosmic acceleration epochs 
in convenient F(R) gravity (see ref. [31J). During inflation, it is assumed that the curvature scalar tends to infinity. 
In this case, the model (IIL288[) . with (|II.289[) . behaves as 

lim F(R) = aR n . (11.290) 

Then, by solving the FRW equation (|II.278[) with C = 0, this kind of function yields a power-law solution of the type 

hi 2u(n- l)(2n- 1) 

ff W = T' Wh6re ^ l-3A + 6,Jn(l-3A + 3M) - ^ 

This solution produces acceleration during the inflationary epoch if the parameters of the theory are properly defined. 
The acceleration parameter is given by - = hi(hi — l)/t 2 ; thus, for hi > 1, the inflationary epoch is well reproduced 
by the model (JTT7289J) . Alternatively, the function (jll.2891) has a minimum at i?o, given by 

Ro~(— ) , /'(£)= 0, /(£)=-2A~-£, (11.292) 

\OTfJ 7 

where the condition /3~//a 3> 1 is imposed. Then, at the current epoch, the scalar curvature acquires a small value 
that can be fixed to coincide with the minimum (|II.292p . such that the FRW equations (|II.276[) and (|II.278[) with 
C = yield 

u 2 K 2A • 2 Pmattor + Pmattcr /TTino\ 

H = 3(3A^T) / ' mattcr + 3(3A^T) = ■ ( IL293 ) 

which look very similar to the standard FRW equations in general relativity, except for the parameter A. (F or a first 
consideration of the FRW dynamics in Hofava-Lifshitz gravity based on general relativity, see refs. [aS-Ell-) As has 
been pointed out, at the current epoch, the scalar R is small, so the theory is in the IR limit where the parameter 
A ~ 1, and the equations approach the usual equations for F(R) gravity. Thus, the FRW equations (|II.293[) reproduce 
the behavior of the well-known ACDM model with no need to introduce a dark energy fluid to explain the current 
universe acceleration. 

As another example of the models described by pi.288[) , one can consider the function [l5|, [53, HH , 

f(R) = (J V^ o)2 " +1+ ^" +1 T = -1 + = ^ T . (11.294) 

fa + h [(R - Ro) 2n+1 + R 2 Q n+1 \ & fo + fi [(R - Ro) 2n+1 + R 2 n+1 

This function could also serve as the unification of inflation and cosmic acceleration; however, in this case, when one 
takes the limit R — > 00, one gets 

lim F(R) =R-2Ai, where Aj = l/2/i , (11.295) 

R— >-oo 

where the subscript i denotes the inflationary. By inserting this into Eqs. (|II.276[) and 278[) with C = 0, the FRW 
equations take the same form as in (jll.2931) . Then, for the function (|II.294[) . the inflationary epoch is produced by an 
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effective cosmological constant, which implies that the parameter A > 1/3, or the equations themselves will present 
inconsistencies, as was discussed in the above section. For the current epoch, it is easy to see that the function (|II.294p 
exhibits a minimum for R = Rq, which implies, as in the model above, an effective cosmological constant for late time 
that can produce the cosmic acceleration. The emergence of matter dominance before the dark energy epoch can be 
exhibited, in analogy with the case of the convenient theory. Thus, we have shown that the model pi.294[) also unifies 
the cosmic expansion history, although with different properties during the inflationary epoch as compared with the 

model (|II.289I) . " 

It is also interesting to explore the de Sitter solutions allowed by the theory (|II.262[) . By taking H(t) = Ho, the 
FRW equation (|II.278|) . in the absence of any kind of matter and with C = 0, reduces to 

= F{R ) - 6i?o (1 - 3A + 3fi)F'(R ) , (11.296) 

which represents an algebraic equation that can be solved to yield the possible de Sitter points allowed by the theory. 
As an example, let us consider the model pi.289[) . where Eq. pi.296[) takes the form 



R n {aR n -f3) 6# 2 (-l + 3A-3^) 1 + najR^' 1 + i?™' 1 (2 7 i? - np) + Rq " (7 #0 + 2na) 



0. (11.297) 



Here, Rq = 3(1 — 3A + 6^)i?g . By specifying the free parameters of the theory, one can solve Eq. (jll.2971) . which yields 
several de Sitter points, as the one studied above. The de Sitter points can be used to explain the coincidence problem, 
with the argument that the present will not be the only late-time accelerated epoch experienced by our universe. In 
standard F(R) gravity, it was found that this model contains at least two de Sitter points along the cosmic history 
(see [ll3| ). In the same way, the second model studied here (jll.294|) provides several de Sitter points in the course of 
cosmic history. Note that when \i = A — i, Eq. (|II.296|) turns out to be much more simple; it reduces to F(Rq) = 0, 
where the de Sitter points are the roots. For example, for (|II.297|) we have Rq(1 + 7-Ro) + R^aR^ — j3) = 0, where 
the number of positive roots (de Sitter points) depends on the free parameters of the theory. 

In summary, it has been demonstrated here that, as in F(R) Hofava-Lifshitz gravity, the so-called viable models, 
as in (jll.289[) or 294[) . can in fact reproduce the whole cosmological history of the universe, in the same way as 
convenient modified gravity. 



G. Covariant power-counting renormalizable gravity 

M otiva ted by the idea of ref. [84j . a model of covariant, p ower -counting renormalizable gravity was proposed 

[ll4l Ill5t. Its prototype using the perfect fluid is given in ref. [l OQt ] . Such a theory seems to be more natural than 
models [79L IHJ because Lorentz invariance is not broken explicitly. Moreover, the theory looks very similar to R 2 
gravity with a scalar. The present section is devoted to review of such covariant model in the FRW universe. 



1. Perfect fluid formulation 

Let us briefly review the covariant, power-counting renormalizable gravity of ref. |l00| . The starting action is 

(11.298) 



S = J d±xV=g^-a(T^R^+(3TR) 2 ^ 



Here, is the energy-momentum tensor of some string-inspired fluid. The action (III.298|) is fully diffcomorphism 
invariant. We consider the perturbation from the fiat background = r]^ + h^ v . The following gauge conditions 
are chosen: hu — hu = hu = 0. The fluid energy-momentum tensor in the flat background has the following form: 



T tt =p, Tij = p% = wpSij 
Here, w is the EoS parameter. Then, one finds 



(11.299) 



T^R 



}1U 



/3TR 



j-i + I + (-1 + 3w) ^| d 2 t (^7,,,) + (w-(3 + 3w(3) d'&hi 
+ {-w + !3-3wf3)d k d k (fthij) . 



(11.300) 
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If 



the second term in the action (|II.298I) becomes 



^ = -2(3^1)' ( IL301 ) 



a {T^R„ V + pTR) 2 = ap 2 ( | + \) {Wh^ - 8 k d k (Phy) }' , (11.302) 

which does not contain the derivative with respect to t and breaks the Lorentz invariance. Imagine p is almost constant. 
Then, in the ultraviolet region, in which k is large, the second term in the action (III. 2981) gives the propagator behaving 
as 1/ |fe| 4 . It renders the ultraviolet behavior (compared with Eq. (1.4) in [Hj]). Note that the form pi.301[) indicates 
that the longitudinal mode does not propagate; only the transverse mode propagates. 

The action (|II.298|) gives 2 = 2 theory. For the theory to be ultraviolet, power-counting renormalizable in 3 + 1 
dimensions, z — 3 theory is necessary. To obtain such a theory, we note that, for any scalar quantity $ with the 
choice 

7 - , (H.303) 

6W — 1 

one obtains 

T^V^V^ + 7TV"Vp$ = p (w + 1) d k d k <Z> , (11.304) 

which does not contain the derivative with respect to the time coordinate t. This is true even if the coordinate frame 
is not the local Lorentz frame. The derivative with respect to the time coordinate t is not contained in any coordinate 
frame where the perfect fluid does not flow. Thus, the proposed action looks as 

S = J d^x^f—g | A _ a (T^R^ + j3TR) (T^V M V, + tTV%) (T^R^ + /3Tii) j , (11.305) 

with P = — 2 (3wli) an d 7 = 3w—i • corresponds to z = 3 theory, which seems to be (power-counting) renormalizable. 
In general, for the case 



S = I (ftx^J—g 



- a {(T^V M V, + yTV p V p ) n (T^R„ U + (3TR)} 2 



(11.306) 



with a constant n, z = 2n + 2 theory, which is super-renormalizable for n > 1 emerges. 

The second terms in the actions (|II.298p . (jll.3051) . and (jll.3061) . which effectively break the Lorentz symmetry, are 
relevant only in the high- energy /UV region because they contain higher-derivative terms. In the IR region, these 
terms do not dominate, and the usual Einstein gravity follows as a limit. 



2. Covariant field theory of gravity with a Lagrange multiplier 

We now show that the perfect fluid, which appeared in (|II.298|) . (|II.305p . and pi.306[) . can be realized by a scalar 
field and the Lagrange multiplier field. As a result, we obtain a covariant, power-counting renormalizable field theory 
of gravity. 

The following constrained action for the scalar field <j> is proposed: 

^ — r . n 



= I (fx^g- { -X ( -d„ m + U(<t>) ) \ . (11.307) 



2 

Here, A is the Lagrange multiplier field, which gives a constraint of 

ia M <OT + £/(0) = O, (11.308) 

that is, the vector (d^cj)) is time like. At least locally, one can choose the direction of time to be parallel to (cL<^>). 
Then, Eq. pi.308|) has the following form: 

1 '-'- 2 



nl^r) ( IL309 ) 
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The equation given by the variation over <fi will be discussed later. 

The scalar field energy-momentum tensor T^ u with an arbitrary scalar potential is defined as follows: 

t% = d^dA - 9lut (^d^dy + v(<t>)) . 

The "energy density" and "pressure" p<p become: 



(11.310) 



n = \{^) -v(« = ^)-vw, p + = \(j[) +n<t>) = u(d>) + v(ct>). (n.311) 

Here, the constraint pi.309[) is used. Note that (unknown, for the moment) V(4>) is not identical with £/(</>): V{4>) ^ 
U((j)). In the case that V{<j>) — U {<]>), Eq. (|II.311|) shows that p^ — 0, which corresponds to dust wi th w j, = p<j>/p<p = 0. 
Of course, the classical and quantum dynamics of such constrained theories are quite non-trivial jl!7j . 
For simplicity, we choose V(0) and U(<fi) to be constants: 

U{<t>) = U Q , V{4>) = V„ . (11.312) 

Then, if Uq = Vq, the EoS parameter w,/, vanishes. In the general case, one has w$ = jf-ry-- Let us now use T^ v 



as 



an energy- momentum tensor in the previous section. Eq. pi.301[) shows /3 = — 2(3^-1) = 2U^-w ' ^ ne can s i m P n fy 

T+^R^ + fiT^R = ePWQR^ + U R . (11.313) 

Here, Eqs. (|II.308|t . (|II.310|) . and (|II.312|I are used. Similarly, 7 in (|II.303j) has the following form: 7 = ^Xv ■ which 
gives, by using (|II.304|) . 

T^V^V^ + 7 T*V P V P $ = d^fid^V^Vv® + 2Z7 VVp$ . (11.314) 
Eq. (|II.313|) enables us to write down the z — 2 total action corresponding to (|II.298|) as 

S = J d 4 x^^ - a (d^d^R^ + U R) 2 - A Qd M <W + U j } . (11.315) 

On the other hand, the z = 3 total action corresponding to (jll.3051) is: 

S = I d4x ^g{^ ~ a (d^d^R^ + U R) (c^c^VpV, + 2U V P V P ) (d^d^R^ + U Q R) 



-A ( -d^4> + U ) \ , (11.316) 



and 



S 2 n+2 = J d i x^\^-a{{d^d^V i y v + 2U^ p V p ) n {d^d v <j ) R^ + U R)} 2 



-A -dutn + Uo) } . (11.317) 



for the z = 2n + 2 model (n = 0, 1, 2, • • • ), and 

S 2 n+ 3 = j d A x^{^-a{{d^(t>d v 4>Wy.Wu + 2U V p V p ) n (d^d^R^ + UoR)} 

x {(0"^0»>V M V„ + 2£/ VV p r +1 {d»</>d v <j>R^ + U R)} - A + ET ) } , (H.318) 

for t he g = 2n + 3 model (n = 0, 1, 2, • • • ) (compare with 79]) 3 . Note that the actions (|II.315j) . (|II.316j) . (|II.317|) , 
and (|II.318|) are totally diffeomorphism invariant and are only given in terms of the local fields. The actions f|II.315[) . 
pi.316|) , pi.317|i . and (|II.318|> do not depend on Vq. 



3 Note that the model l)II.317|l or l)II.318|l has better renor maliza tion properties and no scalar graviton when the action is modified a 
little bit and the scalar projector is inserted in second term [116(1 . 
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By the variation over </>, for example for the z = 2 (n = 0) case in (|II.317|) . one finds 

= Aad" {d v <j>R„ u (d p c/)d a <t)R pr7 + U R)} + (Xd^) . (11.319) 

For the z > 3 (n > 1 case in (1II.31T[) or n > case in (111.3181 ). one obtains a rather complicated equation. 
The dispersion relation of the graviton is now given by 

uj = ac k z , (11.320) 

in the high-energy region. Here, cq is a constant, u) is the angular frequency corresponding to the energy and k is the 
wave number corresponding to momentum. If a < 0, the dispersion relation becomes inconsistent and, therefore, a 
should be positive. 



3. FRW cosmology 



The gravitational terms different from general relativity in (|II.317|) and (|II.318[) are relevant in the high-energy 
region. Such terms might affect the inflationary era. In this section, we briefly study FRW cosmology within the 
theory under discussion. To obtain the FRW equations, the following form of the metric is assumed: 

. 2 



ds 2 = -e 2b ^dt 2 + a(t) 2 ( dx T > 



i=l,2,3 



and it is assumed that the scalar field <f> only depends on time. Then, Eq. (1II.308[) looks as: 



Thus, one gets 



if) =^ 



d^-pd^R^ + U R = 6U c- 2b H 2 , 4>d v '0V M V \ + 2C/ V P V P = -GU e- 2b Hd t . 



and the actions pi.317|) and pi.318|) have the following form 

-b 



S- 



2n+2 



d 4 xa 3 



> 2 ^H+12H 2 -6bH)-(6Uo) 2n+2 e b {(e- 2b Hd t ) n (H 2 e- 2b )y 



-A - 



$2n+3 



d 4 xa 3 



3 — f^t 
2 \dt 

b 



e 
2^2 



+ e b U 
(6H + 12H 2 - 6bH 



_ 2 2, i+ 3 . 32^ U2n+2 e b^ e -2b m y (j^- 31 ) } { (e^tf d t ) " + ' (iJ 2 e - 2b ) } 



dt 



e b U 



The first FRW equation looks as follows: 

= ^H 2 -{6U ( . 



n2m+2 _3 

aa 



{D n (H 2 )) - A{-l) n H 2 D n (a 3 D n (H 2 )) 



n-l 



(-1)" (D n - k (H 2 )) (D k (a 3 D n (H 2 ))) 



— 2XUq — Pmatter j 



for pi.324[) and as follows: 



= ^H 2 - (6CT 



2n + 3 „,„~3 I" 3 I nn 



aa 



[a 3 (D n (H 2 )) (D n+1 (H 2 )) - 2(-\) n H 2 D n (a 3 D n+1 (H 2 )) 



n-l 



-2{-l) n+1 H 2 D n+1 (a 3 D n (H 2 )) - 2^ (-l) fe (D n - k (H 2 )) (D k (a 3 D n+1 (H 2 ))) 



k=0 



-2^2 ("if (D n+1 - k (H 2 )) (D k (a 3 D n (if 2 ))) 



k=0 



2XU a P matter j 



(11.321) 

(11.322) 
(11.323) 



(11.324) 



(11.325) 



(11.326) 



(11.327) 



4G 



for (|II.325[) . We have also put 6 = after the variation over b, where the metric (III.321[) reduces to the standard FRW 
metric and the operations of D and D for a scalar <p are defined by 

D ^ H % ' D ^^r- (IL328) 

On the other hand, by the variation over a, we get 

= -L( 2 H + 3H 2 ) + 2 2n+2 3 2n+1 aU 2n+2 { -3 (D n (H 2 )) 2 + 4(-l) n a' 3 ^ (HD n (a 3 D n (H 2 ))) 

+ 2 i>~ 3 4 f ( t+ ( Dn ~ k ( h2 ))) ^ ( a3Dn ( h2 )))) ) > ( IL329 ) 

for (|II.324[) and 



dt \\dt 

k=l v v 



= ^(2H + 3H^+2 2n + 3 3 2n+2 aU 2n+2 {-3(D n (H 2 ))(D n+1 (H 2 )) 

+2 (-1)" a- 3 ^- (HD n (a 3 D n+1 (H 2 ))) + 2 a' 3 ^- (HD n+1 (a 3 D n (H 2 ))) 

+i a ~ 3 j t ((I ( Dn - k (^ 2 ))) ^° n+i (^ 2 ))) 

+X>~ 3 ^ {{i ( D ™ ( ff2 ))) ^ ( ff2 ))) ) \ (IT - :wl, -» 



fc=i 



for (|II.325|) . The simplest case is n — in pi.324[) when the FRW equations are 

3 



^H 2 = -108a£#ff 4 + 2A[/ +/Wter, (11.331) 



~ (2H + 3H 2 ) = 36aU 2 (3H 4 + 4ff 2 ij) + p matter . (11.332) 

In the early universe, in which the curvature was large, the contribution from the Einstein term, which corresponds 
to the right-hand side in (|II.331|) and (|II.332j) . and the contributions from the matter p ma tter and ^matter, can be 
neglected. Then, a solution of (|II.332j) is given by 

H = j t , (11.333) 

which expresses the (power-law) accelerated universe expansion corresponding to that with a perfect fluid with w = 
-1/2. Eq. pi.33ip gives 

X-^P- (H.334) 

For more complicated versions of the above theory, one can expect to find a more rich class of background cosmological 
solutions. The cosmological reconstruction may also be used for such models. 



H. Dark fluid with an inhomogeneous equation of state 

Usually, a perfect fluid with a constant EoS parameter w is discussed in the cosmological-related literature. One 
may, however, consider a general equation of state fluid: 

p = w(p,a,H,H, H,--A p+ (p,a,H,H,H,---^j . (11.335) 

It is clear that any modified gravity may be presented as such an effective gravitational fluid coupled with general 
relativity as was already mentioned in section A of the present chapter. Such an approach to modified gravity can 
make things simpler when background evolution is investigated. However, it may lead to a number of problems in the 
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study of cosmological perturbation theory or structure formation. Indeed, when developing such an approach, one 
often forgets that the initial modified gravity is a higher-derivative theory and that perturbation theory should also 
be a higher-derivative theory. Non-linearities are also effectively dropped. That is why we prefer to work with the 
explicit modified gravity formulation but not with general relativity coupled to gravitational fluid. 

The above fluid formulation is very general. Moreover, even if we restrict it to the case t hat th e pressure p depends 
only on the energy density p, the EoS fluid, which arbitrarily reproduced in FRW cosmology |ll8j | , may be constructed. 
Needless to say, the unification of inflation with the dark energy epoch may already be achieved for such a fluid. The 
FRW equations pi.3[> show that, for the perfect fluid with the generalized equation of state 



P 




(11.336) 



the solution of the FRW equations is given by H = f(t). Here, f(t) is an arbitrary function of the cosmological time 
t. Therefore, the arbitrary universe evolution given by H = f(t) can be realized via the perfect fluid whose equation 
of state is given by (|II.336|> . 

It is often convenient to use e-folding TV = In ^- instead of the cosmological time t because e-folding is related with 
the redshift z, which is directly given by observations (for more details, see the next chapter). One may consider the 
reconstruction using N instead of the cosmological time t. Then, the second equation (JOj) can be rewritten as 



1 



Ptotal 



(2HH' + 3H 2 ) 



(11.337) 



Here, ' expresses the derivative with respect to N: ' = d/dN. Thus, the fluid with the generalized equation of state 



P=-p--f 

K 




(11.338) 



gives the solution of the FRW equations by H — f(N). 

It was already mentioned in section A that a phantom fluid whose EoS parameter w is less than —1 leads to the 
future Big Rip singularity Several types of softer future singularities may occur for an effective quintessence fluid. 
The classification of four future finite-time singularities is given in section A following ref. |34j |. 

The Type I singularity corresponds to the Big Rip singularity [l|, |33j, which emerges when w < — 1. The Type II 
singularity corresponds to the sudden future singularity [3a ] at which a an d p are finite but p diverges. The Type III 
singularity appears, for instance, for the model with p = —p — Ap a |l!9| . which is different from th e su dden future 
singularity in the sense that p diverges. This type of singularity was discovered in the model of ref. |l!9l |. where the 
corresponding Lagrangian model of a scalar field with a potential was constructed. 

Let us show that all above four types of future singularities may naturally occur for fluid dark energy. As an 
example, one may start from the following fluid dark energy 



P=-P~ f(p) 



(11.339) 



where f(p) is ge nera lly an arbitrary function. The case of / (p) oc p a with a constant a was proposed and investigated 
in detail in ref. [l!9j |. Using the conservation law (jII.I44l) for such a choice, the scale factor a is given by 



Using the first equation (|II.3p . the cosmological time t is 



dp 



In the case 



by using Eq. (III.340j) . it follows that 



Ky/3pf(p) 



f(p) = Ap a , 



(11.340) 



(11.341) 



(11.342) 



a = a,Q exp 



3(l-a)A_ 



(11.343) 
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When a > 1, the scale factor remains finite even if p goes to infinity. When a < 1, a — >■ oo (a— s-0) as p — !► oo for 
A > (A < 0). Because the pressure is now given by 

p = -p - Ap a , (11.344) 

p always diverges when p becomes infinite. If a > 1, the EoS parameter io = p/p also goes to infinity, that is, w — > +oo 
(-oo) for A < (A > 0). When a < 1, we have w -> -1 + (-1 - 0) for A < (A > 0) as p -> oo. 
Using Eq. ljll.341li for (|II.342L one finds 

t = t. + for a*-, (11.345) 



and 



In 



Po / 



1 



t = t s + J: 7 , for a = - . (11.346) 
V3kA 2 

Therefore if a < 1/2, p diverges into an infinite future or past. On the other hand, if a > 1/2, the divergence of p 
corresponds to a finite future or past. In the case of a finite future, the singularity could be regarded as a Big Rip or 
a Type I singularity (compare with the case in section A). 



For the equation of state (|II. 342ft , the following cases were discussed [34{ : 



• In the case a = 1/2 or a — 0, there does not appear any singularity. 

• In the case a > 1, Eq. (|II.345[) indicates that when t — > t s , the energy density behaves as p — > oo and, therefore, 
\p\ —> oo due to (|II.344[) . Eq. (|II.343I) shows that the scale factor a is finite even if p — > oo. Therefore, the a > 1 
case corresponds to a Type III singularity. 

• The a = 1 case corresponds to p — wp with a constant w if we replace — 1 — A with w. Therefore, if A > 0, 
either the Big Rip or a Type I singularity occurs, but if A < 0, no future singularity appears. 

• In the case that 1/2 < a < 1, when t —> t s , all of p, \p\, and a diverge if A > 0; this corresponds to a Type I 
singularity. 

• In the case that < a < 1/2, when t — > t s , we find p, \p\ — > and a — > do, but by combining (III. 3431) and 
(III.345I) . we find 

In a - |f-t s |^=T72 . (11.347) 

Because the exponent (a — l)/(a — 1/2) is not always an integer, even if a is finite, the higher derivatives of H 
diverge in general. Therefore, this case corresponds to a Type IV singularity. 

• In the case that a < 0, when t —> t 8 , we find p — >• 0, a — >• ao, but |p| — > oo. Therefore, this case corresponds to 
a Type II singularity. 

Thus, it is explicitly demonstrated that simple fluid dark energy may bring the universe evolution to one of four 
possible future singularities as well as to an asymptotically non-singular universe. The particular future behavior is 
defined by the choice of the equation of state and of the parameters of this equation. Precisely the same phenomenon 
occurs for any field model of dark energy or modified gravity, so there is no qualitative difference between different 
dark energy models in this respect. The emergence of a finite-time future singularity in different modified gravities 
will be described in the fourth chapter. 

In summary, this chapter provided a review of several popular models of modified gravity is given. Of course, 
the explicit choice of the models is related to our scientific interests. In particular, the following models were dis- 
cussed: F(R) gravity, modified Gauss-Bonnet theory, scalar-Einstein-Gauss-Bonnet gravity, non-local gravity and 
non-minimally coupled gravity. More exotic theories with explicit Lorentz symmetry breaking like Hofava-Lifshitz 
F(R) gravity or with apparent Lorentz symmetry breaking like covariant power-counting renormalizable field gravity 
were also presented. In all cases, the spatially flat FRW cosmology was investigated. The qualitative possibility of 
a unified description of early-time inflation with late-time acceleration was also demonstrated. It turns out that the 
easiest and very realistic unification of these two eras is possible in F(R) gravity. 
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III. COSMOLOGICAL RECONSTRUCTION OF MODIFIED GRAVITY 



To explain the unified cosmic history of the universe, especially the inflation and/or the late-time accelerating 
expansion of the present universe, several models of extended or modified gravity were proposed in the previous 
chapter. Usually, we start from a theory, which is defined by the action, and solve equations of motion to define the 
background dynamics. In this chapter, however, we consider the inverse problem, i.e., the cosmological reconstruction 
of gravitational theories. In other words, using the fact that modified gravity is defined in terms of some arbitrary 
function (or with the help of some potential), we show how the complicated background cosmology, which complies 
with observational data, may be reconstructed. Eventually, this leads to a particular choice of the above arbitrary 
function (or potential). This provides the qualitative difference between general relativity and modified gravity; the 
latter theory (more exactly, some class of it) may have as a solution any given gravitational space-time. The problem 
is, of course, how to define the specific version of the theory that has this particular space-time as an explicit (exact 
or approximate) solution. The reconstruction scheme is explicitly constructed here for several modified gravities 
introduced in the previous c hapter . The general approach to reconstruction in modified gravity and dark energy 
models was developed in refs. pi lilol lilTj . Th is kind of reconstruction may also be used for a spherically symmetric 
solution like black holes (see, for example, [l22j ). 



A. Scalar-tensor gravity 

As the first example, let us consider the reconstruction of scalar-Einstein gravity (or scalar-tensor theory), whose 
action can be written as 

S = J d^x^-g I JL R - W - V(cf>) + L mattcr | . (III.l) 

Here, u>(<p) and V(<p) are functions of the scalar <fi. The function oj((p) is not relevant and can be absorbed into the 
redefinition of the scalar field (j>. In fact, if one redefines the scalar field cj) by 

<p = J^d^^JWi, (in.2) 

the kinetic term of the scalar field in the action (|III.1|) has the following form: 

,quj. / -dutpd^tp when uU>) > njj o\ 

-^)d^ = [ d ^ <p wh en ^) <0 • ( IIL3 ) 

The case of w(</>) > corresponds to the quintessence or non-phantom scalar field, but the case of &((/>) < corresponds 
to the phantom scalar. Although ui(<j>) can be absorbed into the redefinition of the scalar field, we keep ui((f>) for a 
later convenience. 

1. The reconstruction using the cosmological time 



This section is based on the studies 123Ul25 [. 



For the action (jlll.ip . in the FRW equations (|II.3|) . the energy density and the pressure are as follows: 

P = i w (0)0 2 + y(0), P = i w (0)0 2 - v{cj>) . (IH.4) 



Then 



~H , V((f>) = \ (m 2 + H) . (III.5) 



Assuming cu((j>) and V((j>) are given by a single function f(<f>), as follows, 



"(4>) = -^f'W , v{4>) - 4 (W) 2 + /'(£>) . ( IIL6 ) 
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the exact solution of the FRW equations pi.3[) with (|III.4I) (when we neglect the contribution from the matter) has 
the following form: 

cf> = t, H = f(t). (III.7) 
It can be confirmed that the equation given by the variation over 4> 

= oj{4>)4> + ^'(<p)<j> 2 + 3Huj(ct>)<i> + V'{4>) , (111.8) 

is also satisfied by the solution (|III.7I) . Then, the arbitrary universe evolution expressed by H = f(t) can be realized 
by an appropriate choice of w((/>) and V((f>). In other words, by defining the particular type of universe evolution, the 
corresponding scalar-Einstein gravity may be found. 

As mentioned in (jIII.2l) and (|III.3|) . uj(4>) can be absorbed into the redefinition of the scalar field </>. By keeping ui((j>), 
however, we can construct a model that exhibits a smooth transition from the non-phantom phase to the phantom 
phase (and vice versa). 

Let the matter energy density p ma ttor and the matter pressure Pmatter be given by the sum of the contributions from 
the matter whose EoS parameters are u^attcr = Pmatter/Pmatteri 8 a constants. Here, the suffix "i" specifies the type 
of matter. Using the conservation law of the energy density (jll.1441) . we find /? ma ttcr = X) / 9 o a_3 ' 1+tUmattor ' > - Here, p l 
is a constant. If uj(4>) and V((f>) are given in terms of a single function g((f>) as 

K, z 

= - 2 (3 ff '(0) 2 + g"(<f,)) + £ <a T ~ V o^o 3(1+ "" ) e- 3 ( 1 +^— >«W , (III.9) 

i 

the solution of the FRW equations (|II.3|) is 

tf> = t, H = g'(t) , (a = a e 9(t) ) . (111.10) 

Thus, even in the presence of matter, any given cosmology defined by H = g'(t) can be realized by potentials ((IIL9]). 

As an example, the model with only one kind of matter with the EoS parameter w ma ttcr (|III.9j) can be considered. 
It is also assumed that Wmattor > —1- Choosing g((f>) as 



one obtains 



3 (^matter + 1) \t s ~ 



, , _ _^ i 1 1 \ W ma ttcr + 1 -3(l+nwtor) ( f s ~ 4>Y 



V W - 4 " 1 ^ 2/11 



« 2 1 3 (uWtcr + 1) V <P U - <j> J 3 («Wter + 1) ^ V (t s - 

2 



Wmattcr " 1 n „-3(l+«Wtcr) ftj ~ 4>) 

The Hubble rate H follows as 



^matter ^ - J(,l+m„attcr \"s YJ ("TTT 1 

H Poa ~2 • (111.12) 



Because 



^ 3 (iy ma tter + 1) \* + * s -*/ 



g= 3K nattcr + i) )' (IIL14) 
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the EoS parameter w c g (|II.4[) goes to w mat tcr > — 1 when (40 and goes to —2 — w ma ttcr < — 1 at late times. The 
crossing w c g = — 1 occurs when H = 0, that is, t — t s /2. Note that 

a_ 16£ s (3w mattC r + l)t s \ (III 15) 



a 27 (uwter + 1) (t, - ty t 2 

Thus, if u; ma tter > — 1/3, the deceleration of the universe turns to the acceleration at t = t a = (3i0matter + 1) t s /A. 
The second example in (|III.9|) is given by 



g(cj>) = -alnl l-p'ln^ 1 . (111.16) 

Here, a and (3 are dimensionless positive constants. Let us choose (3 ~ O (lO -2 ) for later convenience. Note that 
the parameter order of 10~ 2 is not so unnatural. When only one kind of matter with the EoS parameter iu ma ttcr is 
included, for simplicity, Eq. (jlll. 16[) gives the following expressions for oj((f>) and V((p): 

2ap(l3-l + p\n£) ... , i , / ± \ 3Katte,+i)« 

w(W = ; T2 n Poa 1-pTn- 

ac 2 (l-pMn^J </> 2 

a/3(3a/3 + p , -l + pTn^ N ) ,„ i / , n 3(«, mattet +i)a 

VW = P V P n 2 P ^ + Wmatt o er ~ 1 p a^ 1+w — > ( 1 - fi In g ) . (111.17) 



Ac 2 (l-/31n£J 2 \ 

Supposing pqci 3 ( 1 +"'™** tBr ) ^ ( K ~ 2 )i an y small parameter, compared with the Planck scale does not appear. Now, 
the Hubble rate is given by H = a/3/ (1 — pTn (t/n)) t, which is positive if 

< t < t s = Ace? , (111.18) 
and has a Big Rip-type singularity at t = t s . Because 10 61 ~ e 140 , with the choice (3 ~ 1/140, we obtain t s ~ 

k x 10 61 ~ (lO~ 33 eV) , whose order is that of the age of the present universe. Thus, due to the property of the 
exponential function (or logarithmic function), the small scale like t s appears rather naturally. We should also note 

that if a/3 ~ 0(1O 0-2 ) and t is a present age of the universe t ~ (l0~ 33 eV) , the observed value of the Hubble rate 
H ~ 10~ 33 eV can also be reproduced. Because 

« (1- /?m£) 2 i 2 

the universe changes to an accelerating from a decelerating universe when 



(111.19) 



t = t a = ne^^ 1 < t s . (111.20) 
Because the energy density p of the scalar field <f> and that of the matter p ma tter are given by 

3Q P -3(l+«Wt.r) I i fli * 

p = o Poa n 1 — p In - 

K 2(l_/31 n i) 2 t 2 Wo V 

(, \ 3(-u) m attcr + l)a 
l-/?ln-J , (111.21) 

the coincidence time t c could be given by solving the following equation: 



AC 



1-/3 In- *c= I _,nl „> • ( IIL22 ) 



9 -3(l + tUmatter) 



One may regard p ma tter as the sum of the energy density of usual matter, like baryons, and that of (cold) dark matter. 
If p corresponds to the energy density of the dark energy, the current data indicate that p : p ma ttor ~ 7 : 3. Then, in 
the present universe, it follows that 

. \ 3(u> ma ttcr + l)a+2 Q 2 al 

1-0 In-) t 2 ~%Hr~- (IH.23) 



At 



Thus, in the model (|III.16|) . the acceleration of the present universe and the coincidence problem seem to be explained 
rather naturally. This is achieved via explicit cosmological reconstruction of the scalar potentials. 
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2. The reconstruction of the scalar-tensor models using e-folding 



In this subsection, the formulation using the e-folding TV instead of cosmological time is considered as in Eqs. (|II.337[) 
and (III.338[) . The FRW equations (|II.3[) with pil.4[> (when the matter contribution is neglected) are given in terms 
of N. as follows: 



2 K z ' 

We now identify <fi = TV; then, 



= V(cf>) = -^(HH> + 3H 2 ) . (111.25) 



The above equations (|III.25[) show that if one studies the model 

CJ W = -%7^' ^) = 4(/>)/'W+3/(0) 2 ) , (IIL26) 

K 2 / (0) V ' 

which is given in terms of a single function /, then the exact solution is given by H (N) = f (N). 
Let us consider an example by using the following function /: 

f{4>) = Hi^ . (IH.27) 

Here, Hi and 7 are positive constants. In this case 

w (0 = nnz > V(cb) = ^ (— 7^~ 27_1 + 3-T 27 ) , (111-28) 

and the exact solution is given by if = HiN^ 1 . With a redefinition of by — > <p = 2 v^7^ ; the kinetic term becomes 
canonical, and the potential is given by 

At late times (large N), only the last term is relevant. Conversely, at late times, if the potential is given by V oc <y9~ 47 , 
the solution has the form of H cx TV -7 . The effective EoS parameter (III. 41) is equal to 

w eff = -l+|L, (in.30) 

which goes to — 1 when N goes to infinity. 

Because we naively find C = Plan | o scale - 10 61 ~ e 140 , if ~ C, it follows that N ~ In C and, therefore, 7 ~ 28. 
Thus, the fine-tuning problem might be relaxed, b ut in this case one obtains w c g ~ —0.87, which is a little bigger 
than the observational value -0.14 < 1 + w < 0.12 [l26j . 

In the case that C = Weak scale _ 1Q 44 _ e ioi ; if h± _ ^ it fo n ows that JV ~ InC and 7 - 21. This finishes the 
cosmological reconstruction of scalar-tensor theory using e-folding. 



3. Two-scalar model 



Note that when studying the transition from the non-phantom p hase to th e phantom phase, the instability of the 
single scalar-tensor gravity becomes infinite at the transition point jl23l . Il24| . Thus, we only investigate two scalar- 
tensor theory, which has an instability that can always finite. Moreover, with two scalars it is easier to realize unified 
universe history expansion. 

The action of two scalar-tensor gravity is 

5 = | ^xV^l^jje-iw^^-i^xj^x^x-vw.x)} ■ (111.31) 
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Here, u>(<j>) is a function of the scalar field (j>, and T](x) 1S a function of another scalar field \. Note that the above 
theory parameterizes a large number of popular dark energy models. Indeed, if both scalars could be canonical scalars, 
then a kind of the quintessence model emerges. One of scalars may be phantom (wrong sign of the kinetic term) while 
another scalar remains the canonical one. This is dark energy sometimes called quintom theory. Finally, both scalars 
could be phantom. Moreover, two scalars may be organized so that a charged scalar with real and imaginary parts 
appears. 

It is also assumed that the FRW metric is spatiallyflat (|II.6|) and that <f> and \ only depend on the time coordinate 
t. The FRW equations lead to 

W (0)0 2 + r,( X )x 2 =-\h, V{4>, x) = \ ( 3R2 + h) ■ ( IIL32 ) 



if 



«(*) + v(t) = -4/'(t) , V(t, t) = 4 (3/(t) 2 + /'(*)) , (HI-33) 



the explicit solution follows as: 

d> = X = t, H = f(t). (111.34) 
One may choose uj to always be positive and r\ to always be negative, for example 

<"(0 = —T {/'(*) " V^TTW} >0, 'j(x)=4v / « 2 + /'W 2 <0. (111.35) 
Here, a is a constant. Define a new function /((/>, x) by 



/(0,x) s — (j d<M</>) + y ^(x)J • (HI.36) 

The constant of the integration could be fixed to require 

f(t,t)=f(t). (IH.37) 

If V((/),x) is given using f(<j>,x) as 

x) = ^ (3/(0, x) 2 + + , (IH.38) 

the FRW and the scalar field equations are also satisfied: 

= uj(4>)4> + \u'm 2 + ZHwWi + dV{ ^ x) , = v (x)x + \v'(x)x 2 + 3H V (x)x + dv{ ^ x) . (111.39) 

Thus, for given cosmology the cosmological reconstruction (the calculation of scalar potentials) is made. 

It is now interesting to investigate the (in)stability of the model under discussion. By introducing the new quantities, 
X$, X x , and Y as 

= X X = X , Y=^^, (HI.40) 
the FRW equations and the scalar field equations (jIII.39[) are: 

dx^ . 2 dx x _ i{ x ) , 2 . 

""2i^M (^-l)-3(^-n, ~M--2Hrid {X x -l)-S(X x -Y), 
dY 1 

diV = 2^ - 1) + *x - 1)} • (HI.41) 

Here, d/dN = H~ 1 d/dt. In the solution (|III.34j) . = X x = Y = 1. The following perturbation may be considered 

X = 1 + SXcj, , X x = 1 + 5X X , Y = 1 + SY . (111.42) 
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Thus, 



d 

dN 



sx x 

SY 



= M I SX X 

SY 



M = 



\ 





1 



v'ix) 



- 3 



(111.43) 



The eigenvalues of the matrix M are given by solving the following eigenvalue equation 



A + 



2k 2 H 2 



A + 



Ox) 



v'(x) 

3 



A - 



1 



K 2 H 2 



2k 2 H 2 



A + 



Hu)((j>) 



(111.44) 



The eigenvalues (|III.44|I for the two-scalar model are clearly finite. Thus, the instability could be finite. In fact, right 
at the transition point where H = f'(t) = and, therefore, f'(<f>) = f'(x) = 0; f° r the choice in (|III.35|i . one gets 



w(0) = ~V(X) = ^ , w'(0) = -^?, V'(X)=0. 

K K 

Then, the eigenvalue equation (|III.44|) reduces to 

3 H 
= A 3 + {-A - B + 6) A 2 + (AB -3A-3B + 9)\- -AB + 9B, A = — , B 

2 a 



K 2 H 2 



(111.45) 



(111.46) 



Here, a > is chosen. The eigenvalues are surely finite, which shows that even if the solution (|III.34j) is not stable, 
the solution has a non- vanishing measure and, therefore, the transition from the non-phantom phase to the phantom 
phase can occur. We should also note that the solution ([III. 341) can be stable. For example, for the case A, B — > 0. 
Eq. (|III.46j) further reduces to 



= A(A + 3) 2 



(111.47) 



Then, the eigenvalues are given by and —3. Because there is no positive eigenvalue, the solution (|III.34[) is stable 
in this case. 

It is not difficult to extend the above formulation to the multi-scalar model, whose action is given by 



S = Id xy/—g 



Here, uii(4>i) is a function of the scalar field (pi. We now choose to satisfy 

2 



i 

by a proper function. The potential V(cj)i) is chosen as 



W = 4r(W) + £|£ 



Here 



K 2 r 
/(&) = — 2" 2-/ / WiUiifc) 

The constant of the integration in (|III.51|) is determined to satisfy 



K<k) 



fit). 



all <pi=t 

Then, a solution of the FRW equations and the scalar field equations is given by 

4>i = t, H = f(t). 



(111.48) 



(111.49) 



(111.50) 



(111.51) 



(111.52) 



(111.53) 
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The remaining arguments coincide with those given for the two-scalar model. 

It is interesting to demonstrate that the transition from the matter-dominant period to the acceleration period can 
be realized in the present formulation. In the following, the matter contribution can be neglected because the ratio 
of the matter with the (effective) dark matter may be small. 

The first example is as follows: 

H = fit) = g'(t) =9o + j- (111.54) 

When t is large, the first term in (|III.54|) dominates and the Hubble rate H becomes a constant. Therefore, the 
universe is asymptotically the de Sitter space, which corresponds to an accelerating universe. On the other hand, 
when t is small, the second term in 54[) dominates and the scale factor behaves as a ~ t Bl . Therefore, if g\ — 2/3, 
the matter-dominated period emerges. Here, one of the scalars may be considered as usual matter. 
Because 

f'(t) = -^<0, (111.55) 
instead of (|III.35j) . by using a positive constant a, one may choose 



o#) = - 2(1 | Q) /'W>0, V(x) = ^f(4)<0, (IIL56) 



that is, 



2(1 + a) g x 2a g x 

W (W = — — . ViX) = -j ■ (111.57) 



One should note that in the limit a — > 0, the scalar field x decouples and the single-scalar model follows. Thus, one 
gets 

it a \ , + a 9i \ 1 )o( , 0- + a )9i &9l\ 2 + , «9x \ /TTT 

fW,X)=9o + , V(<p, x) = -o i 3 ffo + — 2 + -Tf' (111.58) 

9 X « [ V 9 X J 9 X J 

Thus, for the scalar-tensor theory with such potentials the matter-dominant era occurs before the acceleration epoch. 
This fact is easily established within the reconstruction scheme. 

Before going to the second example, we consider the Einstein gravity with cosmological constant and with matter 
characterized by the EOS parameter w. The FRW equation is as follows: 

^H 2 =p a- 3 ^ + 4 1 2- ( IIL59 ) 
Here, I 2 is the inverse cosmological constant. The solution of (|III.59|) is given by 

a = a e 9(t) , g(t) = 2 - - In (asivh ^ ± w) (t-f.))V (HL60) 



Here, t s is a constant of the integration and 



2 _ 1 2/2 —3(1+11)) / TTT e1 \ 

a = -k I p a . (111.61) 



Thus, this is the ACDM era as it follows in the commonly accepted approach. 

As the second example, the reconstruction of the ACDM metric (IIII.60[) for scalar-tensor theory is done below. 
Because 

/(*) - ✓(*) - y -th (^0 {t _ tj) ) > m = m = sinh - 2 (30 + «o {t _ t ^ < (m 62) 
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it is convenient to use pil.561) instead of (IIII.35I) . Then, one arrives at 
3(l + tfl)(l + a) . , _ 2 {3(1 + w) 



w(0) 

vix) 
f(<t>,x) 



K 2 l 2 



■ sinh 



21 



(<p-t s ))>0 



3(l + w)a . , _ 2 /3(l + io) . . , 

i . sm h ( v ., ' {x-Q) < , 



n 2 l 2 



l + a , /3(1 + 
— - — coth 1 



21 



\ 21 



p { (1 + a) coth 



M)-fcoth(^±^(x-* s )). 



3(1 + ib) 



2/ 



— i a ) I — a coth 



3(1 + w) 
21 



(x ~ t.) 



(111.63) 



3(l + w)(l + q) 
I 2 



sinh 



3(1 



2/ 



(<f> -t s )) + 



3(1 



I 2 



^sinh" 2 



3(1 



21 



This completes the reconstruction of scalar potentials to reproduce the ACDM epoch. 

Thus, in both examples, (|III.54| and (|III.60j) . the matter dominant stage, the transition from the matter dominant 
phase to the acceleration phase and acceleration epoch occur. In the acceleration phase of the above examples, 
the universe asymptotically approaches the de Sitter space. This does not conflict with the Wilki nson Mic rowave 
Anisotropy Probe (WMAP) data. Indeed, three years WMAP data have been analyzed in ref. jl27l - ll29l |. The 
combined analysis of WMAP data with the Supernova Legacy Survey (SNLS) constrains the dark energy equation 
of state wde, pushing it toward the cosmological constant. The marginalized best fit values of the equation of state 
parameter at the 68% confidence level are given by —1.14 < wde < —0.93. In the prior case of a flat universe, the 
combined data give —1.06 < u>de < —0.90. In the examples (|III.54[) and (|III.60|) . the universe goes asymptotically to 
the de Sitter space, which gives u>de — > — 1, which does not conflict with the above constraints. Note, however, that 
one needs to fine-tune go in pil.54[) and l/l in (|III.60[) to be <7o ~ l/l ~ 10 -33 eV, to reproduce the observed Hubble 
rate H - 70 kms^Mpc- 1 ~ 10" 33 eV. 

As a final remark, there is no problem regarding the inclusion of the usual matter (for example, ideal fluid and 
dust) into the above reconstruction scenario and/or finding the scalar potentials corresponding to realistic cosmology 
in the multi-scalar case. Moreover, in the same way, one may include the radiation-dominated epoch where quantum 
effects may still be neglected in the above scenario. In the next section it is shown how similar reconstruction scheme 
may be developed for modified gravity. 



4- Brans-Dicke gravity 
The action of the original Brans-Dicke model [l3C)j | has the following form: 

which can be generalized as 



1 /V M^±\ 



S 



M4>)r i 



2k 2 



(111.64) 



(111.65) 



This kind of model was applied to the dark energy problem in [13lj . and it was found that even if u>((j>) > 0, that is, 
the scalar field is canonical, the phantom universe can be realized. 
The FRW equations are as follows: 



H 2 = 



This can be rewritten as 



J2 



2e v 

T2~ 



H+-(cp + <p 2 



H = — 



2e v 



H + 3H 2 + - (Cp + ^ 2 ) 



Then, if the following choice of scalar potentials using a function /((/>) [132] is taken: 



/'(0) + -(^"(0)+^(^ 2 ) 



2 e v{<t>) 



/'(0)+3/(0) 2 + -( ( p"(0) + ^(0) 2 ) 



(III.! 



(111.67) 



, (111.68) 
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the explicit solution looks like 

H = f(t), = t. (111.69) 

Note that ip(<j)) can be an arbitrary function of <p. By this reconstruction formalism, the phantom crossing (w = — 1 
crossing) can be realized for uj(4>) > 0, that is, when the scalar fie ld is canonical (non-ghost). A slightly different 
version of the reconstruction formalism for above theory is given in |l33l ] . 



B. The fc-essence model 



Here, based on we consider the reconstruction of the fc-essence model, which represents the evident general- 

ization of quintessence theory. The fc-essence model is a rather general model that includes only one scalar field and 
the action is given by 

S = J <£xj—g ^ A _ K (4>, X) + L raattor ) , X = d^d^ . (111.70) 

Here, <f> is a scalar field. 

Now, the Einstein equation has the following form: 



^2 ( R *»> - \d^ R ) = ~ K (<f>> X ) 9»« + 2K X (cf>, X) d^dA + . ( TTL7 I ) 



Here, Kx (</>, X) = dK (<f>, X) jdX. On the other hand, the variation over <fi gives 

= -K4, (</>, X) + 2V" (K (0, X) d^) . (111.72) 
Here, K§ (</>, X) = dK (0, X) /d(f> and it is assumed that the scalar field 4> does not directly couple with the matter. 



1. The reconstruction using the cosmological time 

Let us start from the FRW universe (|II.6[) and assume that the scalar field tj> only depends on time. The FRW 
equations are given by 

-|ff 2 = 2X dK( ^ x X) -K(tf>,X) + p BaMm (t), -^(2H + 3H 2 ) = K X) + Pm&ttex (t) . (01.73) 
The FRW equations pil.73[) show that for the following model 

K(<j>,X) = jr / (X + l) n KWffl, KW{<t>) = -\{2f{<t>) + m<l>?) . KW(<f>) = \f'(<f,), (in-74) 

K K 

there exists a solution given by 

H = f(t), cf> = t. (111.75) 

Note that in (|IIL74|) . K^ n \(f>) n > 2 can be an arbitrary function. To investigate the role of K( n \(f) (n > 2), the 
following perturbation may be discussed: 

H = f(t) + SH, 4> = t + 6(f», (111.76) 

The stability of the solution (IIII.75|) should now be investigated. In (|III.76[) . it is assumed that 6H and S(j> do not 
depend on the spatial coordinates. Then, 5H is solved as a function of S(f>: SH = SH(S(f>) 7 and one obtains 



d5<j) 

In 



q 9 a \ K / e ^(2) 1 . 



(IIL77) 



If the quantity inside [ ] is negative, the solution pil.751) becomes stable. Because Eq. (jlll. 77[) contains K^ 2 \ if K^> 
is chosen properly, the solution corresponding to an arbitrary evolution law becomes stable. 
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Let us consider the condition that there is a Schwarzschild or Schwarzschild-(A)dS solution. Let <f> be a constant: 

4> = 0o ■ (111.78) 

Eqs. (|III.71|) and (|III.72|) reduce to 

1 1 



yR^~2 g ^ R ) ^- K ^o,0)g^ + T^, (111.79) 
O = ^(0o,O). (111.80) 

When = 0, if -ftT(</>o,0) does not vanish, a solution of (|III.79|) is given by Schwarzschild- (A)dS space-time. On 
the other hand, if K (4>, 0) vanishes, the Schwarzschild solution, which is asymptotically Minkowski space-time, is 
a solution. The equation (jIII.80[) requires that, in general, (</>, 0) has an extremum or K (0, 0) is a constant 
independent of <f>. Especially if K (0, 0) vanishes identically, Eq. (|III.79[) gives the Schwarzschild solution. 
K ((f), X) can be written as 

K((^,X) = ^2(l+X) n K n (4). (111.81) 

Then, if K (<p, 0) = 0, it follows that 

^#„(0)=O. (IH.82) 

In particular, one may choose 

K 3 (4>) = -K {cj>) - K x {ct>) ~ K 2 (<f>) . (111.83) 

Thus, if the condition (|III.82j) or (|III.83|) is satisfied, the Schwarzschild space-time is always a solution independent 
of the value of <f> as long as is a constant. Here, any point source of matter creates the Schwarzschild space-time 
which generates the Newton potential. The correction to Newton's law does not appear. Note that the value of the 
scalar field tj> changes by the evolution of the universe, but as long as the condition (|III.82[) or (|III.83[) is satisfied, in 
a local region where <f> is almost constant, the correction to Newton's law is negligible. 

2. The reconstruction using e-foldmg 

Because it is often convenient to use the redshift z or e-folding N, we rewrite the FRW equations (|III.73[) as 

3_ H 2 = 2X dK ^' X) - K (<f>,X) + Pmattcr (N) , -- L (2HH' + 3H 2 ) = K {&X) + Pmat ter(N) . (111.84) 

Here, H' = dH/dN. If the matter energy density is given by a sum of the contributions from the matter with constant 
EoS parameters u>i, one finds 

Pmatter(A') = > J Paid k ' = P0i e ) 

i i 

2Wter(A0 = E W*" 3 ^) = E ^iPoi^ {1+Wi){N - No) ■ (HI.85) 

i i 

Here, pen's are constants. For general energy density Pmatter(^), because the conservation law 144[) can be rewritten 
in terms of N as 

/Matter W + 3 (/WterW + Pmattcr(iV)) - , (111.86) 

it follows that 

fWter(iV) = - j0mat tcr(A r ) - fatter W ■ (HL87) 
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With the help of the FRW equations (|III.84li . we find 

K(4>,X) = ~ (iH^ + 3H 2 ^j + PmattC r(N) + ^p' m&ttex (N) , 

2 dK(d>,X) _ 1 dH 1 
H dx - ^H— - - Pinattei (N) . (111.88) 

If a new variable G(N) = H(N) 2 is defined, the equations (|III.88[) can be rewritten as 

K{<j>,X) = -±(G'(A0+3G(A0)+p ma t ter (A0 + ^matter (N) , 

H2 dK^X) = l_ G , {N) _ l pLatteAN) (m 8Q) 

Using the appropriate function g<p(4>), if we choose 

K{4>,X) = Y\ * + 1) K (n) (0), 

K^{<j>) = -1 (g'^cp) + 3g^)) = , (111.90) 

the explicit solution for the FRW equations (|III.73P is found as 

G{N) =H(N) 2 = ff0 (7V) + y p matter (iV) , = 7V (X = -ff 2 ). (111.91) 

Here, ir( n )(c/>) with ii > 2 can be arbitrary. The reconstruction is fulfilled. 

Let us now investigate the (in)stability of the solution (|III.91|) . For this purpose, the perturbation from the solution 
(IIII.91|) is taken as follows: 



G(N) = Go (AT) + SG(N) (g (N) = g^{N) + y Matter C^)) , = N + 



5<t>. (111.92) 



Note that ./V-dependence in the energy density p ma tter is usually given by a fixed function as in (|III.85[) . Therefore, 
£/9mattcr = 0- Thus, the equations (|III.89[) give 

(9'm + 3^{N)) 6J>(N) - ^ (g^ + M^(JV) - |^(A0) - + 36G(N)) , 

= ^6G'(N). (111.93) 



Then, we find 



5<j)'(N) \ _ 1 ( A B\( 6(f>(N) 
SG'(N) ) ~ L{N) \ C D ){ SG{N) J ' 

C = (g'l(N) + ag^N)) L(N) - g'^Nf (a + gjj^ , D = -3L(N) + g'^N) (a + . 



Here, 



L(N) = g't(N) + 8k 2 K^{N) . (111.95) 
The stability of the solution requires A + D < and AD — BC > 0, which gives 

1 ^ (9'4N) - G' (N)) < 3 - • ( IIL96 ) 



L(N)Go(N) wv ' uv " 2G (N) 

(-2 5 ;(iV)^(AT) - 9g'^N) 2 - 6g';(N)G Q (N) + 3G' (N)4(N) + 6G' (N)g'^N)) > g'^N) . (111.97) 
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It might be possible to make the system stable by choosing L(N) and, therefore, K^ 2 \N) properly. Neglecting the 
matter, one finds g,p{N) = Go(N) and obtains 

3 > Hj^y , (G' (N)G'{(N) - 3G' (N) 2 - 6G'{(N)G (N)) > C%{N) . (IIL98) 

Thus, as long as the first condition is satisfied, the second condition can be always satisfied by choosing L(N), and 
therefore K^'(N), as an appropriate function. Using the above-developed general reconstruction scheme, any specific 
accelerating universe may be formulated as a solution of such a theory. 

C. Modified F(R) gravity 
1. The reconstruction using the cosmological time 

It is easy n ow to generalize the cosmological reconstruction scheme for modified gravity with F(R) action (|II.2p 
[l8L fl20l fl2ll ]. The action (jlOj) is equivalently rewritten as 

S = J d^x^~g j-L {P{4,)R + Q(0)) + £ mattcr J . (111.99) 

Here, P and Q are proper functions of the auxiliary scalar 4>. The presence of an auxiliary scalar may simplify the 
formulation below, making it more similar to the one for scalar-tensor gravity. In fact, by the variation over (f>, it 
follows that = P'(4>)R + Q'((f>), which may be solved with respect to <f> as <j> = 4>(R). By substituting the obtained 
expression of 4>(R) into (|III.99|) . one arrives again at the F(i?)-gravity: 

S = J d^x^—g | + £ raattor j , F(R) = P{<t>(R))R + Q(4>{R)) ■ (III.100) 

By the variation of the action ([III. 991) with respect to the metric g^, the standard spatially flat FRW equations are 
obtained as: 

= -6H 2 P(ef>) - Q{<t>) - 6H dP( -^ + 2 K 2 p niattel , (III.101) 
= UH + 6H 2 ) P{4>) + Q(0) + 2 d2P{ t {t)) + AH^PI + 2 K 2 Pmattor . (III.102) 



dt dt 
Simple algebra leads to the following equation 

n d 2 P(4>(t)) dPj^fij) ■ 2 mnrvtt 

= 2 — 3 2H h 4HP{4>) + 2k (^matter + ^matter) ■ (III. 103) 

dt J dt 

As one can redefine the scalar field freely, a very natural choice is <fr = t. It is assumed that /O ma ttcr and p ma ttor are 
the sum from the contributions of the matters with constant equation of state parameters Wi- In particular, when a 
combination of the radiation and dust is assumed, one gets the standard expression 

Pmattcr = Proa" 4 + PrfOa" 3 , Pmattcr = ^a~ 4 , (III. 104) 

with constant p r Q and pdo- If the scale factor a is given by a proper function g(t) as a = aoe 9 ^ with a constant ao, 
Eq. 102[) reduces to the second-rank differential equation: 

= - + ^"i^m + 2« 2 5> + P»o«o 3(1+U " ) e^ 3(1+tUl)9(0) ■ (HI.105) 

In principle, by solving (|III.105|) the form of P(4>) may be found. Using (|III. 1011) (or equivalently (|III.102p ). the form 
of Q{4>) follows: 

Q(4>) = -6 {g'm 2 Ptf) - Gg'W^f- + 2k 2 £ ^C^V^+^C*) . (III.106) 
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Thus, in principle, any given cosmology expressed as a = aoe 3 ^ can be realized as the solution of some specific 
(reconstructed) F{R) gravity. Some explicit examples of such a reconstruction are presented below. Of course, the 
reconstruction is applied to cosmological eras that are fundamentally important in modern cosmology. 

Model reproducing ACDM-type cosmology 

Let us investigate whether ACDM-type cosmology can be reconstructed exactly by F(R) gravity in the present 
formulation when we include dust, which can be a sum of the baryon and dark matter, and radiation. 

In Einstein gravity, when there is matter with the EOS parameter w and cosmological constant, the FRW equation 
has the form (|III.59I) . The solution of (|III.59[) is given by (|III.60|> . Let us show how it is possible to reconstruct F(R) 
gravity by reproducing such an epoch (|III.60|) . When matter is included, Eq. (|III.105[> has the following form: 

n d 2 P(<jj) 2 , /3(l + to),, A dPU) 6(1 + to) 0/3(1 + 10),, A „, /TTT1 „„ N 
= - T coth ~2l (< ^ s) ) U Smh 2l ^ U) I m (IILl07) 



^ ^ I ^ I X ' ( M + \ 1 1 „ „ "J I ^ , , -. « U I ^ ' / M + \ 1 1 



+-^PrO% lasinhl (</>-*«) J J + Pdoa (^asinh^ (4>-t s ) 

Because this equation is a linear inhomogeneous equation, its general solution is given by the sum of the special 
solution and of the general solution that corresponds to the homogeneous equation. For the case without matter, by 
changing the variable from <fi to z as follows, 



z = - sinh- 2 ' (t - QJ , (III.108) 

Eq. 10T[) without matter can be rewritten in the form of the Gauss hypergeometric differential equation: 

d 2 P r_ ~ . \ 1 dP 



3(1 + to) 



= z(l-z)— + j-(a + j3 + l)z 



, a(3P, 
dz 



7 = 4+ — 1 & + p + l = 6+ 1 ap = - —^- (III.109) 
3(1 +10) 3(1 +w) 3(1 + 10) 

whose solution is given by Gauss' hypergeometric function: 

/ 7, % T(7) ^ T(a + n)T(B + n) z" 
P = P F(a,P,r,z) = P — ^-+^V V ^7+ v a III.110 

Here, T is the T function. There is one more linearly independent solution like (1 — z)^~ a ~P F(pf — a, 7 — $, 7; z), but 
we drop it here, for simplicity. Using (|III.106|) . one finds the form of Q(<fi): 

Q = - 6 -^^F(a, ft 7; z) ~ 3(1 ; ( y + (1 12 ^ )P V (a + Ij + 1, 7 + 1; z) • (ULlll) 

From (|III.108|) . it follows that z — > when t = <fi —> +00. Then in the limit of t = (j> — > +00, one arrives at 
P(4>)R + Q(4>) — > PqR~6Pq/1 2 . Identifying P = 1/2k 2 and A = 6/1 2 , the Einstein theory with cosmological constant 
A can be reproduced. The action is not singular even in the limit of t — > 00. Therefore, even without the cosmological 
constant nor cold dark matter, the cosmology of the ACDM model can be reproduced by F(i?)-gravity. 

We now investigate the special solution of (|III.107j) . By changing the variable as in 108() . the inhomogeneous 
differential equation looks as follows: 



= z(l - Z )±L + [7 - U + p + l) z] d 4- &PP + r, ( _ z) -2(i+3»)/3(i + ») + £ { _ z y¥& 
dz z L V / J dz 

'^« VWM ' ^5(IT^»'•^ 3 °" ,<,+ •" , ■ < IIL112 » 

It is not trivial to find the solution of (IIII. 1 12|) . Let us consider the case that w — and z — > —00, that is, t —> t s . In 
the limit of t -> t s , Eq. (|III.112[) reduces to 

d 2 P dP 

° = + ^lz~ -vPP + V (-*. ) > (III.113) 



G2 



whose special solution is given by 



P = P {-z)- 2 '\ P = 



9rj 
25 



(III.114) 



In principle, other special solutions of Eq. 1121) could be found. This proves that, in the presence (or in the 
absence) of traditional matter, the standard ACDM cosmology can be reproduced exactly by F(R) gravity. 



2. Reconstructed model realizing the phantom divide crossing 

Using the ab ove-d eveloped method, we reconstruct an explicit model in which a crossing of the phantom divide 
can be realized [l35j . In other words, the modified F(R) theory that describes the transition from acceleration to 
superacceleration (phantom era), will be presented. 

A solution of Eq. 105[) without matter can be given by 



P(<f>), 



P(0) 

= -10 In 

P(<j>) = P+^ + + P 



to 



7+1' 



1± ^1 + 1007(7+ 1) 



where 7 and C are positive constants, to is the present time, and p± are arbitrary constants. 
From Eq. (|III.116[) . it follows that g{4>) diverges at finite </> when 

4> = t s = t C-V^ , 



(III.115) 
(III.116) 
(III.117) 
(III.118) 

(III.119) 



which shows that the Big Rip singularity could occur at t = t s [l|. One only needs to consider the period < t < t s 
because g(<f>) should be a real number. Eq. 1 16[) also gives the following Hubble rate H{t): 



H(t) = 



dg(4>) 



7fe)" 7 " + (7 + l)c(^)' 



7+1 



(III.120) 



where it is taken that <fi = t. 

In the FRW background (|II.6[) . as in (|II.3[) . the effective energy density and pressure of the universe are given by 
p cif — 3H 2 /k 2 and p e ff = — (2H + 3H 2 ^j /k 2 , respectively. Then, the effective EoS parameter w c s is defined as in 
pi.4[) . For the case of H(t) of Eq. (|III.120[) , from Eq. (|II.4|) . we find that w e g is expressed as 



w e s = -1 + U(t) 



where 



-7 + 4 7 (7 + 1) (£) 27+1 + (7 + 1) (^) 



U(t) 



2H 
3IP 



2(27+1) 



15 



27+1 



(III.121) 



(III.122) 



7+(7 + l)(^) 

For the case of Eq. (|III.120[) . the scalar curvature R = 6 LH + 2H 2 ^j is given by 

7 (20 7 - 1) + 44 7 (7 + 1) (^) 27+1 + (7 + 1) (20 7 + 21) (^) 
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R 



2(27+1) 



t 2 



27+1" 



(III.123) 
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In deriving Eqs. (IIII.122I) and (IIII.123I) . Eq. (|III.119I) was used. 

When t ~ > 0, i.e., f<t s , H(t) behaves as 



H(t) 



IO7 



In this limit, it follows from Eq. (|II.4[) that the effective EoS parameter is given by 

-I 



1 

157 



(III.124) 



(III.125) 



This behavior is identical to that in Einstein gravity with matter with an EoS parameter greater than —1. 
On the other hand, when t —¥ t s , one finds that 



H(t) 



10 



-10 



t 



In this case, the scale factor is given by a(t) ~ ao (t s — t) lu . When t — > t s 
singularity appears. In this limit, the effective EoS parameter is given by 



w c ft 



-1 



1 

15 



16 
15 



(III.126) 

therefore, a — > 00, namely, the Big Rip 

(III.127) 



This behavior is identical with that of the case in which there is phantom matter that has an EoS parameter smaller 
than — 1. Thus, we have obtained an explicit model showing a crossing of the phantom divide. 

It follows from Eq. (|II.4I) that the effective EoS parameter w e ft is equal to —1 when H = 0. Solving w c g = — 1 with 
respect to t by using Eq. (|III.121[) . namely, U(t) = 0, we find that the effective EoS parameter crosses the phantom 
divide at t = t c given by 



tc t* 



2 7 + J47 2 + 



7+1 



1/(27+1) 



(III.128) 



From Eq. (|III,122[) . one sees that when t < t c , U(t) > because 7 > 0. Moreover, the time derivative of U(t) is given 
by 



dU(t) 
dt 



2 7 (7 + l)(2 7 +ir 



15 



7+(7 + l)(^) 



27+1 



1 3 



2 7 



1 - 



27+I' 



(III.129) 



Eq. (|III.129[) shows that the relation dU(t)/ (dt) < is always satisfied because only the period < t < t s is considered 
as mentioned above. This means that U(t) decreases monotonically. Thus, the value of U (t) evolves from positive to 
negative. From Eq. (IIII.121I) . we see that the value of w c g crosses —1. Once the universe enters the phantom phase, it 
stays in this phase, namely, the value of w e ft remains less than —1, and finally the Big Rip singularity occurs because 
U(t) decreases monotonically. Note that other types of finite-time future singularities in modified gravity are possible 
as demonstrated after Eq. ([II. 961) in the second chapter (see also ref. [3{|). It follows from Eqs. 1 1 5[) . (jIII.116[) . 
PII.117|) and (|III.119|) that Pit) is given by 



P(t) 



27+1 



3=± 



(III.130) 



Using Eq. (IIII.130I) . one gets 



Q(t) 



-6H 



ay 



\ 27+1 



E 

3=± 



3 -H 
2 



k 

t 



it** 



(III.131) 



If one can solve Eq. (|III.123|> with respect to t as t — t(R), in principle we can obtain the form of F(R) by using 
this solution and Eqs. (IIII.100I) . (|III.130[) and (|III.131|) . In fact, however, for the general case it is difficult to solve 
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FIG. 2: Behavior of t 2 s F(R) as a function of R for 7 = 1/2, p+ = -l/ts + , p_ = 0, /3+ = (l + 2vT9) /2 and t a = 2t . 



Eq. (|III.123p as t = t(R). Thus, as a solvable example, we show the behavior of t 2 F(R) as a function of R = t 2 s R in 
Fig. 1 for 7 = 1/2, p+ = -l/if + , _p_ = 0, /3+ = (l + 2\/l9) /2 and t s = 2t . The quantities in Fig. |HI_C2] are shown 
as dimensionless value. The horizontal and vertical axes show R and t 2 F, respectively. (Here, R = t 2 s R = AR/Rq, 
where Rq is the current curvature. In deriving this relation, we have used t s = 2to, to ~ H^ 1 , where Hq is the present 
Hubble parameter.) From Fig. 1, we see that the value of F(R) increases as that of R becomes larger. 

To examine the analytic form of F(R) for the general case, we will now investigate the behavior of F(R) in the 
limits t -> and t -> t s . When t -> 0, from Eq. (IIII.124I) we find 



607(207- 1) 



R 



(III.132) 



In this limit, it follows from Eqs. (|III.100p . (|III.124|) , (|III.130|i . (|III.131|> and piI.132|) that the form of F(R) is given 
by 



F(R) 



±J60 1 (20 1 -1)R- 1 



/2 



[^V607(207-l)i?- 



1/2 



27+1 



R 



E 

3 = ± 



^ 57_1 ~^^[60 7 (207-l)] ft/:2 /? ' " 



V 207 - 1 



(III.133) 



On the other hand, when t — > t s , from Eq. (|III.126[) we obtain 



t~t a - 3\ 



R ' 



(III.134) 



In this limit, it follows from Eqs. (|III,100p . pil.1261) . (|III.130|1 . (jIII.131|> and flIII.134[> that the form of F(R) is given 
by 

( \±[t s -3vmR-v 2 ]Y\ 5 

f(R) - — t^tt RJ2pj [t. - sVuor-^ 



v- 

x/l 



1 _ 3^140^-1/2 



0i 



[t s + (P 3 - 15) R- 1 ' 2 



t, - 3V140i?- 1 /2 



(III.135) 
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The above modified gravity may be considered as an approximated form of a more realistic and viable theory. For 
large R, namely, t 2 R 1, the expression of F(R) in (IIII. 135|) can be written approximately as 



F{R) » I 



3VT40(27+ 1) V*o 



\t 5 s R 7 ' 2 . (III.136) 



Note also that viable models that lead to a Big Rip singularity as is discussed in the second chapter necessarily 
predict the phantom divide crossing. This gives exact (non-asymptotic) and realistic modified gravities that show such 
phenomena. Moreover, some of the viable models with accelerating solutions that are asymptotically non-singular or 
tend to softer Type II, Type III and Type IV singularities may also show phantom divide crossing. However, in this 
case the phantom era turns out to be transient. 

3. The reconstruction using e-foldmg 

Let us demonstrate the reconstruction of F(R) gravity in terms of e-folding. Moreover, unlike the previous sub- 
section, the slightly different approach is developed without the introduction of an auxiliary scalar. We introduce the 
e-folding instead of the cosmological time t as: N = In ^L. The variable N is related to the redshift z by 

C - N = ^ = l + z. (III.137) 
a 

Because ^ = H-^ and, therefore, ^ = H 2 -£^ + H^-^, one can rewrite T[) by 

= + 3(H 2 + HH') F'(R) - 18 (aH 3 H' + H 2 (H'f + H 3 H"^j F"(R) + K 2 Pmattcr . (111.138) 

Here, H' = dH/dN and H" = d 2 H/dN 2 . As usual, the matter energy density p ma tter is given by the sum of the fluid 
densities with the constant EoS parameter Wi 

E-3(l+tUi) -3(1 + 10;) -3(l + Wi )N /TTT 1 Qn\ 

Pioa y ^ =/ ,PiQ% e • (III. 139) 

i i 

Let the Hubble rate be given in terms of N via the function g(N) as 

H = g{N) = g (- In (1 + z)) . (III.140) 

Then, the scalar curvature takes the form: R = 6g' (N)g(N) + 12g(N) 2 , which can be solved with respect to N as 
N = N(R). Using (IIII.139|) and (IIII. 140 j) . one can rewrite piI.138j) as 



d 2 F(R) 
dR 2 



= -18 (4 5 (N (R)) 3 g' (N (R)) + g (N (R)) 2 g' (N (R)f + g (N (R)f g" (N (R))] 

+3 (g (N (R)f + g' (N (R)) g (N (R)j) ^ - £H + K 2 ^ p^^e^^ , (III.141) 

i 

which constitutes a differential equation for F(R), where the variable is the scalar curvature R. Instead of g, if 
G(N) = g (N) = H 2 is used, the expression (|III.141[) can be simplified as: 

= -9G (N (R)) (AG' (N (R)) + G" (N (R))) ^fjf ) + ( 3G (N (R)) + \g' (N (R))^ dF ^ R) 



dR 2 V 2 / dR 

+ K 2 J2 p i0 ao 3{1+Wi) e-^ 1+ ^ N ^ . (III.142) 

i 

Note that R = 3G'(N) + 12G(N). Thus, when we find an F(R) that satisfies the differential equation (|III.141|) or 
(IIII. 142|) . and such F(R) theory yields a solution (|III.140p . Thus, such reconstructed F(R) gravity realizes the above 
cosmological solution. It is interesting that the above scheme may be applied not only to general relativity but also 
to modified gravity (partial reconstruction). In other words, one can start from modified gravity which consistently 
describes the entire sequence of the universe evolution eras and satisfies local tests; however, some of its predictions 
seem to contradict the observational data. In this case, one can apply the partial reconstruction of such a model (for 
instance, slightly modifying it at the very early universe) so that the above contradiction is resolved but all good 
properties of the theory are preserved. This provides a way to avoid any future exclusion of a non-realistic modified 
gravity to arrive at more narrow and more realistic sub-class of such theories. 
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4- Stability of the cosmological solution 



In this section, we investigate the condition of the stability of the given cosmological solution. The FRW equation 
pi.7l) can be rewritten by using G = H 2 and the e- folding TV as follows piI.142[) : 







-9G(A) (4G'(A) + G"(A)) 



d 2 F(R) 



dR 2 



R=3G'(N) + 12G(N) 



+ ! 3G(A) + ^G'(A) 



dF(R) 



dR 



F(R) 



+ K 2 p mat tcr (A) 



R=3G'(N)+12G(N) 

(III.143) 



R=3G'(N) + 12G(N) 



Assume that the A-dependence of p ma ttcr is known as in tllll . 139|) . Let a solution of 143[) be G = Gq(N), and 
consider the perturbation from the following solution: 



Then, Eq. (|III.143|) gives 
d 2 F{R) 



= Go (AO 



dR 2 



3G (N) - -G (A) 



R=3G' {N) + 12G (N) 

d 2 F(R) 



G(N) = G (A) + 6G(N) . 



5G"(N) + { 3G (A) (4G' (A) + G '(A)) 



(in. 144) 



d 3 F{R) 



dR 3 



R=3G' (N) + 12G (N) 



dR 2 

12G (A) (4G (A) + G '(A)) 



R=3G' (N) + 12G (N) _ 

d 3 F(R) 



5G'(N) 



(-4G (A) + 2G (A)+G '(A)) 



dR 3 
d 2 F(R) 



R=3G' (N) + 12G (N) 



dR 2 



R=3G' (N) + 12G (N) 



1 dF(R) 

3 dR 



5G.(III.145) 



R=3G' (N) + 12G (N) _ 



Note that in this formulation, as it is assumed that the TV-dependence of /3 ma tter is known, we do not need to consider 
the fluctuation of p m attcr- Using the cosmological time t instead of A, we need to include the fluctuation of p m atter- 
Then, the conditions of the stability are given by 



6(4G (A) + G '(A)) 

+6-^>0 
G (AO 

36(4G o (A) + G o '(A0) 
, 3G '(A) , 1 



d 3 F(R) 



dR 3 



d 2 F(R) 



R=3G' (N) + 12G (N) 



dR 2 



R=3G' (N) + 12G (N) / 



(III.146) 



d 3 F(R) 



dR 3 
dF(R) 



d 2 F(R) 



R=3G' (N) + 12G (N) 



dR 2 



12 



R=3G^(iV)+12G (iV) / 



6G (A) 
Go (A) 



Go (A) G (A) dR 



d 2 F{R) 



R=3G' (N)+12G (N) 



dR 2 



> 0. 



(III.147) 



R=3G' (N) + 12G (N) / 



In the case of de Sitter space, where H and, therefore, Go and R = Rq = 12Go are constants, Eq. ()III.146|) gives 6 > 
and is trivially satisfied. Eq. (|III.147|) becomes 



12G 



dF(R) 



dR 



d 2 F(R) 



R=12G 



dR 2 



-Rq 



dF(R) 



R=12G Q 



dR 



d 2 F(R) 



R=Ro 



dR 2 



> 0, 



(III.148) 



R=R 



which is the standard result. 

One may consider the case that F(R) gravity admits two de Sitter solutions. If one is stable but another is unstable, 
there could be a solution in which the unstable de Sitter solution could transit to the stable solution. If the Hubble 
rate H in the unstable solution is much larger than the Hubble rate in the stable one, the unstable solution may 
correspond to the inflation and the stable one to the late-time acceleration. Note that we may directly construct an 
F(R) theory that admits the transition from an asymptotic de Sitter universe with a large Hubble rate to another 
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asymptotic de Sitter universe with a small Hubble rate. Then, if the solution and F(R) satisfy the conditions piI.146[) 
and (|III.147[) . such a transition occurs. 

Consider the stability of a de Sitter universe satisfying the condition (|II.9j) . or equivalently (|II.37|) . Eq. (|II.9j) can 
be rewritten as 

Let R — R be a solution of (|III.149|> . Then, F(R) has the following form: 

F(R) 



R 2 



= fa + f(R) (R ~ Ro) n ■ (III.150) 



Here, /o is a constant, which should be positive if we require F(R) > and n to be integers greater than or equal to 
two: n > 2. Assume that the function f(R) does not vanish at R = i? , f(Ro) ^ 0. When n = 2, 

F'(Rp) f(R )A 
- Ra+ F^R -)--f + f(A )- (IIL151) 

Then, Eq. (|III.148[) shows that the de Sitter solution is stable if 

-/ </(A )<0. (III.152) 

When n > 3, one finds 

- i?0 + lw = - (IIL153) 

It is necessary to conduct a more detailed investigation to check the stability. Using the expression of m\ in (|II.30|) , 
one can investigate the sign of w? a in the region R ~ Rq. Note that the expression (|II.38|) is not used because the 
expression is only valid at the point R = Rq. Thus, one gets 



ml ~ - H^Wo) {R Ra)n -, (in 154) 

4/( 



Eq. (IIII. 154|) indicates that when n is an even integer, the de Sitter solution is stable if /(i?o) < but unstable if 
f(Ro) > 0. On the other hand, when n is an odd integer, the de Sitter solution is always unstable. Note, however, 
that when f(Ro) < (f(R ) > 0), we find m 2 a > (ml < 0) if R > R but m 2 a < (ml > 0) if R < R . Therefore, 
when /(i?o) <0, R becomes smaller, but when /(i?o) > 0, R becomes larger. The stability condition may be used to 
obtain realistic (unstable) de Sitter inflation for specific F(R) gravity. 



5. FRW cosmology in F(R) gravity with a Lagrange multiplier 

Here, we consider F(R) gravity with a Lagrange multiplier field, following ref. [l36] . In usual F(R) gravity, the scalar 
mode called scalaron appears, which often affects Newton's law. In this section, we try to suppress the propagation 
of the scalaron by imposing the constraint under the Lagrange multiplier field. As a result, however, a propagating 
mode in the Lagrange multiplier field seems to appear, which may break Newton's law but in an easier way. The 
solution of this question may require additional modifications of the constraints. Another purpose of this section is 
the reconstruction. In usual F(R) gravity, one needs to solve the complicated differential equation developed in the 
previous sections to realize the reconstruction. It is shown that the reconstruction can be realized more easily in the 
model with the Lagrange multiplier field. 

The starting action is given by 

S = J ^x^S^FiiR)- \(^d fl Rd^R + F 2 (R)^ . (III.155) 
Here, A is the Lagrange multiplier field, which gives a constraint of 



-d a Rd> 1 R + F 2 {R) = 0. 



(III.156) 
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On the other hand, by the variation over the metric g^ v , the equation of motion follows: 

= \g^Fi{R) + \d^Rd p R + {-R^ + V M V„ - .g^V 2 ) {F{(R) - XF^R) - V (AV M i?)) . (III.157) 
If the Ricci tensor is covariantly constant and the scalar curvature is a constant, 

R»u = , R = R , (III.158) 

Eqs. piI.156|) and (|III.157|> reduce to 

= F 2 (R ), (III.159) 

= F^Rq) - ii? (F[(R ) - \F^R )) . (111.100) 

If Eq. (|III.159[) has a solution, Eq. (IIII. 160[) can be solved with respect to the Lagrange multiplier field: 

A Fi(Ro) ■ (IIL161) 

Then, if Rq is positive, the above solution describes de Sitter space-time, which may correspond to dark energy or 
the inflationary epoch. The spatially flat FRW metric Eqs. (|III.150|) and (/i, v) — (0, 0)-component of (|III.157|) have 
the following form: 

= -^R 2 + F 2 (R), (111.102) 
= --Fx(R) + 18A (h + AHH} 2 

+ {3 (H + H>) - 3H±} [FUR) - XF-(H) + (| + 3H) (»«) } . (in.163) 

When F 2 (R) > 0, Eq. (|III.102j) results in 



R dR 



which can be solved with respect to R as a function of t R — Fn(t). Because 

dH 



(111.104) 



R = 6^- + 12H 2 , (111.105) 
one can find the behavior of H = ^ by solving the differential equation 

0^ + 12H 2 = F a (t) . (IIL166) 
dt 

Using the obtained solution for H = H(t) (and R = Fu(t)), Eq. (|III.103|) becomes a differential equation of the 
multiplier field A, and the behavior of A: A = X(t) follows. (Note that FRW cosmology in scalar theory with a 
Lagrange multiplier constraint was studied in refs. [137| ) . 

Conversely, when the behavior of H{t) is known from the observational data, one may reconstruct F 2 {R) to reproduce 
H(t) by using (jIII.1021) . H{t) gives the behavior of R as R — R(t), which can be solved with respect to t as t = t(R). 
Using (|III.102[) . the explicit form of F 2 (R) is found to be 



1 fdR\ 2 



(III.167) 

t=t(R) 



Note that Fi(R) can be an arbitrary function. Thus, the cosmological reconstruction of the model can be performed 
more easily than that in the usual F(R) gravity. 
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As an explicit example, one may consider 



where hp is a constant. Here, 



H(t) 



ho 
t 



(III.168) 



Therefore, we find 



which gives 



Another example is given by 



R 



dR 
lit 



-6hn 



\2h% 



t 2 



or t = 



-6h + 12/ig 



R 



Here, R± and ui are constants. When t - 
Sitter. One may identify the epoch of t 



12 (-h + 2h 2 ) 



F 2 (R) 



R 



2Ri 



^/6(-h + 2h 2 ) 



7? :! 



12 (-ho 



2h 2 ) 



R— / \ R-\- / \ 

— - (1 — tanhwi) H — — (1 + tanhwt) 



(III.169) 



(III.170) 



(III.171) 



(III.172) 



±oo, R — > ±R± and, therefore, the space-time becomes asymptotically de 
— oo as inflation and t — > +oo as late acceleration. Because 



R = 



{R- - R+) w _ (R- - R+) uj 
2 cosh 2 tot 2 



1 - 



(R. 



i?4 



2RY 



(R- - i? 4 



(III.173) 



from Eq. 167|) . one gets 



2 , ,2 



F 2 (R) 



(R- - R+) uj 
8" 



1 - 



+R+- 2RY 
(R- - R+f 



(III.174) 



Thus, the unification of early-time inflation with dark energy epoch is possible also in Lagrange multiplier constrained 
modified gravity. 

Thus, the universe evolution only depends on the constraint equation (|III.162[) and not on Fi(R). Fi(R) can affect 
the correction to Newton's law. In convenient F% (R) cosmology, all of the dynamics are defined by the form of this 
function. With the constraint (jIII.1621) . F±(R) becomes irrelevant. The cosmological dynamics are defined by the 
form of F2(R). In convenient F(R) gravity, the scalar propagating mode appears, which often violates Newton's law. 
One can show that such a scalar does not propagate in the theory under investigation. Indeed, Eq. (|III.157j) contains 
the second derivative of the multiplier field A. We now need to solve the second-order differential equation to find A, 
which can indicate that A could propagate and that a correction to Newton's law can be found. The magnitude of 
the correction depends on the choice of Fi(R) and/or ^(-R)- 

To investigate the regime of Newton's law, Fi(R) is chosen as 



Fi(R) 



R 
2^2 



The matter presence is assumed. Eq. 157|) has the following form: 



= 



2k 2 



z9ij,vR - R/j 



-V - (-R^ + V^V, - 5 ^V 2 ) (XF^(R) - (AV^i?)) 



For the solution with A = 0, Eq. 157|) reduces to the Einstein equation, 



1 fl 



— — ^ [ T,9iivR — Rfiu + Tftv 



(III.175) 



(III.176) 



(III.177) 



70 



Here, T^„ is the matter energy-momentum tensor. Without matter T^ v — 0, the Schwarzschild space-time with 
R = R^v = is a solution, which also satisfies the constraint equation 156|) if F 2 (0) = 0. In the case with matter 
T^v 7^ 0, however, the Einstein equation (|III.177j) gives 

R = ~k 2 T. (III.178) 
Here, T is the trace of the energy-momentum tensor. The constraint equation 156[) is rewritten to 

K 4 

= —d^Td^T + F 2 (-k 2 T) , (III.179) 

which is not always satisfied. Thus, in the presence of matter, the constraint equation (|III.156[) should be modified 
to be 

= ±d li Ra"R + F 2 (R) - yfyTflT - F 2 (-k 2 T) . (III.180) 

This indicates that the total constrained action with matter should be, instead of ([III. 155ft , 

R . fl „ „„„ „ „ k 4 



S = Id xyj—g 



^ - A <j ^Rd^R + F 2 (R) - '—d^Td^T - F 2 (-k 2 T) } + C n 



(III.181) 



In this case, Newton's law can be easily reproduced. Of course, some other form of the constraint may also solve this 
problem. 

We also note that the form of the constraint could be correct when i*\ (i?) is given by (|III.175I) . For general form 
of Fi (R) , the constraint should be changed. 

Because T vanishes in the vacuum, the constraint (|III.180[) reduces to 



= -d^Rd^R + F 2 (R) ~ F 2 (0) . (III.182) 



If F 2 (0) = 0, the constraint pil,180[) gives (|III.156|) and (IIII.162I) and the cosmological evolution can be generated. 
Note that F 2 (R) piI.171|) satisfies the condition F 2 (0) = but F 2 (R) in (|IH.15ip does not. When the condition 
F 2 (0) = is satisfied, there are two classes of solutions for the constraint 162|) . One is a trivial solution R = 0, 
and another corresponds to the non-trivial cosmological evolution given by (|III. 164ft . Near solar systems and galaxies, 
the solution could correspond to the trivial one with R = to reproduce Newton's law, but at large scales, the solution 
should correspond to (|III.164[) . so that the evolution of the universe can be generated. It is not so trivial to prove 
or to deny that the two solutions are connected in the intermediate region between the galaxy scales and large-scale 
universe. More careful (possibly numerical) analysis should be done in order to solve this problem. 



D. F(Q) gravity 

The cosmological reconstruction scheme developed above can be applied to an arbitrary modified gravity. As one 
more example of its application, we discuss the reconstruction in modified Gauss-Bonnet gravity. The action (III. 123ft 
can be rewritten by using the auxiliary scalar field <f>: 



d x-\/—g 



R 

2k 2 



By the variation over the scalar field <p, an algebraic equation is obtained: 



(III.183) 



(III.184) 



It can be solved, at least locally, with respect to (f> as 4> = 4>{Q). Inserting the obtained expression of 4>(G) into the 
action 183ft . we find the action (|II. 123ft : 



S = I d A x^—g 



(III.185) 



Here. 



f(Q) = -V{<t>{Q)) + ti(<j>{Q))G. 



(III.186) 
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In the following, we use the results obtained in refs. |52H54j and consider the reconstruction of scalar-F(CJ) gravity. 
By the variation over the metric in the FRW universe (|II.6p . one obtains the following equations: 

= *H* + VU,)+ 24fl**^, 
K at 

= -1 f 2H + 3H 2 ) - V{4>) - M^m)) _ mH H^m - 16H sdm)} . (IH.187) 
k V / dt 1 at at 

Using (|III.187|) . it follows that 



Vm)) = ^H(t) 2 - 3a(t)H(t)W(t) , W(t) = - / -^rH(h) . (III.188) 

Because there is no kinetic term of cf>, one can redefine <f> properly and even identify <j> with the cosmological time: 
4> = t. In the same way as in the previous section, if we consider the following V(<j)) and £ given in term of a single 
function g, 

V(<P) = ^g' (0) 2 -3</ (<p)e^U(<j>), 



8 7 r g'(<t>i) 

the explicit solution follows: 



i r<t> e g{4>i) 9 

o/ *hzKTv U (fa) , U(0) = -j dcPxe-^g"^) , (III.189) 



a = a e ff(t) (i? = g'(t)) . (III.190) 



In terms of e- folding N 



1 r N , e N i 



V(j>{N)) = -^H(N) 2 - 3e N H(N)W(N) , W(N) = / ^^-H(Ni) . (III.191) 



fx 



N 



dNi 



Now we identify <j> with the e-folding N: cf> — N instead of 4> = t- Then, if V(<j>) and £ are given by 

V{4>) = \h{4>f -2>h(<t>) e*U{cj>), 



£(<£) = -/ #iT7TT3^i) . U(cf>) = - I dfo.e-*ti fa) , (111.192) 



e^ 1 , , „, 2 

we obtain the solution H = h(N). This general scheme may now be applied to explicitly reconstruct the modified 
Gauss-Bonnet gravity by realizing any given FRW metric with a particular scale factor. 



E. Scalar-Einstein-Gauss-Bonnet gravity 

In the same way as above, let us reconstruct string-inspired scalar-Einstein-Gauss-Bonnet gravity in III CI The 
reconstructed theory realizes the solution with the arbitrary Hubble rate H . Combining (III.173[) and (III. 174|) and 
deleting V(<p), we obtain 

2 • - 9 
= -^H + c/> 2 



ilt 1 dt dt 

H 2 dmt)) \ 
dt V a dt J 



8 d /^^m _ 
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The following expressions of £ (4>(t)) and V (<f>{t)) may be obtained: 

7(<K*)) = lff(t) 2 - i^(t) a - Mt)H(t) J* ( J^(*0 + <^i) 2 ) • On-194) 

Therefore, for the model where V(0) and £(</>) are given by adequate functions git) and /(</>) as follows, 

«« = I / «■ %,)■ / "WW*"*"*" O' e<«> + tw) ' (nL195) 

the solution of (|II.173|> and pi. 174ft is given by 

<\> = f~\t) (t = /(<£)) , a = a e"W (if = </(*)) . (III.196) 

Similarly, the reconstruction of the model in terms of N , instead of the cosmological time t may be presented. Using 
the same steps that led to Eg. QUI. 1940 . we find 

1 1°^ nN\ pNi j/y / ' n \ 

= s y y ^ (^(^j+^yw'j , 

F(0(iV)) = ±H(N) 2 -~H(N) 2 ct>'(N) 2 -3e N H(N) J ^_ J ff'(7V 1 ) +jff (7V 1 )0'( iVl )2j . (ln .197) 

In terms of functions h(cj>) and /(</>), if V(^>) and £(0) are given by 

3 / n2 h(f ((/>)) , . . / o , . ft' (/((/-l)) 

w = ^(/(0)) -^f - 3 K^) e/w y 



m = o I #i 7 n3 / dfcfiMe-'^ —h' ( /(0a) + v / [ , (III.198) 
h(f(<f>i)) ' V '' ' 

one obtain the following solution: 

= /- 1 (7V) (iV = /(0)), H = h(N). (III.199) 

This may now be applied to obtain an explicit reconstruction of scalar-Gauss-Bonnet gravity, which admits any given 
FRW cosmology as a solution. 

F. F(R) Hofava-Lifshitz gravity 

In this section, we consider the reconstruction of F(R) Hofava-Lifshitz gravity. It is based on the results of ref. [79| . 
To start, let us analyze the simple model F(R) = R, of which the cosmology was studied in 90l - |ll2j |. In such a 
case, the FRW equations are similar to those of general relativity, 

k 2 ■ k 2 

H = g/g^ _ ] \ ^matter j H = — 2^3^ _ j\ (^matter + Pmattcr) , (III. 200) 

where, for A — > 1, the standard FRW equations are recovered. Note that the constant /i is now irrelevant because, 
as pointed out above, the term in front of /i in (|II.274I) becomes a total derivative. For such a theory, one has 
to introduce a dark energy source as well as an inflaton field, to reproduce the cosmic and inflationary accelerated 
epochs, respectively. It is also important to note that, for this case, the coupling constant is restricted to be A > 1/3; 
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otherwise, Eqs. (|III.200|) become inconsistent. It seems reasonable to think that, for the current epoch, where R has a 
small value, the IR limit of the theory is satisfied by A ~ 1, but for the inflationary epoch, when the scalar curvature R 
goes to infinity, A will take a different value. It has been realized that for A = 1/3, the theory develops an anisotropic 
Weyl invariance (see [Hj]), and thus it takes a special role, although for the present model this value is not allowed. 

We now discuss some cosmological solutions of F(R) Hofava-Lifshitz gravity. The first FRW equation, given by 
pi.278j) with C = 0, can be rewritten as a function of the number of e- foldings N = In instead of the usual 

time t. Here, we extend the reconstruction formalism to the Hofava-Lifshitz F(R) gravity. Because 4? — H-Ar and 



<'o 

'■'nil i'\ 

dt ~ A± dN 

H2 7m+ H 7MJW' the first FRW equation (|II.278|) is rewritten as 



dN 2 " dN 2 ' " dN dN 

= F(R) - 6 [(1 - 3A + 3p)H 2 + pHH'] dF ^ + 36pH 2 [(1 - 3A + 6fj,)HH' + pH' 2 + pH"H] d ^B- - /^matter , 

dR d 2 R 

(III.201) 

where the primes denote derivatives with respect to N. Thus, in this case, there is no restriction on the values of A or 
p. By using the energy conservation equation and assuming a perfect fluid with an EoS of 

Pmattcr ^matter Pmattcr ; 

the energy density yields 

Pmattor — POO ~ P0 a e • (lll.ZUZ) 

As the Hubble parameter can be written as a function of the number of e-foldings, H = H(N), the scalar curvature 
in (|II.274j) takes the fo mi 

R = 3(1 - 3A + 6^)H 2 + GfiHH' , (III.203) 

which can be solved with respect to N as N = N(R), and one obtains an expression (|III.201[) that gives an equation 
of F(R) with the variable R. This can be simplified by writing G(N) = H 2 instead of the Hubble parameter. In such 
a case, the differential equation (|III.201|) yields 

= F(R) - 6 [(1 - 3A + 3^)G + + 18/* [(1 - 3A + 6/t)GG' + »GG"] - K 2 p a~ 3{1+w) e- 3 ^ N , 

2 J dR d 2 R 

(III.204) 

and the scalar curvature is now written as R — 3(1 — 3A + 6/i)G + 3/iG'. Thus, for a given cosmological solution 
H 2 = G(N), one can resolve Eq. (|III.204|) . and the F(R) that reproduces such solution is obtained. 

As an example, we consider the Hubble parameter that reproduces the ACDM epoch. It is expressed as 

H 2 = G(N) = H 2 + ^ Po a- 3 = H 2 + ^ Po a^e- 3N , (III.205) 

where Hq and po are constant. In general relativity, the terms on the r.h.s. of Eq. (IIII.205I) correspond to an effective 
cosmological constant A = 3Hq and to cold dark matter with an EoS parameter of w — 0. The corresponding F(R) 
can be reconstructed by following the same steps as described above. Using the expression for the scalar curvature 
R = 3(1 — 3A + 6/z)G + 3/xG', the relation between R and iV is obtained, 

where k = ^poCIq 3 . Then, substituting (IIII.205I) and (|III.206[> into Eq. (|III.204|) . one obtains the differential expression 

= (1 - 3A + 3p)F(R) - 2 ( 1 - 3A + ^ J R + 9/x(l - 3\)H 2 ^S- 

V 2 ) dR 

d 2 F( R) 

-6pi(R - 9pH 2 )(R - 3H 2 (1 -3X + 6/t)) \U - R - 3(1 - 3A + 6p)H^ , (III.207) 

d 2 R 

where, for simplicity, we have considered a pressureless fluid w = in Eq. (IIII.204[) . Performing the change of variable 
the homogeneous part of Eq. pil.207|) can be easily identified as a hypergeometric differential 



3H 2 (l+3(p-\)) 

equation 



d 2 F dF 
= x(l - x)-^z + (7 - (a + + l)x) — - apF , (III.208) 
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with the set of parameters (a, ft, 7) being given by 

The complete solution of Eq. (jIII.208|) is a Gauss hypergeometric function plus a linear term and a cosmological 
constant from the particular solution of Eq. (|III.207j) . namely 

F(R) = CiF(a, p, 7; x) + C 2 x 1 - 1 F{a - 7 + 1, /3 - 7 + 1, 2 - 7; x) + —Ft - 2A , (IIL210) 

where C\ and C2 are constants, ki = 3A — 1 and A = — 3 ^°i]^x +tu,) \ ■ Note that for the metric (|III.205|) . the 
convenient F(R) gravity was reconstructed and studied in refs. [H, Il20l Il2ll Il38l Il39j . In this case, the solution 
(IIII.210P behaves similarly to the convenient F(R) theory, except that now the parameters of the theory depend on 
(A, /z), which are allowed to vary, as was noted above. Thus, the cosmic evolution described by the Hubble parameter 
(|III.205j) is reproduced by this class of theories. 

One can also explore the solution pil.205[) for a particular choice of the parameters /i = A — |, which plays a 
special role in the cosmological expansion as is shown below. In this case, the scalar R turns out to be a constant, 
and Eq. (|III.207j) . in the presence of a pressureless fluid, has the solution 

F{R) = — R-2A, with A = ^(3A-l)i^. (IIL211) 

Thus, for this constraint on the parameters, the only consistent solution reduces to the Hofava linear theory with a 
cosmological constant. 

As a further example, we consider the phantom (super)accelerating expansion. Such a system can be easily expressed 
in general relativity, where the FRW equation reads as H = ^Pph- Here, the subscript "ph" denotes the phantom 
nature of the fluid, which has an EoS given by p p h = w p hP P h with w p h < —1. By using the energy conservation 
equation, the solution for the Hubble parameter turns out to be 

H(t) - 77^7 . ( IIL212 ) 

where Hq = —1/3(1 + w p h), and t s is the Rip time. As in the above example, one can rewrite the Hubble parameter 
as a function of the number of e-foldings; this yields 

G(N) = H 2 (N) = H%c 2N/H " . (III.213) 

Then, by using the expression of the scalar curvature, the relation between R and N is given by 

e 2N/H _ R fm 2U) 

H (AH + 6»)- (m - 2U) 
By inserting (|III.213|) and (|III.214|) into the differential equation (|III.204|) . we get (without matter) 

R^flSp-+k 1 R^B- + k F(R) = 0, (III.215) 
dR? dR 



where 



6fi(AH a + 12/i) ' 1 12p{AH + 12 n) v ; 



Eq. (|III.215[) is an Euler equation, whose solution is well-known 



F(R) = C 1 R m + +C 2 R m - , where m± = 1 ^ - ^| l? ^ . (III.217) 

Such a theory belongs to the class of models with positive and negative powers of the curvature introduced in ref. [2| as 
is discussed in the second chapter. Thus, an F(R) Hofava-Lifshitz gravity has been reconstructed that reproduces the 
phantom dark epoch with no need of any exotic fluid. In the same way, any given cosmology may be reconstructed, 
including the unified inflation within dark energy era. 
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G. Non-minimal Yang-Mills theory 



We now develop a reconstruction scheme of the Yang-Mills theory by starting with the following action 

R 



S = dx 



-fJ 



fJ,V\ 



(III.218) 



where = d^A^ — d^A 1 ^ + f abc A b A^, and J" may be assumed to be a conti nuously differentiable function. 
Here, we review the reconstruction of the modei piI.218|) based on ref. [l40j . 



We concentrate here on the SU(2) case where / = e . Because 



» \ r Hv r ) _ _. hbc A b p c~tP 

8A h ~ ~< ' 

5(F°F a i M ') , „ 



S (d a A h ^ 

the equation of motion for the field potential A a turns into 

SS 



(III.219) 
(III.220) 



SS 
6A* 



= 0. 



and, from here, 

d v [V^gF' (F^F aaf} ) F auf *] + J~gT' (F^F aa/3 ) e abc A b l/ F cv ' 1 = . 
The variation of pil.218|) with respect to g^" yields the following equation of motion 



2k 2 



(III.221) 



(III.222) 



(III.223) 



where the following relation is used " - = 2J 7 ' (f^F 110 ^ F^F^ . Considering now a FRW universe (1TL61) . 

and the following Ansatz for the gauge field, 



A" 



0, n = 0, 



the [i = component of (|III.221[) becomes an identity, the fi = i component yields 



di 



a{t)F' (F^ p F aal3 ) A(i)e A « + ^£L T' (F% F aa P) e 3A ^ =0, 

J a(t) 



and the (i, t) component of (|III.223[) is 

3H 2 (t) 1 



2k 2 ' 2 

and the component of (|III.223|) reduces to 



~ T (F^F aafS ) + 2a 2 T' (F£ F aa ^ 



A 2 (i> 



s 2A(t) 



1 

2^ 



2H(t) + 3H 2 (t)\ - ~JF {F^F aa ?) - 2a 2 ^y^' (F^F aal3 ) 



=2A(t) 



a 2 (t) 



X 2 (t)-2a 2 



g 2A(t) 



= 0. 



Adding (|III.226[) to pil.227|i . one arrives at 

a 4 (i)i?W -4k 2 a 4 F'(F^F aa P) e 4X ^ = , 

and then 

a\t)H(t) 

„-4A(t) 



F'(F^F a 



4n 2 a 4 



(III.224) 



(III.225) 



(III.226) 



(III.227) 



(III.228) 



(III.229) 
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Using pil.229p . Eq. (|III.225|) reduces to: 

2a 2 H{t)e 2m + ha(t) a(t) H(t) + a 2 (t)H{t)] \(t) - 'Sa 2 (t)H(t)\ 2 (t) + a 2 (t) H (t)X(t) = 0, (III.230) 

which constitutes a differential equation for X(t). Thus, by using Eq. (|III.224[) . once we have the function A(i), given 
by (IIII.230|) . we can obtain the corresponding Yang-Mills theory that reproduces the selected cosmology. The Ansatz 
considered above actually leads to a mathematical solution of the problem. 

As an example, we consider the case of a power-law expansion: a(t) — (j^j ? where ti and hi are constant, and 

assuming \{t) = (hi — 1) In (j^J + Ai, where Ai is again a constant, Eq. (|III.230| reduces to the following algebraic 
equation 

hi(hi - 1) + a 2 t\e 2Xl = 0, (III.231) 



thus, 



and 



A(() = (ftl _ 1)ln (l) + i ln (M^)). (IIL233) 

With the help of this reconstruction scheme, the function A(i) given by (jIII.2331) is obtained. Using pil.224[) . one is 

able to reproduce the cosmology given by the power- law expansion: a(t) — (^j^j 

In summary, a cosmological reconstruction scheme was developed in the current chapter. It was presented in a very 
general form and was applied to a number of modified gravities introduced in the second chapter. Several examples 
of important cosmological epochs are considered. The cosmological reconstruction of modified gravities that realize 
such epochs as solutions was successfully made. Despite the fact that we limited ourselves to only a spatially flat 
FRW universe, this scheme may be easily used for the generation of other background solutions, i.e., an anisotropic 
universe, black holes, wormholes, etc. 



IV. FINITE-TIME FUTURE SINGULARITIES IN MODIFIED GRAVITY 



This chapter is devoted to the study of finite-time future singularities, which often occur in the current dark 
energy models. It is demonstrated that this is a quite general phenomenon; the effective phantom and some effective 
quintessence dark energy models including fluid, scalar, and modified models may show such behavior in the finite 
future. We start from a simple dark fluid universe and demonstrate that coupling with dark matter cannot cure such 
future singularities. The introduction of a specific form of modified gravity may naturally solve the problem of future 
singularities. The occurrence of singularities in modified gravity is considered in detail, as is the possible avoidance 
of the singularities via the addition of extra gravitational terms. 



A. Future singularities in a dark, fluid universe: dark matter versus modified gravity 

Let us reiterate several simple facts about coupled phantom dark energy. First, we assume that the energy density 
Pde and the pressure pde of the dark matter satisfy the standard conservation law 144|) . Let the EoS parameter 
wde: wde = Pde/ Pde < — 1. The phantom dark energy solution of the FRW equation (jll.3[) (after replacing p ma tter 
and pmattor with p DE and p D E, respectively) is given by 

2 

H = 3(1+ " DE) , (IV.l) 
t s t 

with the famous Big Rip singularity, at t = t s . 

Following ref. [47j . the model where phantom dark energy couples with dark matter is considered. Then, the 
conservation law is modified as 



Pde + 3H (1 + w DE ) Pde = 



—Qpde , Pdm + 3-ffpDM = Qpde ■ 



(IV.2) 
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Here, Q is assumed to be a constant. The first equation can be solved as 

Pde = PDE(o)a" 3(1+ " DE) e- Qt ■ (IV.3) 
Here, Pde(o) is a constant of the integration. The second equation (|IV.2[) gives 



Pdm = Qa{ty 3 I dt'pvE^a-^e-Q 1 . (IV A) 



Thus, the second FRW equation is 



K 

An exact solution of (|IV.5|) is the de Sitter space 



4 (2H + 3H 2 ) = PDE = wdePde = wp D E(o)a- 3{1+WOE] e- Qt . (IV.5) 



Q 



where ao is given by 



a(t) = a e H = --^ r , (IV, 

3(1 + w D e)/ 



^ ( q ^ -3(i+tu DE ) friT „x 



n 2 \3 (1 + w D e) 

Note that because we are considering the phantom dark energy with u>de < — 1, if (IIV.6|) is positive and Eq. PV.7P 
has a real solution. 

Thus, the coupling of the dark matter with the phantom dark energy may give the de Sitter solution instead of 
the Big Rip solution (IIV.1|) . This does not always mean that the Big Rip singularity can be avoided, but it opens a 
possibility that the universe could evolve into a de Sitter universe instead of the Big Rip. 

Let us show that the above construction may help to solve the coincidence problem. One may identify the Hubble 
rate H with the present value of the Hubble rate H = — 3 ( 1+ ^, DE ) = Hq ~ 10~ 33 eV. Eq. (|IV.3|) shows that the dark 
energy density is a constant 

— 3(l+lO DE ) /TX T Q\ 

PDE = /0 DE (o)a • ( IV - 8 ) 

Eq. (|IV.4|) can be integrated as 

Pdm = PDMofl -3 - (1 + wde) /ODE(o)a(7 3<1+ ™ DB) ■ (IV.9) 
Here, pdmo is a constant of integration, but the first FRW equation 

3 

-^H 2 = p DM + PDE , (IV.10) 

shows pdmo — and, therefore, the dark matter density pdm is also constant: 

PDM = - (1 + WDE) PDE(0) a 3(1+U,DE) = - (1 + w De) PDE ■ (IV. 11) 

Then, if the de Sitter solution (|IV.6j) is an attractor, by choosing 

1 4 
- (1 + wde) ~ g , i-e. wde ~ ~ , (IV.12) 

the coincidence problem could be solved. That is, even starting with a wide range of initial conditions, the solution 
approaches the de Sitter solution, in which the ratio of the dark energy and the dark matter is approximately 1/3 
and is almost independent of the initial condition. 

To check whether or not the de Sitter solution pV.6|) is an attractor, the following perturbation is considered 

-t+A(i) 



a = a e 3 ( 1+ »™) . (IV.13) 

Here, A(t) is assumed to be small. The second FRW equation (IIV.5|) gives 

1 foA 2( ? A^ Q/i , <\ -3(1+u>de)a 3(1 + w D e)/ Q V. n\nA\ 

~— 2A_ 77-r rA =-3(l + u;DE)wpDE(o)a ' ; A = - ttt— H A , (IV.14) 

k 2 V (1 + ^de) / k \3(1 + wde)/ 
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which is a very simple linear differential equation with a constant coefficient. Here, Eq. (IIV.7|) is used. Assuming 
A = c At , it follows that 



= A 2 -— ^— A + 3(1 + ™ DE )(-— 9. ) (iv.ir,) 

1 + ^DE \3 (1 + Wde), 



that is, 



a = a ± = , - , ± - { r_«_y _ in±^i f 7 g _vr . (Iv .i6) 

2(l + w DE ) 2|V 1 + W W 3 \(1 + wde)/ J 

Then, A_ < but A+ > and, therefore, de Sitter solution (|IV.6I) is not stable. Because A_ < 0, however, the 
solution is saddle point and, therefore, with the appropriate initial condition, there is a solution that approaches the 
saddle point de Sitter solution ()IV.6|) . Although we need to study the global structure of the space of the solutions, 
such an initial condition could correspond to the direction toward A_. Thus, generally speaking, the coupling of the 
dark energy with dark matter does not prohibit the existence of the (Big Rip) singular solution. If the singularity 
corresponds to the stable solution, because the de Sitter space is not completely stable, the solution will finally evolve 
into the singular solution. Thus, the appropriate choice of initial conditions may help to realize the non-singular de 
Sitter cosmology, which also solves the coincidence problem. At least, one may expect that the future singularity 
would occur at more later times compared with the case in which the dark mat ter d oes not couple with dark energy. 
Let us consider the simple example of a perfect fluid with the following EoS |l4l| : 



p=-p + A P a , (IV.17) 

with a constant A and a. The EoS (|IV.17|) is a special case of (|II.339| ). In the spatially flat FRW space-time (|II.6|1 . 
the Hubble rate is found to be 



( 3 



H 



A 



when a = 1 , A > 
when a = 1 , A < 



Be-^, whence i A<0 ( IV - 18 ) 



2 



Ct l/(l-2a) 

{ or(7(t s -*) 1/(1 - 2a) whena ^ 



Here, B, C , and C are positive constants. Now, one can describe the future, finite-t ime singularities of the universe 
filled with the above dark fluid for different choices of theory parameters (see ref. [122]). When a < 0, a Type II 
or sudden future singularity occurs. When < a < 1/2 and 1/(1 — 2a) is not an integer, a Type IV singularity 
occurs. When a = 0, there is no singularity. When 1/2 < a < 1 or a = 1 and A < 0, a Type I or Big Rip-type 
singularity occurs (see also the last section of the second chapter). When a > 1, a Type III singularity occurs (for 
the classification of all four types of future singularities, see [34|). 

In the case of a Type II singularity, where a < 0, H vanishes as H ~ (t s — t) 1 /^ 1-2 ") when t — > t s and, therefore, p 
vanishes as follows from the FRW equation. Then, near the singularity the EoS (1II.339[) is reduced to 

p - Ap a . (IV.19) 

On the other hand, in the case of a Type I singularity, where 1/2 < q < 1 or a = 1 and A < 0, H and, therefore, p 
diverge when t — > t s . Then, the EoS (|II.339|) reduces to 

p porp~-(l-A)p. (IV.20) 

In the case of a Type III singularity, where a > 1, H and p diverge when t — > t s and, therefore, the EoS (|II.339[) 
reduces to 

p - Ap a . (IV.21) 

Let us now consider the dark energy (|IV.17[> coupled with the dark matter as in (|IV.2|) . The conservation law is 
given by 

Pue + 3i/Apg E = -Qpde , Pdm + 3Hp BM = Qpue ■ (IV.22) 
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The solution of pv.22|) is 

Pde (t) = e~ Qt (-3A (I -a) J dt' H{t')e^-^ Qt '\ ' ° , p DM = Qa^ 3 ^ df 'a(t') 3 p DE (t') . (IV.23) 
If A < (and Q > 0), there is a de Sitter solution = Hq with a constant Hq 



and pde and pbm are constants 



P ™ = {—Q-) ' PDM= 3H {— Q-) ■ (IV ' 25) 

This demonstrates that if Hq ~ Q, we find pdm/Vde ~ 1/3, and the coincidence problem may be solved. If Hq = Q, 
Eq. (|IV.24[) determines the value of A: 

•^-KsK)'"- (iv ' 26) 

One can now investigate the (in)stability of the de Sitter solution H = Hq by putting H = Hq + 5H. The 
perturbation of the energy density is 

«5 PDE = -3A f-^jp) 1 " e- (1 ~ a)Qi ^ dt'SHit')^ 1 -^' . (IV.27) 
The second FRW equation becomes 

_ J_ (sH + 6H 5h) = 3A (l + (-^p) ^ /' dt'SHit')^ 1 -^*' . (IV.28) 

By differentiating both sides of (|IV.28|) . one gets 



2 i ( , aQ\( 3AH 



= ^ + {6-ff + (l-a)Q}^+|6if (l-a)Q + 3^ + ^ J ^ j> <Lff 

= (5 J ff + {6i7 + (l-a)Q}(5ij + 3ffo(l-2a)g ( 5i7. (IV.29) 
In the second equality, Eq. pV.24p is used. Assuming SH oc e At , one gets 

= A 2 + {6H a + (1 - a) Q} A + 3H (1 -2a)Q, (IV.30) 

whose solution is given by 

A = A± = _ 6g o + (l-«)Q ± 1 | {giJo + (1 _ a) Q}2 (i _ ^ Q y (iy 31) 

Because i/o, Q > 0, if 

q < \ , (IV.32) 

both A± are real and negative if 

D = (6-ffo + (1 - a) Qf - 12H (1 - 2a) Q > , (IV.33) 
or complex, but the real part is negative if 

D = (6H + (1 - a) Q) 2 - 12H (1 - 2a) Q < . (IV.34) 
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Thus, as long as a < \ pV.32p . the de Sitter solution is stable and, therefore, the singularity can be resolved. On the 
other hand, if 

6H + (l-a)Q >0, 3iJ (l - 2a)Q < 0, (IV.35) 

that is, 

l + ^>«>i, (IV.36) 

we find A+ > and A_ < 0, as at the beginning of this section. If 

6i? + (l-a)Q<0, 3H (1 - 2a) Q < 0, (IV.37) 

that is, 

a > 1 + ^7T ' (IV ' 38) 

we find A± > and the dc Sitter solution is completely unstable. 

Thus, the Type II singularity, where a < 0, and the Type IV singularity, where < a < 1/2 and 1/(1 — 2a) is not 
an integer, can be resolved by the coupling of the dark energy with the dark matter. As mentioned above, even with 
a coupling of the dark matter with the dark energy, there could be a singular solution. The de Sitter solution with 
a < i is, however, at least a local minimum. Then, if the universe started with an appropriate initial condition, the 
universe would evolve into the de Sitter universe (asymptotically evolving into a de Sitter universe). A Type I (Big 
Rip) singularity, where 1/2 < a < 1, and a Type III singularity, where a > 1, could not be removed by the coupling 
of the dark matter with the dark energy. Even if a solution goes near the de Sitter point, the solution could evolve 
into the singular solution if the singular solution is stable. Nevertheless, the ability to avoid some future singularities 
due to the coupling of the dark energy with the dark matter looks quite promising. 

Because the coupling of the dark energy with the dark matter does not always remove the singularity, we now 
consider what kind of fluid could cure the future singularity. In the case of the Big Rip singularity, for example, the 
energy density of the dark energy diverges like pde ~ 1/ (t s — t) 2 ~ R when t — > t s . Here, R is the scalar curvature. 
Then, one becomes interested in a fluid with a positive pressure (and a positive energy density) that grows more 
rapidly than the dark energy pressure. There is no such fluid with a constant EoS parameter. However, one can 
consider a pressure that is proportional to a power of the curvature; for example, 

Pfluid « R 1+t , (IV.39) 

with e > 0. Then, the total EoS parameter becomes greater than —1 for a large curvature, and the Big Rip does not 
occur. 

This kind of inhomogeneous effective fluid [l!8| could be realized by quantum effects (for instance, taking account 
of conformal anomaly) or by modified gravity as was explained in the second chapter. As an example, we consider 
F(R) gravity, where F(R) = R + f(R) behaves as f(R) oc R m . In the case of rn = 2, the power-law solution, like 
the Big Rip singularity is prohibited. It can happen that the (asymptotic) de Sitter solution occurs, instead of the 
power-law solution and, therefore, the singularity would be resolved. Therefore, the introduction of the special form 
of f(R) (precisely, the R 2 term) prevents the future singularity in the universal way. Note that as with other dark 
energies, the modified gravity itself (see next section) may lead to all four possible future singularities, which may 
again be resolved by the addition of the R 2 term 37]. Thus, curing all future singularities of fluid dark energy may 
be done via modifying the gravity. This seems to be quite a fundamental reason to start from modified gravity as a 
gravitational alternative for dark energy. 



B. Finite-time singularities in F(R) gravity 



Let us show that F(R) gravity generates all four known types of finite-time singularities. Moreover, the resolution of 
future singularities may be done again by adding the R 2 term [37J,|40|. This section is based on results of refs. |36ll39l|. 

As the first example, we consider the case of the Big Rip singularity (|IV.1[) . To find the F(R) gravity that generates 
the Big Rip-type singularity, it is convenient to use the reconstruction scheme developed in subsection IIII C ll Let 
h = 3 ( 1+ ^ DE ) ■ Li the case of (|IV.1|) . the general solution of (IIII. 105[) when the matter contribution is neglected, 
Pio = 0, is given by 



(IV.40) 
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when ho > 5 + 2\/6 or ho < 5 — 2-v^ and 
P(0) = (t s - ^)-^ +1 )/ 2 



, ... -/i 2 , + 10/io - 1 
A cos (t s - <j)) In — y — 



j ? sin((t^0)ln- fe ° + m °- 1 



when 5 + 2-\/6 > ho > 5 — 2-\/6- One finds the form of F(R) when i? is large as 

F(R) cx i? 1 ""-/ 2 , 

for the ho > 5 + 2\/6 or /io < 5 — 2\/6 case and 

F(R) oc R^o+^/i x (oscillating parts) , 



for the 5 + 2^/6 > h > 5 — 2y6 case. 

Let us investigate a more general singularity 



H~h (t s - t) 



-P 



(IV.41) 



(IV.42) 



(IV.43) 



(IV.44) 



Here, ho and /3 are constants, ho is assumed to be positive, and t < t s , as it should be for an expanding universe. 
Even for a non-integer /3 < 0, some derivative of and, therefore, the curvature become singular. Assume j3 ^ 1. 
Furthermore, because /3 = corresponds to the de Sitter space, it is proposed that j3 ^ 0. When > 1, the scalar 
curvature i? behaves as 



i? - 12iJ 2 - 12/i 2 (t s - t) 2/3 . 
On the other hand, when fj < 1, the scalar curvature i? behaves as 

R ~ QH ~ Qhof3 {t s - ty' 3 ' 1 . 

If we write 

P(0) = e- feo ^-«"" +1 / 2 ( 1 -/9)5(^), 
Eq. (|III.105p . without the matter contribution, has the following Schrodinger equation like form: 



d 2 S(<p) , {Who ,,-P-i_hl ,. 

+ { 2 ( s 0) 4 [ s 



-26 



s. 



When <f> = t — > ts, in the case that /3 > 1, one finds 



-/9-1 



On the other hand, in the case that /3 < 1, it follows that 



5/3/io 



(«« ~ <t>) 



-/3-1 



> 



ft 2 

^ (t. ~ 0) 



-2/3 



^(ts-^y 20 



In either case, Eq. (|IV.48|) reduces to the following form: 







when <h = t —> t s . Here 



With a further redefinition 



Vo = , a = /3 + 1 when /3 < 1 

Vo = ^ , a = 2/3 when fi > 1 



(IV.45) 
(IV.46) 
(IV.47) 

(IV.48) 

(IV.49) 

(IV.50) 

(IV.51) 

(IV.52) 
(IV.53) 



82 



Eq. (|IV.51[) has the following form: 



dy 2 H 16 4; y 2 (2 -a) 2 ' 



Note that y -> when <j> -> t s if 1 - a/2 > 0, but y -> 00 when ^ ->• i s if 1 - a/2 < 0. Then, if 1 - a/2 > 0, 
Eq. (|IV.54p reduces to the following form when <f> — > t 3 : 

(Pip (a 2 a\ 1 

dy 2 \ 16 4 J y 2 



= ^-(^-t)^- (IV.55) 



whose general solution is given by 

if = Ay 01 ' 4 ' 1 + By- a/i . (IV.56) 

Here, A and B are constants of integration. On the other hand, if 1 — a/2 < 0, Eq. (|IV.54[) reduces to the following 
form when <b — > t s : 



d 2 ^ | Wp 
dy 2 (a - 2) 

When Vq > 0, the general solution of (|IV.57p is given by 



= Ti + • (IV.57) 



y = A cos (wy) + B sin (wt/) , cj = -. (IV. 58) 

a — 2 

Here, A and _B are constants of integration. On the other hand, if Vq < 0, the general solution has the following form 



y = Ac^ + Be-^ , u = 2 ^^° . (TV. 59) 

a — 2 

From the above analysis, one may obtain the asymptotic solution for P when cf> t s . 

• /3 > 1 case: From (|IV.52|1 . a = 2/3 > 2 and, therefore, 1 — a/2 — 1 — (3 < 0, which corresponds to (|IV.57|) . 
Because we also find V > 0, the solution is given by (IIV.58I) . Combining (|IV.44|I . (|IV.47I) . (IIV.53I) . and (jIV.58j) . 
we find the following asymptotic expression of P(4>): 



P 



us 2^ry (IV - 60) 

When (j> — > t s , P(4>) tends to vanish very rapidly. For large R, F(R) looks like (at large R) 

f Ft NO 3-1 )/ 2 * 3 

F(R) oc e (/lo/2(/3_1) HT5hFj i?-V4 x (oscillating part) . (IV.61) 

• 1 > ft > cas e: From (|IV.52p . it follows that a = j3 + 1 and, therefore, 1 - a/2 = 1/2 - ,5/2 > 0, which 
corresponds to (|IV.55|) . Because 

f-l-(-f)=/3>0, (IV.62) 

the second term in (|IV.56|) dominates when <j> -> < s if £ 7^ 0. Then, by combining (|IV.44|) . (|IV.47j) . (jIV.53|) . and 
()IV.56j) . the following asymptotic expression of P(4>) is obtained: 

P(0) ~ Se-^/ 2 ^-^^-^ 1 ^ (*. - 0) ( ^ +1)/8 . (IV.63) 

Therefore, 

- e -(^o/2(l-/3))(-6/3/ lo i?) ( ' 3 - 1)/ ^ + 1) ^7/8 _ (IV.64) 

Eq. PV.46P shows that when <f> = t — >■ t s , R — >■ 00 in the case that j3 > — 1, but i? — > in the case that /3 < — 1. 
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< case: As in 1 > /3 > case, one gets a = /3 + 1 and, therefore, 1 — a/2 > 0; however, because 



=0<O, (IV.65) 

the first term in (|IV.56|) dominates when <f> — > t s if A ^ 0. Thus, the asymptotic expression of P{4>) is as follows: 

P{4>) ~ ^-(fco/ad-fflJC*.-*) 1 -" (t s _ . (IV . 66) 

Then, F(R) is given by 

~ (-6/ lo /3 J R)( 02+2 ' 3+9 ) /8 ^ +1) e -( fe o/2(i-/3))(-6fto^)^- 1 >/^ 1 ' . (IV . 67) 

Note that ~6h (3R > when h Q , R > 0. 

When /3 > 1 in (|IV.44|1 . i? behaves as in (jIV.45|) . and when j3 < 1, the scalar curvature i? behaves as in (jiV.46[) . 
Conversely, when R behaves as 

R ~ 6H ~ i?o (/3 + 1) (* s - t) -7 , (IV.68) 
if 7 > 2, which corresponds to /3 = 7/2 > 1, H behaves as 

H~^(t s -t)^ /2 . (IV.69) 

On the other hand, if 2 > 7 > 1, which corresponds to 1 > /3 = 7 — 1>0, H is given by 

H-^^its-t)-^, (IV.70) 

and if 7 < 1, which corresponds to j3 — 7 — 1 < 0, one obtains 

H~H Q + (t s - . (IV.71) 

Here, Hq is an arbitrary constant, which is chosen to vanish in pV.44[) . Then, because H = a(t)/a(t), if 7 > 2, we 
find 

a(t) oc exp f(--l) ( * s " * } " 7/2+1 ) ' (IV - 72) 



when 2 > 7 > 1, a(t) behaves as 



and if 7 < 1, 



^'M'l 67( ^°_i) (ts-tr 7 ) - 'JV.T31 



a(t) cx exp [H t + — (t, - t) 7 ) . (IV.74) 

In any case, there appears a sudden future singularity at t = t s . 

Because the second term in (|iV.71|) is smaller than the first term, the asymptotic solution of Eq. (IIII.105[) looks as 
follows: 

P~fl.(l+j^(*.-^), (IV.75) 

with a constant Pq, which gives 

F(R) ~ F R + F x R 2f) /^ + V . (IV.76) 
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Here, Fq and F\ are constants. When > /3 > — 1, we find that 2/3/ (/3 + 1) < 0. On the other hand, when (3 < — 1, 
it follows 2/3/ (/3 + 1) > 2. As is seen from (|IV.42p . F(R) generates the Big Rip singularity when R is large. Thus, 
even if R is small, the F(R) generates a singularity where higher derivatives of H diverge. 

Let us also investigate how the effective EoS parameter w c ff for pV.71|) behaves when t ~ t s . In Einstein gravity, 
the FRW equations are given by 

\h 2 = P , ~ (2H + 3H 2 ) =p. (IV.77) 
As was discussed in section A of the second chapter, for F(R) gravity, we may introduce the EoS parameter w e e by 



Then, if ft = j/2 > 1, it follows that 



2U 

«>eff = -1 - 3jya . (IV.78) 



w e s~-l-^-(t s -t)- 1+f! (IV.79) 



when t -> t s . If 1 > /3 = 7 - 1 > 0, we find 



w cS ~-^-{t s -ty 1+f3 -+-QO. (IV.80) 
3h 

Finally, if j3 = 7 — 1 <0, one gets the expression (|IV.80[) when H = 0. When H 7^ 0, on the other hand, we obtain 



«feff~-i-^(t.-tr 1 -", (iv.M) 



when i?o vanishes. Then, if — 1 < /3 < 0, u^s — ► +00 when t —> t s . On the other hand, if /3 < — 1, u> e g — > — 1. 

Eq. (|IV.77[) also shows that, even for F(R) gravity, we may define the effective energy density p e g and the effective 
pressure p e g by 

p cff = \h 2 , PcS = -\ (2H + 3H 2 ) . (IV.82) 

We now assume that H behaves as (jIV.441) . Then, if /3 > 1, when t — > t s , a ~ exp(/io (i s — t) 1 ^ j {fi — 1)) — > 00 and 
PcS, \Pcs\ — > 00. If /3 = 1, we find a ~ (t s — i) fto — > 00 and p c ff, If off | —^00. If < j3 < 1, a goes to a constant 
but p, |p| — >• 00. If — 1 < (3 < 0, we find that a and p vanish but |p e ff| — > 00. When j3 < 0, instead of (IIV.44I) . as in 
PV.70I) . one may assume 

H ~ H + h (t s -ty P . (IV.83) 

Thus, p c ff has a finite value 3Hq/k 2 in the limit £ — > t s when — 1 < /3 < 0. If /3 < — 1 but /3 is not any integer, a is 
finite and p c g and p c ff vanish if H = or both of p e g and p e ff are finite if H =/= but higher derivatives of H diverge. 

In the classification of the singularities in (34{, Type I corresponds to the /3 > 1 or /3 = 1 case, Type II to the 
— 1 < /3 < case, Type III to the < (3 < 1 case, and Type IV to the j3 < — 1 case (where j3 is not an integer). Thus, 
we have constructed F(R) gravity examples that show all four types of the above finite-time singularities. Viable 
F(R) gravities (see second chapter) may also show (soft) future singular behavior. It is interesting to mention that 
future singularities manifest themselves in different presentations of F(R) gravity. For instance, if the scalar-tensor 
versi on sh ows the Big Rip-type singularity, then its F(R) counterpart becomes a complex theory precisely at the rip 
time (HI]. 

Near the future singularity at t = t s , the curvature becomes large in general. As a result, near the singularity, the 
quantum fields/quantum gravity effects become very important again. All classical considerations are not valid, and all 
speculations about a future cosmic doomsday cannot restrict the classical theory structure because quantum effects can 
stop (or shift) the future singularity. Moreover, the quantum corrections usually contain the curvature powers, which 
become important near the singularity. Thus, any claim about the appearance of the effective phantom/quintessence 
phase in modified gravity that subsequently enters the future singularity is not justified without quantum effects 
occurring near the singularity. One may include the massless quantum effects by taking into account the conformal 
anomaly contribution as a back-reaction near the singularity. The conformal anomaly Ta has the following well-known 
form: 

T A = b ( T + jjoR J + b'Q + b"UR , (IV.84) 
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where T is the square of the 4D Weyl tensor, and Q is the Gauss-Bonnet invariant, which are given by 



T= -R 2 -2R 11V R> J - V + R llvpu RW , Q = R 2 - 4R^R^ + R^ prT R^ pa . (IV.85) 



In general, with N scalars, Nx/2 spinors, Ni vector fields, N 2 (= or 1) gravitons and A^hd higher-derivative conformal 
scalars, 6 and 6' are given by 

N + 6N 1/2 + 12iVi + 611AT 2 -87V H d , AT + llAT 1/2 + 62AT X + 1411AT 2 - 28AT HD 
b= 120(4^ ' b= 360(4^ • (IV ' 86) 

As is seen, b > and b' < for the usual matter, except the higher-derivative conformal scalars. Notice that b" can 
be shifted by the finite renormalization of the local counterterm R 2 , so b" can be an arbitrary coefficient. 
By including the trace anomaly, the trace of (|II.5|) is modified as 

R + 2f(R) - Rf'(R) - 3Df'(R) = ~ (T mattor + T A ) . (IV.87) 
For the FRW universe, we find 

T = , g = 24 (HH 2 + H 4 ) . (IV.88) 



We now assume that H behaves as (|IV.44[) and neglect the contribution from matter by putting T mat tor = 0. Then, 
in the case that f(R) behaves as f(R) ~ — ^r, Eq. PV.87I) becomes 

r+3D {r^)=-y Ta - (iv - 89) 

First, we consider the case that 2b/3 + b" = and, therefore, Ta — G- The solution (|IV.83|) is assumed. If — 1 < /3 < 0, 
R behaves as (t s - t)" ^ 1 and Q behaves as Q - 24HH 2 - (t s - ty 30 ^ 1 . Because -3/3-1 >-/3-l when /3 < 0, T A 
is negligible compared with R. Here, □ (an/i? n+1 ) ~ (t s — ^)( ,3+1 ^ n+1 ) _2 . Then, the curvature singularity appears 
in a finite time, and the quantum correction does not prevent the singularity when 26/3 + b" = 0. One may, however, 
consider the case that 26/3 + b" 7^ 0. In this case, Ta behaves as 

T A - 06 + b") OR ~ (t a - ty 13 - 3 . (IV.90) 

Because R - (t s - t)'^ 1 and □ (an/R n+1 ) - (t s - i)(/ 3 + 1 )(«+ 1 )- 2 ; fa e terms in tne l.h.s of PV.89P are always less 
singular than Ta because 

-/3-1, (/3 + 1) (n + 1) - 2 > -/3 - 3 . (IV.91) 

This indicates that Eq. pV.89|) does not allow the singular solution and the curvature singularity does not appear. 
Therefore, in the case that 26/3 + 6" ^ 0, the quantum effects prevent the singularity appearance. 

In the above analysis, the OR term acts against the singularity. The OR term is generated by a local term R 2 , 
which shows that if one modifies F(R) by adding the R 2 term as 

F(R) F(R) + 7 i? 2 , (IV.92) 

with a constant 7, the curvature singularity cannot be generated. Thus, the addition of R 2 term which is relevant 
at the very early universe gives the universal tool for resolving the finite-time future singularities. Moreover, as was 
mentioned in section A of the second chapter, this term is the source of early-time inflation. Thus, removal of the 
singularity gives the natural prescription for the unification of the inflation with the dark energy epoch in frames of 
modified gravity. 



C. Future singularities in F(Q) gravity 

Let us show that the occurrence of future singularities is a general phenomenon in modified gravity. In this section, 
which is based on ref. (4lj , it is demonstrated that modified Gauss-Bonnet gravity can also drive the universe evolution 
to the singularity. 
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We start with the F(Q) models, in which the finite-time future singularities can occur when the form of H is taken 
as Eq. (|IV.83[) . To find such F(Q) gravities, the reconstruction method is used in a form which is different from that 
in subsection IIII Dl 

Using proper functions P(t) and Q{t) of a scalar field t which is identified with the cosmic time, the action in 
Eq. (jll.1231) can be rewritten as 



d X\J — g 



2k 2 



(R + P{t)g + Q(t))+C 

r. 



The variation with respect to t yields 



dP(t) g | dQ(t) _ Q 



(IV.93) 



(IV.94) 



dt dt 

from which one can find t = t(Q). Substituting t = t(G) into Eq. (|IV.93[) . the action in terms of F(Q) follows 

F(Q) = P(t)G + Q(t) . (IV.95) 

The scale factor is presented as 

a{t) = a exp (g(t)) , (IV.96) 

where clq is a constant and g(t) is a proper function. Using Eqs. (III.126|) . (|II.173I) . (|II.L74|) . pV.96|) . the matter 
conservation law (III.144[) and then neglecting the contribution from matter, we get the differential equation 



2± ( 9 ^) d -m)-2aHt) d -^ 



dt J 



2g i (t)^+g(t)=Q. 



(IV.97) 



With the help of Eq. (|II.I26p . Q(t) is given by 

Q(t) = -2±g*{t) 



(IV.98) 



First, the Big Rip singularity can be analyzed. If ft = 1 in Eq. pV.83|) with H = 0, H and Q are given by 



H = 

The most general solution of Eq. (|IV.97[) is 

P(t) 



h 



, G 



; (i + M 



Ah (ho - 1 

where c\ and c-i are constants. From Eq. (|IV.98j) . one gets 



(ts ~ t) (t. - t) 



1 (t — t) 3 ~ h ° 

(2t s -t)t + ci^ — '- + c 2 . 



3 — ho 



(IV.99) 



(1V.100) 



Q(t) 



6h 2 



24h% 



2h a (ho- 



Furthermore, from Eq. (IIV.94I) . we obtain 

t = 



(t s - ty 



24{hl + ht) 



1/4 



+ t s 



(1V.101) 



(IV.102) 



Solving Eq. (|IV.95|i . we find the most general form F(G), which realizes the Big Rip singularity 



F(Q) = — --— VQ + ciy i +c 2 y. (lV.f03) 
ft (i — ti ) 

This is an exact solution in the case of Eq. (|IV.99|) . In general, if, for large values of Q, F{Q) ~ aQ 1 ^ 2 , where a(^= 0) 
is a constant, the Big Rip singularity occurs for any value of ho ^ 1- In the case of ho = 1, the solution of G(H, H ■■■) 
is zero for F(Q) — aQ 1 / 2 . Note that aQ^ 1+h °^ 4 is an invariant with respect to the Big Rip solution. 
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In the case of ho = 1, it is possible to find another exact solution for P(t) 

P(t) = a{t s - tf In [ 7 (t s - t) z ] , (IV.104) 

where 7(> 0) is a positive constant and q and z are constants. Eq. (|IV.97|) is satisfied for the case of Eq. (|IV.99j) if 
q = 3 — ho = 2 (and, therefore, ho — 1) and za = —1/4. From Eq. (|IV.98|1 . we have 

Q(t) = - {ts U t) 2 ln - 1)} ■ (iv.105) 



The form of F(Q) is given by 



F(£) - ^_Vain( 7 a). (IV.106) 



In general, if, for large values of Q, F(Q) ~ a\/^hi(7C/) with a > and 7 > 0, the Big Rip singularity occurs. The 
same result is found for F(Q) ~ a^/Q\n{"fQ z + Go) with a > 0, 7 > and z > 0, where Go is a constant. 

Next, we investigate the other types of singularities. If /3 ^ 1, Eq. ()IV.83|) with H = implies that the scale factor 
a(i) behaves as 



a(t) = exp 

We consider the case in which H and Q are given by 



h (t s -t) 



1-/3 



0-1 



A solution of Eq. PV.97|) in the limit t ~^ t s is given by 



with z = — 2(3 and a = — 1/4/ig. The form of F(^) is expressed as 



Thus, if, for large values of Q, F(Q) ~ —a\fQ with a > 0, a Type I singularity could appear. 
When j3 < 1, the forms of H and 5 are given by 



H =w^w 0< ' 3<1 ' < IV - m > 



An asymptotic solution of Eq. (|IV.97[) in the limit t — ^ £ s is given by 



with z = -(1 + /?) and a = 1/2/j (1 + /3). The form of becomes 
Thus, if for large values of G, F(Q) has the form 



(IV.107) 



P(t) * w — (IV.109) 



F(G) = -12^. (IV.110) 



Pit) - (IV.112) 



(IV. 113) 



F(0)~a|0P, 7=3^3;, (IV-114) 
with a > and < 7 < 1/2, we find that < (3 < 1 and a Type III singularity could emerge. 
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If, for Q -t —oo, F(Q) has the form in Eq. (IIV.114[) with a > and — oo < 7 < 0, we find -1/3 < < and a 
Type II (sudden) singularity could occur. Moreover, if for Q — > 0~, F(Q) has the form in Eq. (|IV.114p with a < 
and 1 < 7 < 00, it is obtained — 1 < < — 1/3 and a Type II singularity could occur. 

If, for Q — !> 0~, F(G) has the form in Eq. (|IV.114|) with a > and 2/3 < 7 < 1, we obtain -00 < < -1 and a 
Type IV singularity could emerge. It is also required that 7 ^ 2n/ (3n — 1), where n is a natural number. 

We can generate all of the possible Type II singularities as shown above except in the case of — —1/3, i.e., 
H = ho/(t s — t) 1 / 3 . In this case, the form of Q is: 

Q = 2ihl0 + 2Ahf ) {t s - £) 4/3 < . (IV.115) 

To find t in terms of Q, one must consider the whole expression of Q by also taking into account the low term of 
(t s — t). It follows that 

F{Q) * ^^g (g + 8h 3 ) 1/2 + ±(g + sh^ . (iv.ne) 

As a consequence, the specific model F{Q) — o\g{g + +o'2(0 + C3) 1 / 2 , where ax, 0% and C3 are positive constants, 
can generate a Type II singularity. 

Let us end this section by considering the following realistic model, again for n > 0, 

A(g) = g a g Cxr - (1V.117) 

a 2 g n + b 2 

Because for large g, one has F^Q) ~ ai/a2^ Q and for small g, one has f±{Q) — ^1/62^", we do not find any type of 
singularities for 

\ < a < I . (IV.118) 

Eq. (|IV.117|) provides an example of a realistic model free of all possible future singularities when Eq. (|IV.118|) is 
satisfied, independently of the coefficients. Moreover, this model suggests the universal scenario for resolving future 
finite-time singularities. Adding the above model to any singular dark energy results in a combined non-singular 
model. Thus, unlike convenient dark energy, which may be singular or not, (non-singular) modified gravity may 
suggest the universal recipe for resolving the future finite-time singularity. In this respect, modified gravity seems to 
be a more fundamental theory than convenient dark energies. 



D. Late-time dynamics in Hofava-Lifshitz gravity 

Let us show that future singularities may also occur in Hofava-Lifshitz dark energy models (79j . Assume the 
behavior of the Hubble rate is the same as in (IIV.44I) . Let us now consider some simple singularity examples in the 
power-law Hofava-Lifshitz F(R) gravity and compare them with the standard case. 

In the case that 

F(R) = R + xR m , (IV.119) 
one finds the necessary conditions for the presence of the singularities: 

Type I 0>1, m>y- 

Type II m < 0, -1< < (w 12Q) 

Type III < < 1, m ^ V ' ' 

Type IV 0&Q-Z, < -1 . 

It is important to stress that the conditions (|IV.120[) are only necessary for the existence of the singularity. In the 
following, only the case that A > | is considered for simplicity. Then, one sees that a solution 

9 2(m-l)(2m-l)/i „ . „„, 

a = a t\ p = pot- 2 , } \ J 1 , po = 0, IV.121 

(2m — 1) (3A — 1) — 6(m — l)/z 
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presents a singularity of Type I for 



< m < \ 
\ < m < 1 

1 < m < 2 
m = 2 
m > 2 



-6mA+2m+3A-l 
4m 2 -12m+8 
— 6mA+2m+3A — 
-> 2 -12m+8 



6mA-2m-3A+l 
8m 3 -16m 2 + 16m-8 ' 



< yU < 

1 ^ ^ 6mA-2m-3A+; 
^ P ^ 6m-6 



4m 2 

6mA-2m-3A+l , , -6mA+2m+3A-l 
" """ ^ A* ^ 4m 2 -12m+8 



6m — 6 

-6mA+2m+3A-l ^ 6mA-2m-3A+l 



(IV.122) 



>±(9A-3), 



r ^ 4m 2 -12m+8 



6m— 6 



a singularity of Type II for 



m < -1 
m = —1 
-1 < to < ^ < 



6mA — 2m- 



4m 



<M<0, 



! -4 

M < 0, 

^ > 6mA 



-2m-3A+l 



(IV.123) 



and a singularity of Type III for 



m < 



M < ^ > 



4in 2 



6mA — 2m — 3A+1 
6m — 6 



(\ ^ m ^ 1 ft ^ n ^ -6mA+2m+3A-l 

U<m<2 U < ^ < 4m 2 -12m+8 ' 

i <T m <T 1 // <T -6mA+2m+3A-l „ > f) 

„ <. TO <^ 1 ^ ^ 4m 2 -12m+8 ^ •> u ; 



1 <T m < 9 ;/<-f) ;/> -6mA+2m+3A-l 
1 <. TO <. Z /Z < U Zi > 4m 2 -12m+8 



(IV.124) 



to = 2 
m > 2 



4m 2 -12m+8 
— 6 m* 
4m 

M < 

+3A- 
4m 2 -12m+8 

This is very different from the standard case, where we have a singularity of Type I only for m > 2, and a Type III 
for i < to < 1 and to < 0. 
For the solution 



-6mA+2m+3A — 1 < ^ < Q 



a = a i 7 , p = pot 



2m 



7 



3(l + io) ' 



Po 



we have a singularity of Type III when 



x[3/x(m- l)(2m(io + 2) 



1) -m(2m- 1) (3A- 1)] 



k 2 (to [3A - 6/x - 1) + 3(w + l)/z] 
4to 2 (-3A + 6/i + 1) - 12m(tu + l)p 
(w + 1) 2 



-3 < to < 



|(-2m-3) <w< 1, 



<to< i(l-\/l3) 0<tu< 1. 



(IV.125) 



(IV.126) 



For the solution 



a = oot 7 



2/z 



Po = 0, 



' -3A + 6it + l 

which exists only in the Hofava-Lifshitz version of F(R) gravity, we have instead a singularity of Type I for 



(IV.127) 



and of Type III for 



A < ± i(3A - 1) < p < |(3A - 1) to > 1 
A > | |(3A - 1) < p < |(3A - 1) to > 1 



A < | |(3A- 1) < p < m> 1, 
A > | < p < |(3A- 1) to > 1 . 



(IV.128) 



(IV.129) 



The results above clearly show that the presence of singularities is deeply altered in the Hofava-Lifshitz version of 
F(R) gravity. In particular, it seems that the additional parameters make it much easier to realize the singularities. 
The intervals we have presented above for the parameters can then be interpreted as constraints on this type of 
Hofava-Lifshitz F(R) gravity. Thus, we have demonstrated that modified Hofava-Lifshitz gravity can be incorporated 
into the phantom like or quintessence like accelerating cosmologies, which might lead to singularities of Type I, II, 
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or III. Of course, as in the case of conventional theory one can resolve these future singularities by adding a higher- 
derivative gravitational term, which is relevant at very early times. However, the structure of this term is different 
from the case of convenient modified gravity. For example, the addition of the R 2 term does not necessarily remove 
all types of singularities in the theory under discussion. 

In summary, in this chapter, we demonstrated that an effective phantom as well as some quintessence dark energy 
models may drive the universe evolution to future finite-time singularities. This is a universal property, regardless of 
whether dark energy is a fluid, field model or modified gravity. To resolve such future singularities, one has to add 
extra gravitational terms to the theory action. In other words, resolving the singularities invites the modification of 
gravity, even if the initial dark energy was not modified gravity. The additional gravitational term, which is usually 
relevant at very early times, also helps in the realization of the inflationary phase. In this way, the additional support 
for unification of early-time inflation with late-time acceleration in modified gravity follows. 

V. DISCUSSION 

In summary, we reviewed a number of popular modified gravities including F(R) theory, modified Gauss-Bonnet 
and scalar-Gauss-Bonnet gravities, non-minimal models, non-local gravity, modified F(R) Hofava-Lifshitz theory and 
(power-counting renormalizable) covariant gravity. General properties and different representations of such theories 
were considered. Their spatially flat FRW equations, which often have higher-derivative terms and/or non-canonical 
and non-standard couplings were derived. The accelerating cosmological solutions including the de Sitter universe and 
an attempt to obtain the unified description of inflation with dark energy were proposed. It was demonstrated that 
the qualitative possibility of such a unification is a very natural property for this large class of alternative gravities. 
This is also the case for the theories with broken Lorentz symmetry like Hofava-Lifshitz gravity. 

To describe the background evolution of modified gravity, we developed the cosmological reconstruction method 
(inverse problem) in terms of cosmological time or e-folding. Using the freedom in the choice of scalar potentials and of 
the modified term function, which depends on geometrical invariants, such as curvature and the Gauss-Bonnet term, 
we arrived to master differential equations whose solution solves the problem. It was demonstrated that eventually, 
any given cosmology may be reconstructed in this way. We explicitly considered the reconstruction in scalar-tensor 
theory, Brans-Dicke gravity, the k-essence model, F(R) theory and the Lagrange multiplier F(R) theory, modified 
Gauss-Bonnet and scalar-Gauss-Bonnet gravities, Hofava-Lifshitz gravity and non-minimal Yang-Mills theory. As 
an example of the number of universe evolution eras including ACDM one, the effective quintessence or phantom 
acceleration, the unified early-time inflation with late-time acceleration was shown to be the background solution of 
the above theories. It should be stressed that even if specific modified gravity cannot capture the entire sequence of 
the universe evolution (namely, the inflation, radiation and matter dominance, dark energy era) with good accuracy, 
it may be additionally reconstructed so that it can become a viable theory. Thus, such a partial reconstruction is 
an extremely powerful tool that may help the theory to pass the known observationally related tests. Of course, the 
price for that is some extra modifications of the action. Fortunately, such extra modifications may be attributed to 
the gravitational sector relevant at the very early universe. 

Special attention was paid to late-time dynamics of the effective quintessence/phantom dark energy of arbitrary 
nature: fluid, particle model or modified gravity. It is known that all or some of the energy conditions for such a 
(phantom or quintessence) dark energy universe are violated. As a result, the universe may end up with its evolution 
in one of four known types of finite-time singularities. The appearance of all four types of future singularities in 
coupled fluid dark energy, F(R) theory, modified Gauss-Bonnet gravity and modified F(R) Hofava-Lifshitz gravity 
was demonstrated. The universal prescription for resolving such singularities that may lead to bad phenomenological 
consequences was proposed. It consists of the addition of an extra higher-derivative gravitational term relevant in the 
very early universe and helping in the realization of the inflation. In the case of F(R) theory or standard fluid dark 
energy, this term is proportional to R 2 at large curvatures. 

The advantage of the approach proposed in this review is its very general character. Even if one of the particular 
gravity models under discussion is outside of mainstream research, its properties and background solutions will remain 
mathematically and physically correct. Moreover, the developed reconstruction scheme proposes the way to change 
the properties of any particular theory in the desirable way. Thus, we did not discuss the cosmological perturbation 
theory or structure formation in modified gravity. Indeed, the standard way to study the cosmological perturbations 
in modified gravity is based on the close analogy with general relativity. As a result, higher-order differential equations 
are approximated by second-order equations. Clearly, this is not a consistent approach. To analyze the cosmological 
perturbations one has to develop a powerful numerical technique that should be applied to higher-derivative differential 
equations. We also did not discuss the local tests of the gravity (especially, post-Newtonian regime tests) in detail. 
This is due to existence of a number of observationally related reviews [|[ from a point of view. From another point 
of view, only the small curvature regime is relevant for the purposes of the local tests. Clearly, all modified gravities 
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under consideration may be parametrized by some effective theory with smaii deviations from general relativity. In 
other words, a very large class of models with qualitatively different properties may quite easily pass the local tests. 

The next important step in the investigation of modified gravity should be related to a detailed proposal on the 
description of the sequence of the entire evolution of the universe: inflation, radiation/matter dominance, and dark 
energy. As we saw in this work, the equations of motion for such theories are extremely complicated, and including 
all of the universe matter content makes them even more involved. Presumably, such a study can be done only 
numerically. This is an extremely hard task that should give the final answer to the question: is the evolution of the 
observable universe governed by modified gravity? 
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